Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 212,730 C→A 5.2% P134T (CCT→ACT)  tilS → tRNA(Ile)‑lysidine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913212,7300CA 5.2% 329.3 / 5.2 96P134T (CCT→ACT) tilStRNA(Ile)‑lysidine synthetase
Reads supporting (aligned to +/- strand):  ref base C (50/41);  new base A (3/2);  total (53/43)
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.73e-01

AAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGGACG  >  NC_000913/212584‑212877
                                                                                                                                                  |                                                                                                                                                   
aagGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTg                                                                                                                                                  <  2:386709/150‑1 (MQ=255)
 agGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGc                                                                                                                                                 >  1:19485/1‑150 (MQ=255)
    aCTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGGCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCgg                                                                                                                                              <  1:465891/150‑1 (MQ=255)
      tGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGc                                                                                                                                            <  2:19485/150‑1 (MQ=255)
      tGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGc                                                                                                                                            <  1:353305/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:1016673/150‑1 (MQ=255)
             tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                     <  2:1412790/150‑1 (MQ=255)
             tGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGc                                                                                                                                     <  2:254123/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:273896/117‑1 (MQ=255)
                                                                               ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                    <  2:279853/117‑1 (MQ=255)
                                                                               ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAAcgcg                                                                                                    >  1:279853/1‑117 (MQ=255)
                                                                               ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                                   >  1:737012/1‑150 (MQ=255)
                                                                               ggTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCgggg                                                                   >  1:1124652/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:163021/1‑150 (MQ=255)
                                                                                  cACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                                >  1:1235893/1‑150 (MQ=255)
                                                                                  cACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGa                                                                >  1:960277/1‑150 (MQ=255)
                                                                                    ccGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAAc                                                              >  1:667846/1‑150 (MQ=255)
                                                                                       cgcAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTgg                                                           <  2:486923/150‑1 (MQ=255)
                                                                                       cgcAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTgg                                                           >  2:167069/1‑150 (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:398580/1‑150 (MQ=255)
                                                                                                        cAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAg                                          >  1:213549/1‑150 (MQ=255)
                                                                                                         aaTGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGt                                         >  1:699073/1‑150 (MQ=255)
                                                                                                         aaTGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGt                                         >  1:942174/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            >  1:1164826/3‑150 (MQ=255)
                                                                                                                                      cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt            >  1:666425/3‑150 (MQ=255)
                                                                                                                                      cccgGCAGGGCCACTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGt            >  2:238580/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)
                                                                                                                                                  |                                                                                                                                                   
AAGGACTGGGCATTGAGGCCCAGGCGCGGCAGGCACGTTATCAGGCATTTGCCCGCACCTTGTTGCCCGGTGAAGTGCTGGTCACCGCGCAACATCTCGACGATCAATGTGAAACCTTTCTGCTGGCGCTAAAACGCGGCAGTGGCCCTGCCGGGCTTTCGGCTATGGCGGAAGTCTCGGAGTTTGCCGGAACGCGGCTTATTCGCCCGTTGCTCGCCCGCACGCGGGGGGAACTGGTGCAGTGGGCGCGTCAGTATGATTTACGCTGGATTGAAGACGAAAGTAATCAGGACG  >  NC_000913/212584‑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: