Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,703,345:1 +A 5.2% coding (42/906 nt) era ← 30S ribosomal subunit maturation GTPase Era

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,703,3451.A 5.2% 797.4 / 36.8 193coding (42/906 nt)era30S ribosomal subunit maturation GTPase Era
Reads supporting (aligned to +/- strand):  ref base . (98/85);  new base A (6/4);  total (104/89)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

Only 100 of 200 total aligned reads displayed.

GTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGG  >  NC_000913/2703199‑2703494
                                                                                                                                                   |                                                                                                                                                     
 tATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCcgacga                                                                                                                                                   >  1:882500/1‑150 (MQ=255)
 tATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCcgacga                                                                                                                                                   >  1:464205/1‑150 (MQ=255)
  aTCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGAt                                                                                                                                                  >  1:383877/1‑150 (MQ=255)
    cGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCATCGTCGGACGTCcgaacgttgg                                                                                                                                               <  1:223814/150‑8 (MQ=255)
    cGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCATCGTCGGACGTCcgaacgttgg                                                                                                                                               <  1:2341169/150‑8 (MQ=255)
    cGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATgg                                                                                                                                                >  1:2100275/1‑150 (MQ=255)
    cGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATgg                                                                                                                                                >  2:1235347/1‑150 (MQ=255)
     gACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGc                                                                                                                                               >  2:1804435/1‑150 (MQ=255)
      aCGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGTAAGTGATGGAGATTTTCTGCCCAAGGAGTTTGTTCAAGAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCa                                                                                                                                              <  2:2183523/150‑1 (MQ=255)
         tAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCATCGTCGGACGTCCGAACGTTGGCaaat                                                                                                                                           <  2:2075516/150‑3 (MQ=255)
         tAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATa                                                                                                                                           >  1:966500/1‑150 (MQ=255)
         tAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATa                                                                                                                                           >  2:867443/1‑150 (MQ=255)
         tAGAGCGCCTGATACCCGCCTTCAGTCTGGATCCCCACAATGCTGTGCCTAGTTGTCTGCGCCTTGCGGGAAGTGCTGGATATTTTCTGCCCCAGCATTTTGTTCAACAATTTTGATTTGCCAACGTTCGGACGTCCGACGATGGCAATa                                                                                                                                           <  2:649884/150‑1 (MQ=255)
          aGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATaa                                                                                                                                          >  2:1084993/1‑150 (MQ=255)
           gATCGCCTGATACGCGCCTTCACTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCAAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATCGCAATaaa                                                                                                                                         >  2:907194/1‑150 (MQ=255)
             tCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATc                                                                                                                                       <  1:30836/150‑1 (MQ=255)
               gCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCg                                                                                                                                     >  2:62425/1‑150 (MQ=255)
               gCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCg                                                                                                                                     >  1:956733/1‑150 (MQ=255)
               gCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCg                                                                                                                                     <  2:2309084/150‑1 (MQ=255)
                ccTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGc                                                                                                                                    >  2:253923/1‑150 (MQ=255)
                ccTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGc                                                                                                                                    >  2:1272621/1‑150 (MQ=255)
                  tGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAg                                                                                                                                  <  1:1503500/150‑1 (MQ=255)
                       cgcgCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACt                                                                                                                             >  2:950753/1‑150 (MQ=255)
                       cgcgCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACt                                                                                                                             >  2:739576/1‑150 (MQ=255)
                        gcgcCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACtt                                                                                                                            >  2:157496/1‑150 (MQ=255)
                           ccTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTa                                                                                                                         <  1:785204/150‑1 (MQ=255)
                             ttCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATc                                                                                                                       >  1:1718117/1‑150 (MQ=255)
                             ttCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATc                                                                                                                       <  2:493440/150‑1 (MQ=255)
                                aGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATTGAGATTTTCTTCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGAt                                                                                                                    <  1:694628/150‑1 (MQ=255)
                                aGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGAt                                                                                                                    >  1:483584/1‑150 (MQ=255)
                                    tGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTc                                                                                                                >  1:2463609/1‑150 (MQ=255)
                                      gATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCAt                                                                                                              >  1:809957/1‑150 (MQ=255)
                                         ccccACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTcc                                                                                                           <  1:363834/150‑1 (MQ=255)
                                          cccACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCa                                                                                                          <  1:2424978/150‑1 (MQ=255)
                                           ccACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAg                                                                                                         <  1:1223398/150‑1 (MQ=255)
                                           ccACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAg                                                                                                         >  2:2258459/1‑150 (MQ=255)
                                           ccACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAg                                                                                                         >  2:7279/1‑150 (MQ=255)
                                           ccACAATGCGGTGACGAGTTGGCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAg                                                                                                         <  1:1214119/150‑1 (MQ=255)
                                               aaTGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTcc                                                                                                     <  1:869767/150‑1 (MQ=255)
                                                aTGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCa                                                                                                    >  2:849466/1‑150 (MQ=255)
                                                     gTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGtttt                                                                                               <  1:1624020/150‑1 (MQ=255)
                                                     gTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGtttt                                                                                               >  1:1793331/1‑150 (MQ=255)
                                                      tGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTTATTCCAGCTCCAGttttt                                                                                              <  2:1976296/150‑1 (MQ=255)
                                                            gTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACg                                                                                        >  2:2290164/1‑150 (MQ=255)
                                                            gTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACg                                                                                        >  1:1472214/1‑150 (MQ=255)
                                                            gTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACg                                                                                        <  1:1381275/150‑1 (MQ=255)
                                                             ttGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAAGAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGc                                                                                       <  1:429209/150‑1 (MQ=255)
                                                              tGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGcc                                                                                      <  2:2100275/150‑1 (MQ=255)
                                                              tGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGcc                                                                                      >  2:329920/1‑150 (MQ=255)
                                                               gTCTGCGCCTTGCGGGAAGTGATTGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCt                                                                                     <  1:1972300/150‑1 (MQ=255)
                                                                 cTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGt                                                                                   >  2:1186448/1‑150 (MQ=255)
                                                                   gcgcCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTc                                                                                 >  2:1941445/1‑150 (MQ=255)
                                                                   gcgcCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTc                                                                                 >  1:601600/1‑150 (MQ=255)
                                                                   gcgcCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTc                                                                                 <  1:2334111/150‑1 (MQ=255)
                                                                     gcCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTcgg                                                                               >  2:775097/1‑150 (MQ=255)
                                                                     gcCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTcgg                                                                               <  2:1637665/150‑1 (MQ=255)
                                                                      ccTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTcggc                                                                              <  1:740292/150‑1 (MQ=255)
                                                                      ccTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTcggc                                                                              >  1:2317983/1‑150 (MQ=255)
                                                                        ttGCGGTAAGTGATGGAGATTTTCTGCCCTAGCTGTTTGTTCAACAATGTGGATTTGCCAACGTGCGGACGTCCGACGTTGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTcggcgg                                                                            <  2:2223149/150‑1 (MQ=255)
                                                                          gCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCa                                                                          >  2:1347816/1‑150 (MQ=255)
                                                                             ggAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGcc                                                                       <  2:1161068/150‑1 (MQ=255)
                                                                             ggAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGcc                                                                       <  2:1885905/150‑1 (MQ=255)
                                                                                 gTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCt                                                                   >  1:1181054/1‑150 (MQ=255)
                                                                                     tGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGc                                                               <  2:1705178/150‑1 (MQ=255)
                                                                                      ggAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATTCTCATTCAATCTCCAGTTTTTTCAACGTCTGTTCGGTGGCATCCTGGTCAggg                                                              <  2:483584/150‑3 (MQ=255)
                                                                                      ggAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGcc                                                              <  1:7279/150‑1 (MQ=255)
                                                                                         gaTTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTa                                                           <  1:849466/150‑1 (MQ=255)
                                                                                                 gCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCt                                                   >  2:1462732/1‑150 (MQ=255)
                                                                                                 gCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCt                                                   >  1:1901501/1‑150 (MQ=255)
                                                                                                  ccccAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCtt                                                  >  1:1436906/1‑150 (MQ=255)
                                                                                                  ccccAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCtt                                                  >  1:1357875/1‑150 (MQ=255)
                                                                                                   cccAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTg                                                 >  1:1817931/1‑150 (MQ=255)
                                                                                                   cccAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTg                                                 >  2:835778/1‑150 (MQ=255)
                                                                                                   cccAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTg                                                 >  2:2143510/1‑150 (MQ=255)
                                                                                                   cccAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTg                                                 >  2:2416476/1‑150 (MQ=255)
                                                                                                    ccAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGa                                                <  1:2062395/150‑1 (MQ=255)
                                                                                                           tttGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCtgtg                                         <  2:392527/150‑1 (MQ=255)
                                                                                                            ttGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACgg                                                     <  1:1050187/137‑1 (MQ=255)
                                                                                                            ttGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACgg                                                     >  2:1050187/1‑137 (MQ=255)
                                                                                                             tGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGcc                                       >  1:1129982/1‑150 (MQ=255)
                                                                                                              gTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCa                                      >  2:2226900/1‑150 (MQ=255)
                                                                                                                tCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAac                                    >  2:1527684/1‑150 (MQ=255)
                                                                                                                tCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAac                                    >  1:1596972/1‑150 (MQ=255)
                                                                                                                tCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAac                                    <  1:261525/150‑1 (MQ=255)
                                                                                                                 caacaaTGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAacc                                   <  1:1239234/150‑1 (MQ=255)
                                                                                                                    caaTGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAaccacc                                >  2:1388696/1‑150 (MQ=255)
                                                                                                                      aTGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCgg                              <  1:446807/150‑1 (MQ=255)
                                                                                                                       ggatttattGCCATCGTCGGACGTCCGAACGTTGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGt                             >  2:2058001/8‑150 (MQ=38)
                                                                                                                        gtgGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGtt                            <  1:2297606/150‑1 (MQ=255)
                                                                                                                         tgGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTc                           <  1:2056873/150‑1 (MQ=255)
                                                                                                                              ttGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCa                      <  1:1404841/150‑1 (MQ=255)
                                                                                                                               tGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAg                     >  2:96645/1‑150 (MQ=255)
                                                                                                                               tGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAg                     <  2:456041/150‑1 (MQ=255)
                                                                                                                               tGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAg                     >  2:2443263/1‑150 (MQ=255)
                                                                                                                                gCCAACGTTCGGAAGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCAAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAgg                    >  2:604424/1‑150 (MQ=255)
                                                                                                                                       ttCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTgc             >  1:1307672/1‑149 (MQ=255)
                                                                                                                                        tCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGAc            >  2:2043223/1‑150 (MQ=255)
                                                                                                                                         cGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGAcc           <  1:1272621/150‑1 (MQ=255)
                                                                                                                                           gACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTg         <  2:1436906/150‑1 (MQ=255)
                                                                                                                                                  gacgaTGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTgg  >  1:2075711/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAATGCGGTGACGAGTTGTCTGCGCCTTGCGGGAAGTGATGGAGATTTTCTGCCCCAGCAGTTTGTTCAACAATGTGGATTTGCCAACGTTCGGACGTCCGACGATGGCAATAAATCCGCAGTAACTTTTATCGATGCTCATTCCAGCTCCAGTTTTTTCAACGCCTGTTCGGCGGCAGCCTGCTCAGCCTTACGACGGCTTGAACCTGTGCCAACCACCGGTTCACTCAGGCCGCTGACCTGGCAGTGG  >  NC_000913/2703199‑2703494

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: