Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,720,367 A→C 5.7% intergenic (+24/+23) 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,3670AC 5.7% 229.7 / 6.2 88intergenic (+24/+23)slyB/slyAouter membrane lipoprotein/global transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (59/23);  new base C (3/2);  total (62/26)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

ATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTC  >  NC_000913/1720221‑1720515
                                                                                                                                                  |                                                                                                                                                    
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                                                                                               >  1:489995/1‑150 (MQ=255)
                                                    cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                               >  2:1755781/1‑150 (MQ=255)
                                                    cTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATg                                                                                               >  2:1331063/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:1769620/1‑150 (MQ=255)
                                                                                                                                   gCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTg                >  2:1221367/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)
                                                                                                                                                  |                                                                                                                                                    
ATTCGTAAAGACGATGGTAATACCATCATGGTGGTACAGAAACAAGGCAACACTCGTTTCTCTCCGGGCCAACGTGTCGTACTGGCCAGCAATGGCAGTCAGGTGACCGTTTCTCCGCGCTAAATAAGTTTGCGTGTGGTCAGGTTACTGACCACACGCCCCCTTCATTTCACCCTTTGGCCTGTAACTCAATGATATTATGCTCAAGTTTTGCGATGAGCGTAATCAGTTGCTCCAGTTCCTCTGCGGAGATGCCATGTAATATTTCCGCGCGGGTTTTGTTAATAACAGCTTC  >  NC_000913/1720221‑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: