Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 212,728 G→C 5.3% G133A (GGC→GCC)  tilS → tRNA(Ile)‑lysidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913212,7280GC 5.3% 327.2 / 7.0 95G133A (GGC→GCC) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base G (50/40);  new base C (3/2);  total (53/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 = 7.13e-01

CAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGGACG  >  NC_000913/212580‑212877
                                                                                                                                                    |                                                                                                                                                     
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:273896/1‑117 (MQ=255)
                                                                                   ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                    <  2:279853/117‑1 (MQ=255)
                                                                                   ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                    <  2:273896/117‑1 (MQ=255)
                                                                                   ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                    >  1:279853/1‑117 (MQ=255)
                                                                                   ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                                   >  1:1124652/1‑150 (MQ=255)
                                                                                   ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                                   >  1:737012/1‑150 (MQ=255)
                                                                                     tCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                                 >  1:952266/1‑150 (MQ=255)
                                                                                     tCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggggg                                                                 >  1:946389/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:223474/1‑150 (MQ=255)
                                                                                                    tcGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTggg                                                  >  2:1237857/1‑150 (MQ=255)
                                                                                                        cgaTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGcgcg                                              >  1:1538605/1‑150 (MQ=255)
                                                                                                          aTCAATGTGAAACCTTTCTGCTGGCGCTAAACCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGGTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTc                                            >  2:65578/1‑150 (MQ=255)
                                                                                                          aTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTc                                            <  1:570801/150‑1 (MQ=255)
                                                                                                          aTCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTc                                            <  2:1028222/150‑1 (MQ=255)
                                                                                                            cAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAg                                          >  1:398580/1‑150 (MQ=255)
                                                                                                            cAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAg                                          >  1:213549/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                   <  2:301334/150‑1 (MQ=255)
                                                                                                                                   gcTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGa                   >  1:1572341/1‑150 (MQ=255)
                                                                                                                                    cTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGAc                  >  1:1209554/1‑150 (MQ=255)
                                                                                                                                     tAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACg                 <  2:1124652/150‑1 (MQ=255)
                                                                                                                                      aaaaCGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGa                <  2:942174/150‑1 (MQ=255)
                                                                                                                                          cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt            >  2:238580/3‑150 (MQ=255)
                                                                                                                                          cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt            >  1:1164826/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)
                                                                                                                                                   ggCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGgac   >  1:208319/1‑150 (MQ=255)
                                                                                                                                                    gCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGgacg  <  2:1235893/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAAGAAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGGACG  >  NC_000913/212580‑212877

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: