Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 212,726 T→G 5.2% S132R (AGT→AGG)  tilS → tRNA(Ile)‑lysidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913212,7260TG 5.2% 292.9 / 6.9 96S132R (AGT→AGG) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base T (51/40);  new base G (3/2);  total (54/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CGCAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAG  >  NC_000913/212578‑212873
                                                                                                                                                    |                                                                                                                                                   
cgcAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTg                                                                                                                                                    >  1:1279068/1‑150 (MQ=255)
  cAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGc                                                                                                                                                  >  2:728080/1‑150 (MQ=255)
      aagGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTg                                                                                                                                              <  2:386709/150‑1 (MQ=255)
       agGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGc                                                                                                                                             >  1:19485/1‑150 (MQ=255)
          aCTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGGCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCgg                                                                                                                                          <  1:465891/150‑1 (MQ=255)
            tGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGc                                                                                                                                        <  1:353305/150‑1 (MQ=255)
            tGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGc                                                                                                                                        <  2:19485/150‑1 (MQ=255)
             gggCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTTGCGCTAAAACGCGGCAGGGCCACTGCCgcgtt                                                                                                                                       <  2:1126172/150‑5 (MQ=255)
             gggCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCt                                                                                                                                       >  1:486923/1‑150 (MQ=255)
                cATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTc                                                                                                                                    >  2:31480/1‑150 (MQ=255)
                cATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTc                                                                                                                                    <  1:1664562/150‑1 (MQ=255)
                 aTTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCg                                                                                                                                   <  2:82935/150‑1 (MQ=255)
                 aTTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGGGCCACTGCCGCGttttag                                                                                                                                   <  2:803337/150‑3 (MQ=255)
                   tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                 <  2:254123/150‑1 (MQ=255)
                   tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                 <  2:1016673/150‑1 (MQ=255)
                   tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                 <  2:1412790/150‑1 (MQ=255)
                    gAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGTCCTGCCGGGCTTTCGGCt                                                                                                                                >  2:628892/1‑150 (MQ=255)
                    gAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCt                                                                                                                                >  2:1231205/1‑150 (MQ=255)
                    gAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTCGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCt                                                                                                                                >  1:301334/1‑150 (MQ=255)
                      ggCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTAt                                                                                                                              >  1:1540318/1‑150 (MQ=255)
                      ggCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTAt                                                                                                                              <  2:420068/150‑1 (MQ=255)
                        cccAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCCCCTTGTTGCCCGGTGAAGTGCTGGTCCCCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATgg                                                                                                                            >  1:546383/1‑150 (MQ=255)
                            ggCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGa                                                                                                                        >  1:1444663/1‑150 (MQ=255)
                                         cGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGGCTCGGAGttt                                                                                                           <  2:958601/150‑1 (MQ=255)
                                          gTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTg                                                                                                          <  2:818633/150‑1 (MQ=255)
                                          gTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTg                                                                                                          <  2:978748/150‑1 (MQ=255)
                                            tATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGcc                                                                                                        <  2:1106320/150‑1 (MQ=255)
                                            tATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGcc                                                                                                        >  1:698283/1‑150 (MQ=255)
                                             aTCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCg                                                                                                       <  2:449180/150‑1 (MQ=255)
                                               cAGGCATTTCCCCGCACCATATTGCCCCCTGAAGTGCTGGTCACCGCGCAACAACTCGACGAACAATCTTAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGa                                                                                                     >  2:676795/1‑150 (MQ=255)
                                                aGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGaa                                                                                                    >  2:1632486/1‑150 (MQ=255)
                                                       tGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCt                                                                                             >  1:321807/1‑150 (MQ=255)
                                                        gCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCtt                                                                                            <  2:1540318/150‑1 (MQ=255)
                                                               ccTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGTAACGCGGCTTATTCGcc                                                                                     <  1:728080/150‑1 (MQ=255)
                                                                   gttgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGtt                                                                                 >  2:952318/1‑150 (MQ=255)
                                                                    ttgCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTg                                                                                <  1:501107/150‑1 (MQ=255)
                                                                      gCCCGGTGAAGTGCTGGTCACCGCGCAAAATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTAAGGAGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCt                                                                              >  2:265167/1‑150 (MQ=255)
                                                                       cccGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTc                                                                             >  2:336027/1‑150 (MQ=255)
                                                                       cccGGTGAAGTGCTCGTCACCGCTCAACATCTCGAAGATCAATGTGAAACCTTTCTGATGGCGCTAAAACGCGGCAGTGGCCCTGCCCGCCTTTCGGCTATGGCGGAACTCTCGCAGTTAGCCCGAACGCGGCTTATTCGCCCGTTGCTc                                                                             >  2:1239205/1‑150 (MQ=255)
                                                                        ccGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCg                                                                            <  1:31480/150‑1 (MQ=255)
                                                                             gAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGc                                                                       >  2:1588422/1‑150 (MQ=255)
                                                                               aGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCAc                                                                     >  2:48792/1‑150 (MQ=255)
                                                                                   cTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgg                                                                 >  1:105009/1‑150 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                >  1:279853/1‑117 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                <  2:279853/117‑1 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                >  1:273896/1‑117 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                <  2:273896/117‑1 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                               >  1:1124652/1‑150 (MQ=255)
                                                                                     ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                               >  1:737012/1‑150 (MQ=255)
                                                                                       tCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                             >  1:946389/1‑150 (MQ=255)
                                                                                       tCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                             >  1:952266/1‑150 (MQ=255)
                                                                                       tCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                             <  2:698283/150‑1 (MQ=255)
                                                                                       gcACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                             <  1:502558/149‑1 (MQ=255)
                                                                                        cACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                            >  1:960277/1‑150 (MQ=255)
                                                                                        cACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                            >  1:1235893/1‑150 (MQ=255)
                                                                                        cACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                            >  1:163021/1‑150 (MQ=255)
                                                                                          ccGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAAc                                                          >  1:667846/1‑150 (MQ=255)
                                                                                             cgcAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTgg                                                       >  2:167069/1‑150 (MQ=255)
                                                                                             cgcAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTgg                                                       <  2:486923/150‑1 (MQ=255)
                                                                                               cAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTg                                                     <  1:1231205/150‑1 (MQ=255)
                                                                                               cAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTg                                                     >  1:339347/1‑150 (MQ=255)
                                                                                                aaCATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGc                                                    <  2:1279068/150‑1 (MQ=255)
                                                                                                      tcGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTggg                                              >  2:1237857/1‑150 (MQ=255)
                                                                                                      tcGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTggg                                              >  2:223474/1‑150 (MQ=255)
                                                                                                          cgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGcgcg                                          >  1:1538605/1‑150 (MQ=255)
                                                                                                            aTCAATGTGAAACCTTTCTGCTGGCGCTAAACCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGGTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTc                                        >  2:65578/1‑150 (MQ=255)
                                                                                                            aTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTc                                        <  2:1028222/150‑1 (MQ=255)
                                                                                                            aTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTc                                        <  1:570801/150‑1 (MQ=255)
                                                                                                              cAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAg                                      >  1:213549/1‑150 (MQ=255)
                                                                                                              cAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAg                                      >  1:398580/1‑150 (MQ=255)
                                                                                                               aaTGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGt                                     >  1:942174/1‑150 (MQ=255)
                                                                                                               aaTGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGt                                     >  1:699073/1‑150 (MQ=255)
                                                                                                                 tgtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTAt                                   <  1:1632486/150‑1 (MQ=255)
                                                                                                                  gtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATg                                  <  2:1444663/150‑1 (MQ=255)
                                                                                                                  gtgAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGGTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTacg                                  <  1:265167/150‑3 (MQ=255)
                                                                                                                    gAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGAt                                >  2:11355/1‑150 (MQ=255)
                                                                                                                          tttCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGc                          <  1:1290140/150‑1 (MQ=255)
                                                                                                                             ctgctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGATCGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTgg                       <  1:1099248/150‑1 (MQ=255)
                                                                                                                             ctgctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTgg                       <  2:952266/150‑1 (MQ=255)
                                                                                                                             ctgctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTgg                       <  2:946389/150‑1 (MQ=255)
                                                                                                                               gctgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGAt                     >  1:852382/1‑150 (MQ=255)
                                                                                                                                 tgGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTg                   <  2:1538605/150‑1 (MQ=255)
                                                                                                                                  ggCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGa                  <  1:381685/150‑1 (MQ=255)
                                                                                                                                    ctcTAAAACACAGCAGTGGCCCTGCCGGGATTTCGGCTATGGCGGAAGTATAGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAg                >  2:467088/3‑150 (MQ=255)
                                                                                                                                    cgcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAg                >  2:1017783/1‑150 (MQ=255)
                                                                                                                                     gcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGa               >  1:1572341/1‑150 (MQ=255)
                                                                                                                                     gcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGa               <  2:301334/150‑1 (MQ=255)
                                                                                                                                      cTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGAc              >  1:1209554/1‑150 (MQ=255)
                                                                                                                                       tAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACg             <  2:1124652/150‑1 (MQ=255)
                                                                                                                                        aaaaCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGa            <  2:942174/150‑1 (MQ=255)
                                                                                                                                            cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt        >  1:1164826/3‑150 (MQ=255)
                                                                                                                                            cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt        >  2:238580/3‑150 (MQ=255)
                                                                                                                                            cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt        >  1:666425/3‑150 (MQ=255)
                                                                                                                                               ggCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAAt     >  2:511719/1‑150 (MQ=255)
                                                                                                                                                gCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATc    >  2:832735/1‑150 (MQ=255)
                                                                                                                                                  aGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAg  >  2:562793/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGCAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAG  >  NC_000913/212578‑212873

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