Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,366 T→A 9.6% intergenic (+23/+24) slyB → / ← slyA outer membrane lipoprotein/global transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,720,3660TA 9.6% 212.7 / 4.1 86intergenic (+23/+24)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (57/20);  new base A (3/6);  total (60/26)
Fisher's exact test for biased strand distribution p-value = 1.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTC  >  NC_000913/1720217‑1720515
                                                                                                                                                     |                                                                                                                                                     
ggAAATTCGTAAAGACGATGGTAGTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGtt                                                                                                                                                       >  2:677401/1‑150 (MQ=255)
    aTTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACAGTTTCTCCGCGATAAATAACTTTGAGTGTGGTCAGGTTACTg                                                                                                                                                   >  1:106702/1‑150 (MQ=255)
        gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                                                                                                                                               <  2:1646722/150‑1 (MQ=255)
        gTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACca                                                                                                                                               >  1:607465/1‑150 (MQ=255)
         tAAAGACGATGGTAATACCATAATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcac                                                                                                                                              <  2:1443984/150‑1 (MQ=255)
          aaaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTAACCTGACcaca                                                                                                                                             <  1:1787266/150‑1 (MQ=255)
          aaaGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACcaca                                                                                                                                             >  2:1847/1‑150 (MQ=255)
            ccatgATGGTATTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACg                                                                                                                                           >  2:669827/5‑150 (MQ=255)
             gACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGc                                                                                                                                          >  2:148969/1‑150 (MQ=255)
                gATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGTAACCTGACCACACGcaaa                                                                                                                                       <  1:1205161/150‑4 (MQ=255)
                   ggTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCtt                                                                                                                                    >  2:977569/1‑150 (MQ=255)
                   ggTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGACGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCtg                                                                                                                                    >  2:730114/1‑149 (MQ=255)
                     tAATACCATCATGGTGGTACAGAAACAAGGCATCACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACAACACTCCCCCTTCa                                                                                                                                  >  1:1143113/1‑150 (MQ=255)
                     tAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCa                                                                                                                                  >  1:218340/1‑150 (MQ=255)
                      aaTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCAt                                                                                                                                 >  2:1688235/1‑150 (MQ=255)
                       aTACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGCCCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGACCCCTTCAtt                                                                                                                                >  2:1626294/1‑150 (MQ=255)
                        tACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCCCCCGCCCCCTGCAttt                                                                                                                               >  1:489809/1‑150 (MQ=255)
                          ccATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGCCCACACGCCCCCTTCATTTCa                                                                                                                             >  2:709932/1‑150 (MQ=255)
                               atGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACAAGCCCCCTTCATTTCACCCtt                                                                                                                        >  1:1653875/1‑150 (MQ=255)
                                  gtggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGGGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTgg                                                                                                                     <  1:1080735/150‑1 (MQ=255)
                                    ggtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGcc                                                                                                                   >  1:789703/1‑150 (MQ=255)
                                     gtACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGGCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCt                                                                                                                  <  1:110288/150‑1 (MQ=255)
                                        cAGAAAAAAGGCAACACTCTTTTCTCTCCGGGCCAACGGGTGGTAGTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGGGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTa                                                                                                               <  2:745205/150‑1 (MQ=255)
                                         aGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACCCGCCCCCTTCATTTCACCCTTTGGCCGGtta                                                                                                              >  1:74639/1‑148 (MQ=255)
                                          gAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAAc                                                                                                             >  1:578142/1‑150 (MQ=255)
                                          gAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTAGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGCCCTGTAAc                                                                                                             >  2:33470/1‑150 (MQ=255)
                                          gAAACAAGGCAACACTCGTTTCTCTCCGGGCAAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTATCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAAc                                                                                                             >  1:586271/1‑150 (MQ=255)
                                                aGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGGGTGGTCTGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATg                                                                                                       <  2:294289/150‑1 (MQ=255)
                                                 ggCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGa                                                                                                      >  1:984509/1‑150 (MQ=255)
                                                 ggCAACACTCGTTTCTCTCCGGGCCAACGTGTCATACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGa                                                                                                      >  2:49721/1‑150 (MQ=255)
                                                        cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATTATATTATg                                                                                               >  2:1655196/1‑150 (MQ=255)
                                                        cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                               >  2:1331063/1‑150 (MQ=255)
                                                        cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                               >  2:1755781/1‑150 (MQ=255)
                                                        cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                               >  1:489995/1‑150 (MQ=255)
                                                         tCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCGTCATTTCAACCTTTGGCCTGTTACTCAATGATATTATGc                                                                                              >  2:715893/1‑150 (MQ=255)
                                                          cGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCt                                                                                             >  2:313516/1‑150 (MQ=255)
                                                              tctctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATTTTATGCTCAAg                                                                                         >  1:298027/1‑150 (MQ=255)
                                                              tctctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAg                                                                                         >  1:493002/1‑150 (MQ=255)
                                                              tctctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGCTGACCGATTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATGTCACCCTTTGGCCTGTAACTCAAAGATATTATGCTCAAg                                                                                         >  2:506485/1‑150 (MQ=255)
                                                               ctctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGt                                                                                        >  2:961064/1‑150 (MQ=255)
                                                                 ctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTACCTCAATGATATTATGCTCAAGttt                                                                                      >  1:1030717/1‑150 (MQ=255)
                                                                 ctcCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGttt                                                                                      >  1:1675289/1‑150 (MQ=255)
                                                                   ccGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGATAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTg                                                                                    <  2:218340/150‑1 (MQ=255)
                                                                   ccGGCCCAACGTGCCTTAGTGGCCACCAATGGCAGTCAGGTGACCGTTTCACCGCGCTAAATATGTTTGCGTGTGGTCTGTTTACTGACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTg                                                                                    <  1:134826/150‑1 (MQ=255)
                                                                               gtCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTaa                                                                        >  2:829333/1‑150 (MQ=255)
                                                                               gtCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTaa                                                                        >  2:1696269/1‑150 (MQ=255)
                                                                                    aCTGGCCAGCATGGGCTGTCGGGTTACGGTTTCTCCGCTAAAAATAAGTTTGGGTGTGGTCTGGTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGt                                                                   <  1:677401/150‑1 (MQ=255)
                                                                                     cTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGATAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGtt                                                                  <  1:1018978/150‑1 (MQ=255)
                                                                                       ggCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGc                                                                >  1:1444566/1‑150 (MQ=255)
                                                                                          cAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                             >  2:1595086/1‑150 (MQ=255)
                                                                                          cAGAAATGGCAGTCAGGTGCCCGTTTCTCCGCGCAAAAAAAGTTTGCGGGTGGTCAGGTAACTCACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTcc                                                             <  2:1804752/150‑1 (MQ=255)
                                                                                            gCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAg                                                           >  1:1803897/1‑150 (MQ=255)
                                                                                            gCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAg                                                           >  2:833356/1‑150 (MQ=255)
                                                                                             taaTGGCAGTCAGGTGACGGTTTCTCCGCGTTAAATTAGTTTGCGTGTGGTCAGGTATCTGAGCAGACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGt                                                          <  1:73586/149‑1 (MQ=255)
                                                                                              attGGCTGTCTGGTGACGTTTTCTCGGCGCTAAATATGTTTGGGTGGGGTCTGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGtt                                                         <  2:649643/148‑1 (MQ=255)
                                                                                              aaTGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGtt                                                         >  2:710405/1‑150 (MQ=255)
                                                                                                tgccAGTCTGGTGACCGTTACGCCGCGCAAAATATGTTTCCGTGGGTTCAGGTAACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTcc                                                       <  2:149881/147‑1 (MQ=255)
                                                                                                tGGCAGGCTGGTGACCGTTTCTCCGCGCTAAAAAAGTTTGCGTGGGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTcc                                                       <  2:1227633/150‑1 (MQ=255)
                                                                                                   cAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCctct                                                    >  2:1661829/1‑150 (MQ=255)
                                                                                                   cAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCctct                                                    >  2:365531/1‑150 (MQ=255)
                                                                                                    aGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTg                                                   >  1:1396634/1‑150 (MQ=255)
                                                                                                     gTCTGGTGACCTTTTCACCGCTATAAATATTTTGGGGTGGGGTCTGGTTACTCCCCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGc                                                  <  2:1418446/150‑1 (MQ=255)
                                                                                                      tctggTTACCTTTTCTCCGCGATAAATACGTTTGGGTGTGGTCCGTTAACTTACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCg                                                 <  1:148969/147‑1 (MQ=255)
                                                                                                              ccGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGcc                                         >  1:1734387/1‑150 (MQ=255)
                                                                                                               cGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGAAGCCa                                        >  2:1263714/1‑150 (MQ=255)
                                                                                                                gTTTCTCCGCGCTAAATAAGTTTGCGTGTTGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTAAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGAGGAGATGCCAt                                       >  2:845737/1‑150 (MQ=255)
                                                                                                                      ccGCGTAAAATATTTTTGTGTGTGGTAAGTTTACTCACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAata                                 <  2:710212/150‑1 (MQ=255)
                                                                                                                      ccGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCCCCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAata                                 >  2:22744/1‑150 (MQ=255)
                                                                                                                       cgcgCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAatat                                >  2:408293/1‑150 (MQ=255)
                                                                                                                          gcTAAAAAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTc                             <  1:313516/150‑1 (MQ=255)
                                                                                                                              aaTAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCcgcg                         >  1:666700/1‑150 (MQ=255)
                                                                                                                                    tttGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCCCCCTTTGGCCTGTAACTCAATGATATTATGGTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGtt                   >  2:24005/1‑150 (MQ=255)
                                                                                                                                      tGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTACGCGCGGGtttt                 >  2:900702/1‑150 (MQ=255)
                                                                                                                                       gCGTGTGGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg                >  1:1437873/1‑150 (MQ=255)
                                                                                                                                       gCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg                >  2:1221367/1‑150 (MQ=255)
                                                                                                                                       gCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg                >  2:1769620/1‑150 (MQ=255)
                                                                                                                                        cGTGTGGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGt               >  1:1276955/1‑150 (MQ=255)
                                                                                                                                         gtgtGGTCAGGTTACTTACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGtt              <  2:205351/150‑1 (MQ=255)
                                                                                                                                         gtgtGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGtt              <  2:194612/150‑1 (MQ=255)
                                                                                                                                          tgtgGTCAGTAACCTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTta             >  2:1534825/1‑150 (MQ=255)
                                                                                                                                          tggggTAAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCGCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTta             <  1:1847/147‑1 (MQ=255)
                                                                                                                                             ggTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaata          >  2:1569792/1‑150 (MQ=255)
                                                                                                                                             ggTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaata          >  2:923711/1‑150 (MQ=255)
                                                                                                                                             ggTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaata          <  1:1199152/150‑1 (MQ=255)
                                                                                                                                             ggTCAGGTTACTGACCACACGCCCCCTTCATTTCACCATTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaata          >  2:792732/1‑150 (MQ=255)
                                                                                                                                              gTCTGGTTACTTACAACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTtaataa         <  2:578142/150‑1 (MQ=255)
                                                                                                                                                  ggTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCGCTGCGGAGATGCCATGTAATATTTCCGCGCGGGGTTTGTTAATAACAGc     >  2:955741/1‑150 (MQ=255)
                                                                                                                                                     tACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTc  >  2:1007275/1‑150 (MQ=255)
                                                                                                                                                     tACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTc  >  2:1140189/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGAAATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTC  >  NC_000913/1720217‑1720515

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