Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,483,309 T→C I86I (ATT→ATC)  unknown → Type II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,483,3090TC100.0% 167.7 / NA 59I86I (ATT→ATC) unknownType II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/33);  total (26/33)

AACAAAAAATCTTACTTAGAAGTGAAGTACTAAGTAATCTTGATATTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCACTGGTAAACGAACTGGTTCAGATTGTTAACCGTGTTGCCAATAAAAGCGGTGCTTATTTGACCTCAGCCGATACTCTAGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGACCATTACGAGATCTCATGGCTGACGATACTATAAATGATATTTTAGTTAATGGCCCTAATGACATCTGGGTTGAACGTGCTGGTATTCTGGAAAAAACTGATAAAGAATTTGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGCCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCA  >  1/1483065‑1483550
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
aacaaAAAATCTTACTTAGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGtt                                                                                                                                                                                                                                             <  1:536723/251‑1 (MQ=255)
              cttaGAAGTGAAGTATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGa                                                                                                                                                                                                                               <  1:840605/251‑1 (MQ=255)
                     gtgaagtATTGAGCAATCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGa                                                                                                                                                                                                                               <  1:152822/244‑1 (MQ=255)
                          gtattgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGa                                                                                                                                                                                                                   >  1:590920/10‑251 (MQ=255)
                          gtattgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGa                                                                                                                                                                                                                   >  1:676058/10‑251 (MQ=255)
                            attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGGCGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                                                                                                                                                                 <  1:36229/244‑1 (MQ=255)
                            attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                                                                                                                                                                 <  1:860441/244‑1 (MQ=255)
                            attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                                                                                                                                                                 <  1:177570/244‑1 (MQ=255)
                               gagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGt                                                                                                                                                                                                                      >  1:188727/5‑243 (MQ=255)
                                    aTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATg                                                                                                                                                                                                                                >  1:816649/1‑228 (MQ=255)
                                        tGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTc                                                                                                                                                                                                     <  1:211039/251‑1 (MQ=255)
                                                                   aaGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGtt                                                                                                                                                                                                                                             >  1:308875/1‑184 (MQ=255)
                                                                     gTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATc                                                                                                                                                                                                                           >  1:264397/1‑200 (MQ=255)
                                                                     gTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTc                                                                                                                                                                                                     <  1:23339/222‑1 (MQ=255)
                                                                              gctcattagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTg                                                                                                                                                                                                                         <  1:645237/185‑1 (MQ=255)
                                                                              gctcattagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAAc                                                                                                                                                                 >  1:378392/9‑249 (MQ=255)
                                                                                   ttagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGt                                                                                                                                                                                                         <  1:473317/201‑1 (MQ=255)
                                                                                   ttagTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTaa                                                                                                                                                           >  1:719901/4‑250 (MQ=255)
                                                                                       tAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAGTTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACAGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGa                                                                                                                                                       <  1:91403/250‑1 (MQ=255)
                                                                                       tAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGat                                                                                                                                                      >  1:203429/1‑251 (MQ=255)
                                                                                             aaCTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCa                                                                                                                                                >  1:865880/1‑251 (MQ=255)
                                                                                                       gATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGc                                                                                                                                                                                                                                       <  1:547307/154‑1 (MQ=255)
                                                                                                       gATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACt                                                                                                                                                         <  1:338096/232‑1 (MQ=255)
                                                                                                                          aataaaagtaacacTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCgtcg                                                                                                                   >  1:15157/14‑248 (MQ=255)
                                                                                                                                         tATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATtgatg                                                                                                           >  1:835611/1‑244 (MQ=255)
                                                                                                                                         tATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATtgatg                                                                                                           >  1:287448/1‑244 (MQ=255)
                                                                                                                                           tttGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTc                                                                                                  <  1:865849/251‑1 (MQ=255)
                                                                                                                                                           gcactggTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATc                                                                                                                                                                                                                           <  1:415348/108‑1 (MQ=255)
                                                                                                                                                                ggTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATc                                                                                                                                                                                                                           <  1:1732/108‑1 (MQ=255)
                                                                                                                                                                ggTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGAc                                                                                       <  1:845249/240‑1 (MQ=255)
                                                                                                                                                                ggTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTa                                                                              >  1:521936/2‑250 (MQ=255)
                                                                                                                                                                  ttATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTAccc                                                                           <  1:758068/251‑2 (MQ=255)
                                                                                                                                                                           aaaTTGTTGCTGATGAGATTGAGGGTTATAGTCCGTTTCGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTAt                                                                                                                                                                                                       <  1:184739/118‑1 (MQ=255)
                                                                                                                                                                            aaTTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATc                                                                                                                                                                                                                           <  1:54687/97‑1 (MQ=255)
                                                                                                                                                                              ttgttgCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATtgatg                                                                                                           >  1:353882/1‑207 (MQ=255)
                                                                                                                                                                                           gaTTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATtgatg                                                                                                           >  1:742115/1‑194 (MQ=255)
                                                                                                                                                                                            aTTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAAtgttgtaat                                                 <  1:188624/251‑4 (MQ=255)
                                                                                                                                                                                                 gggTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGc                                                                  <  1:594915/229‑1 (MQ=255)
                                                                                                                                                                                                 gggTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAAt                                                         <  1:329677/238‑1 (MQ=255)
                                                                                                                                                                                                  ggTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCAc                                            >  1:698670/1‑250 (MQ=255)
                                                                                                                                                                                                   gTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCg                                          >  1:698551/1‑251 (MQ=255)
                                                                                                                                                                                                         gtccgttACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGcg                                    >  1:41424/6‑250 (MQ=255)
                                                                                                                                                                                                               tACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGAt                              <  1:705393/251‑1 (MQ=255)
                                                                                                                                                                                                               tACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGAt                              <  1:823350/251‑1 (MQ=255)
                                                                                                                                                                                                                     aTCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGAc                                                                                       <  1:630038/188‑1 (MQ=255)
                                                                                                                                                                                                                         cgtGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCa                    <  1:863456/249‑1 (MQ=255)
                                                                                                                                                                                                                             gCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACt                                                                               <  1:42072/188‑1 (MQ=255)
                                                                                                                                                                                                                                  cGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACt                                                                                                                                                         <  1:754147/109‑1 (MQ=255)
                                                                                                                                                                                                                                  cGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGAt                              <  1:780385/232‑1 (MQ=255)
                                                                                                                                                                                                                                    aTACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCAc                                                                                               >  1:113331/1‑165 (MQ=255)
                                                                                                                                                                                                                                     tACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTa        <  1:834489/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATtg                                                                                                              <  1:453546/147‑1 (MQ=255)
                                                                                                                                                                                                                                       cTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTaag      >  1:852340/1‑250 (MQ=255)
                                                                                                                                                                                                                                       cTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTaag      >  1:366335/1‑250 (MQ=255)
                                                                                                                                                                                                                                         ataAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAg                                                                                            >  1:63955/1‑163 (MQ=255)
                                                                                                                                                                                                                                           aaaTGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCa  >  1:182602/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aaTGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGAc                                                                                       <  1:832628/165‑1 (MQ=255)
                                                                                                                                                                                                                                               gATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGTTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAAtg                                                        >  1:37101/1‑193 (MQ=255)
                                                                                                                                                                                                                                                atatCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTaagt     >  1:46413/1‑241 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 
AACAAAAAATCTTACTTAGAAGTGAAGTACTAAGTAATCTTGATATTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCACTGGTAAACGAACTGGTTCAGATTGTTAACCGTGTTGCCAATAAAAGCGGTGCTTATTTGACCTCAGCCGATACTCTAGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGACCATTACGAGATCTCATGGCTGACGATACTATAAATGATATTTTAGTTAATGGCCCTAATGACATCTGGGTTGAACGTGCTGGTATTCTGGAAAAAACTGATAAAGAATTTGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGCCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCA  >  1/1483065‑1483550

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑