Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 309,381 C→T D216N (GAT→AAT)  ecpB ← putative fimbrial chaperone EcpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913309,3810CT97.9% 137.8 / NA 47D216N (GAT→AAT) ecpBputative fimbrial chaperone EcpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (26/20);  minor base G (0/1);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATG  >  NC_000913/309156‑309619
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
ataGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGCATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGCACTTATTACCCTGCCAAAGTGCCACCCGTCCt                                                                                                                                                                                                                       >  1:16285/1‑251 (MQ=255)
        ggACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTTCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTATTATc                                                                                                                                                                                                               >  1:463905/1‑251 (MQ=255)
         gACGCTCATGCCGGCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTAGTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATcc                                                                                                                                                                                                              <  2:52060/251‑1 (MQ=255)
                aTGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAgtgt                                                                                                                                                                                                       <  1:205331/251‑1 (MQ=255)
                  gccttCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAAGCACCCCACCTAGTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCATACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGGCCTTTGTTATCCGCAGTGTc                                                                                                                                                                                                      <  2:272122/246‑1 (MQ=255)
                    cGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATGACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCAc                                                                                                                                                                                                   <  2:64905/251‑1 (MQ=255)
                            aGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTa                                                                                                                                                                                            >  1:115256/1‑250 (MQ=255)
                                   gactgaAGGCCTGCGGAATAATCTCCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGGCCACGCGGGTAAAacgac                                                                                                                                                                                     <  2:177626/250‑1 (MQ=255)
                                           gtcggCGGAATAAGCCCCCCACCTATTCGCTGCGCCAGTATGTCAGCAGCGCTGGCCTCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTCTCAGTCCTGGGGAGAACCGTCGTAAAGGCATGGTCGTTATCCATCTGCACAGCTATCTATTTCACGTGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGGGTCCACGCGGGTAAAACGACGCGACTTGc                                                                                                                                                                            <  2:118538/247‑1 (MQ=255)
                                             cTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGccc                                                                                                                                                                          >  1:232420/1‑251 (MQ=255)
                                                    ataatCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACGTGCCCGGCATCa                                                                                                                                                                   <  1:353891/251‑1 (MQ=255)
                                                      aatCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAgg                                                                                                                                                                 >  2:68819/1‑251 (MQ=255)
                                                                 cTATTTGCTGAGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAAAGACAAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATTAACTTATTACCATGCCAAAGTGCCACCCGCCCTTTGTTAGCCGCAGTGTCCACGCGGGTAAAACGACGCAACTTGCCCGTCATCAGGTAGTAATtcac                                                                                                                                                      >  2:110675/1‑249 (MQ=255)
                                                                  tATTTGATGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAACCTGGGCTTTAAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTCTTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTctct                                                                                                                                                     >  2:325593/1‑251 (MQ=255)
                                                                          gcgcACTTATGTCAGCAGCGCTGGCATCGGGCTTAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTacac                                                                                                                                               <  2:257803/249‑1 (MQ=255)
                                                                                  aTGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCg                                                                                                                                      <  1:42703/250‑1 (MQ=255)
                                                                                    ggcagcagCGCTGGCATCGGGATGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGGTATCCATCTGCACAGCTATCTATTTCACGGGAAGGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTg                                                                                                                                   <  1:315282/249‑1 (MQ=255)
                                                                                          agcGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTAAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCt                                                                                                                             >  1:482786/1‑251 (MQ=255)
                                                                                                   aTCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGAGGCTTTCAGGCa                                                                                                                     >  1:431838/1‑250 (MQ=255)
                                                                                                    tCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCa                                                                                                                     >  1:310365/1‑249 (MQ=255)
                                                                                                                    aaaaGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGAt                                                                                                   >  1:454558/1‑251 (MQ=255)
                                                                                                                         aCCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGAGCCGCa                                                                                              <  2:361877/251‑1 (MQ=255)
                                                                                                                            aTGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCg                                                                                           >  2:458524/1‑251 (MQ=255)
                                                                                                                               ccGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCg                                                                                        <  1:166841/251‑1 (MQ=255)
                                                                                                                                 gAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGTGAATGAACGTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTcc                                                                                                                                                                                                     <  2:297651/140‑1 (MQ=255)
                                                                                                                                 gAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAACGTGCAACCCGTCCTTTGTTATCCGCAGAGTcc                                                                                                                                                                                                     >  1:297651/1‑140 (MQ=255)
                                                                                                                                      cAAACTGGGGCTGCAGTCCGGGGGCGAGCCCTCGTAAAGGCATGGTCGTTATCCATCTGCACAGCTATCTATTGCCCGGGAATGAAATTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGGGTCCACGCGGGTAAAACGACGCGACGTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGGCCGTCGGCGAGGATCCGCAGCGACGCATTTcc                                                                                 <  2:84459/251‑1 (MQ=255)
                                                                                                                                          cTGGGCTTTCAGTCCTGGGGCGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTACCCGTTggc                                                                                                                                 >  2:313922/1‑199 (MQ=255)
                                                                                                                                          cTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTggc                                                                                                                                 <  1:313922/199‑1 (MQ=255)
                                                                                                                                           tGGGCTTTCAGTCCTGGGGAGAACCGTCGTAACGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGCCGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGTTTTCAGGCCAGGCCCGTAGGCGAGGATCCGCAGCGTCGCCATTCCTgcat                                                                            >  2:393011/1‑248 (MQ=255)
                                                                                                                                            gggCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTAtt                                                                           <  1:486210/251‑1 (MQ=255)
                                                                                                                                             gacTTTCGGTCCTGGGGAGAACCGTCGTCAAGGCATTGTCGTTATCCATCTGCACAGATATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTTAAACGACGCGGCTTGCCCGGCATCGGGAAGTAATTCTCTTTACGCTCCTTACGGCTGGCGGCTGTCAGGGTGGGTCCGTAGGCGAGGCTCCGCAGCGTCGCCTTTCCTGTAttt                                                                          >  2:199749/3‑251 (MQ=255)
                                                                                                                                                            ggAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTg                                                           <  1:338834/251‑1 (MQ=255)
                                                                                                                                                              agaACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGc                                                         <  1:35516/251‑1 (MQ=255)
                                                                                                                                                                            ggCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGt                                           >  1:91752/1‑251 (MQ=255)
                                                                                                                                                                             gCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAACGTGGt                                           >  1:292368/1‑250 (MQ=255)
                                                                                                                                                                             gCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGt                                           >  1:207234/1‑250 (MQ=255)
                                                                                                                                                                             gCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGt                                           >  2:171479/1‑250 (MQ=255)
                                                                                                                                                                               aTTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGATCTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGTCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGt                                        <  2:199332/251‑1 (MQ=255)
                                                                                                                                                                                  gTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGTTAGTCATTCTCTTTACCCTCCTTACCGTTGGCGGCTTTCAGGCAAGGACCGTAGGCGAGGATCCGCAGCATCGCATTTCCTGTATTTGGCAGGGAGCCGTTGGCGGAATGAAAGTGGTAGTTCg                                     >  2:414430/1‑251 (MQ=255)
                                                                                                                                                                                      ttATCCATCTGCACAGCTATCTATTTCAAGGGAATGAACTTATTCCCCTGCCCAAGTGCAACCCGTCCTTTGCTATAAGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGCAGTAATTCTCTTTCCACTCATTACAGTTGACGACTTTCAGGCAAGGTCCGTAGGCGCGGATCCGCAGCGTCGCCTTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGCAAGTGGAAGTTCGCCTg                                 >  1:29285/1‑251 (MQ=255)
                                                                                                                                                                                                    aGCTATCTATTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGcgacg                    >  1:178546/1‑250 (MQ=255)
                                                                                                                                                                                                            atTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGCCGACCAg               <  2:78730/247‑1 (MQ=255)
                                                                                                                                                                                                            atTTCACGGGAATGAACTTATTACCCTGCCAAAGTGCAACCCGGCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAAt            <  2:310365/250‑1 (MQ=255)
                                                                                                                                                                                                                   gggAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGCCTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCTAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCCCCTGACGGGGGGCAACGCCCGGAATGGTGCCg     >  2:476488/1‑250 (MQ=255)
                                                                                                                                                                                                                    ggAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGAt   >  2:376377/1‑251 (MQ=255)
                                                                                                                                                                                                                     gAATGAACTTATTACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGGGGTAGTTCGCCTGACGAGGGGCGGCGACCAGAATGGTGCCgttg  >  2:170242/1‑248 (MQ=255)
                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                              
ATAGAGCGGGACGCTCATGCCGTCCTGAAGCGCCTGACTGAAGGCCTGCGGAATAATCACCCCACCTATTTGCTGCGCACTTATGTCAGCAGCGCTGGCATCGGGCTGAACGAACAAAAAGACCATGCCGAAGGCAAACTGGGCTTTCAGTCCTGGGGAGAACCGTCGTAAAGGCATTGTCGTTATCCATCTGCACAGCTATCTATTTCACGGGAATGAACTTATCACCCTGCCAAAGTGCAACCCGTCCTTTGTTATCCGCAGTGTCCACGCGGGTAAAACGACGCGACTTGCCCGGCATCAGGTAGTAATTCTCTTTACACTCCTTACCGTTGGCGGCTTTCAGGCAAGGTCCGTAGGCGAGGATCCGCAGCGTCGCATTTCCTGTATTTGTCAGGGAGCCGTTGGCGTACTGAAAGTGGTAGTTCGCCTGACGAGGGGCGACGACCAGAATGGTGCCGATG  >  NC_000913/309156‑309619

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