Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,483,330 C→T D93D (GAC→GAT)  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,3300CT100.0% 178.1 / NA 65D93D (GAC→GAT) unknownType II/IV secretion system ATP hydrolase TadA/VirB11/CpaF, TadA subfamily
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/35);  total (30/35)

GTACTAAGTAATCTTGATATTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCACTGGTAAACGAACTGGTTCAGATTGTTAACCGTGTTGCCAATAAAAGCGGTGCTTATTTGACCTCAGCCGATACTCTAGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGACCATTACGAGATCTCATGGCTGACGATACTATAAATGATATTTTAGTTAATGGCCCTAATGACATCTGGGTTGAACGTGCTGGTATTCTGGAAAAAACTGATAAAGAATTTGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGCCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCAGTAAGAATAAAAAAACATTGCAAG  >  1/1483091‑1483574
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
gtattgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGa                                                                                                                                                                                                                                           >  1:590920/10‑251 (MQ=255)
gtattgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGa                                                                                                                                                                                                                                           >  1:676058/10‑251 (MQ=255)
  attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGGCGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                                                                                                                                                                                         <  1:36229/244‑1 (MQ=255)
  attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                                                                                                                                                                                         <  1:177570/244‑1 (MQ=255)
  attgagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAAc                                                                                                                                                                                                                                         <  1:860441/244‑1 (MQ=255)
     gagcaaTCTTGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGt                                                                                                                                                                                                                                              >  1:188727/5‑243 (MQ=255)
              tGATGTTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCATTAGTAAACGAACTGGTTCAGATTGTTAATCGTATTGCCAATAAAAGTAACACTTATTTGACGTCATCCGATGCACTGGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGTCCGTTACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTc                                                                                                                                                                                                                             <  1:211039/251‑1 (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)
                                                                             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:823350/251‑1 (MQ=255)
                                                                                                                                                                                     tACGGGATCTCGTGGCTGACGATACTATAAATGATATCTTAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGAt                                                      <  1:705393/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)
                                                                                                                                                                                                                        atcttAGTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGag                   >  1:258190/4‑250 (MQ=255)
                                                                                                                                                                                                                              gTTAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTc                           >  1:649862/1‑237 (MQ=255)
                                                                                                                                                                                                                               ttAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTaaaaaa            <  1:286497/251‑1 (MQ=255)
                                                                                                                                                                                                                               ttAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTaaaaaa            <  1:834553/251‑1 (MQ=255)
                                                                                                                                                                                                                               ttAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTaaaaaa            >  1:247948/1‑251 (MQ=255)
                                                                                                                                                                                                                               ttAATGGCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTaaaaa             <  1:170150/250‑1 (MQ=255)
                                                                                                                                                                                                                                     gCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAAtgttgtaat                                                                         <  1:644040/184‑4 (MQ=255)
                                                                                                                                                                                                                                     gCCCTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATggtac                                                 <  1:6628/208‑4 (MQ=255)
                                                                                                                                                                                                                                      cccTAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTaaaaaaagcttg      >  1:813685/1‑245 (MQ=255)
                                                                                                                                                                                                                                         tAATGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCATTTCCATTCGTAAGTTCAGTAAGAGTAAAAAAAGCTTGCAag  <  1:136815/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaTGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAg                                                                                                                    >  1:520387/1‑136 (MQ=255)
                                                                                                                                                                                                                                          aaTGATATCTGGGTTGAACGTGCCGGTATTCTGGAAAAAACTGATAAAGAATTCGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGTCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCCGATGGTAGCCGTTTAAATGTTGTAATTGCACCGATTGCGTTAGATGGTACATCCa                                            >  1:534121/1‑208 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
GTACTAAGTAATCTTGATATTGAAAAGATTGATGAACTTCAAAGTGAAAGAAGCTCACTGGTAAACGAACTGGTTCAGATTGTTAACCGTGTTGCCAATAAAAGCGGTGCTTATTTGACCTCAGCCGATACTCTAGTTATGGCTGAAATTGTTGCTGATGAGATTGAGGGTTATGGACCATTACGAGATCTCATGGCTGACGATACTATAAATGATATTTTAGTTAATGGCCCTAATGACATCTGGGTTGAACGTGCTGGTATTCTGGAAAAAACTGATAAAGAATTTGTTAGTAATGAACAATTAACTGATATTGCAAAAAGGTTAGTTGCTCGGGTTGGTCGCCGTATTGATGATGGTAGTCCACTAGTGGACTCCCGACTACCTGATGGCAGCCGTTTAAATGTTGTTATTGCACCGATTGCATTAGATGGCACATCCATTTCTATTCGTAAATTCAGTAAGAATAAAAAAACATTGCAAG  >  1/1483091‑1483574

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