Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 5,103,291 C→T G251E (GGA→GAA)  ftsZ ← cell division protein FtsZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084635,103,2910CT99.8% 2186.7 / NA 504G251E (GGA→GAA) ftsZcell division protein FtsZ
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (270/232);  minor base A (1/0);  total (272/232)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 511 total aligned reads displayed.

GTCGATCACGGTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCAGGGGACAGGTCCGGACCCGCGGTGATGTTCACCATGATGCCGCGCGCGCCCTGCAGGTTGACGTCTCCCAGCAGCGGGGTGCGGATTGCCGCCTCGGGGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCATGCCCGGACGCT  >  NC_008463/5103143‑5103438
                                                                                                                                                    |                                                                                                                                                   
    aTCACGGTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCcagca                                                                                                                                               <  2:3799670/151‑1 (MQ=255)
     tCACGGTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCcagcag                                                                                                                                              >  1:7783411/1‑151 (MQ=255)
     tCACGGTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCcagcag                                                                                                                                              >  2:3426274/1‑151 (MQ=255)
      cACGGTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCagcagc                                                                                                                                             >  2:4487887/1‑151 (MQ=255)
        cGGTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCgg                                                                                                                                           >  2:2790422/1‑151 (MQ=255)
         ggTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCggg                                                                                                                                          <  2:1619407/151‑1 (MQ=255)
           tGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCgggtt                                                                                                                                        <  2:6366167/151‑3 (MQ=255)
            gCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTg                                                                                                                                       <  2:4198138/151‑1 (MQ=255)
             cccACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGc                                                                                                                                      >  1:4271369/1‑151 (MQ=255)
              ccACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCg                                                                                                                                     <  2:11159546/151‑1 (MQ=255)
              ccACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCg                                                                                                                                     >  1:120962/1‑151 (MQ=255)
               cACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCgg                                                                                                                                    <  1:6777602/151‑1 (MQ=255)
               cACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCgg                                                                                                                                    >  1:7167687/1‑151 (MQ=255)
                aCCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGa                                                                                                                                   >  1:10717935/1‑151 (MQ=255)
                 ccTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGAt                                                                                                                                  <  1:4783808/151‑1 (MQ=255)
                  cTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGAtt                                                                                                                                 <  2:608306/151‑1 (MQ=255)
                        aGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTgccgcc                                                                                                                           <  1:9967851/151‑1 (MQ=255)
                         gTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCt                                                                                                                          >  2:1694747/1‑151 (MQ=255)
                           ggCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCg                                                                                                                        <  1:1462517/151‑1 (MQ=255)
                            gCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACCTCGGAGTACTCGCCCAGGGACCGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCgg                                                                                                                       >  1:4626063/1‑151 (MQ=255)
                                 cTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGcc                                                                                                                  >  2:921829/1‑151 (MQ=255)
                                  tCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCt                                                                                                                 <  2:7810309/151‑1 (MQ=255)
                                      aaGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTcgcg                                                                                                             >  2:6799583/1‑151 (MQ=255)
                                      aaGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTcgcg                                                                                                             <  2:8431229/151‑1 (MQ=255)
                                      aaGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTcgcg                                                                                                             <  1:7965784/151‑1 (MQ=255)
                                       agagatcTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCgg                                                                                                            <  1:11132907/145‑1 (MQ=255)
                                       aGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCgg                                                                                                            <  1:1952790/151‑1 (MQ=255)
                                        gCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCggg                                                                                                           >  1:539471/1‑151 (MQ=255)
                                         cGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGc                                                                                                          >  2:6894062/1‑151 (MQ=255)
                                             cTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACgg                                                                                                      <  2:10413174/151‑1 (MQ=255)
                                              tGTTCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGt                                                                                                     <  2:10527693/151‑1 (MQ=255)
                                                ttCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGGTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTc                                                                                                   >  1:1062527/1‑151 (MQ=255)
                                                ttCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGGGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTc                                                                                                   >  2:8935229/1‑151 (MQ=255)
                                                 tCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCTGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCg                                                                                                  >  2:8736968/1‑151 (MQ=255)
                                                 tCGATGATGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGCTGGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCg                                                                                                  <  2:967989/151‑1 (MQ=255)
                                                   gatgatGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGa                                                                                                >  2:7300704/1‑151 (MQ=255)
                                                   gatgatGTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGGTCGGa                                                                                                >  1:8518554/1‑151 (MQ=255)
                                                     tgatgTTGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGAcc                                                                                              <  1:8718566/151‑1 (MQ=255)
                                                           tGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGgc                                                                                        <  1:2712769/151‑1 (MQ=255)
                                                           tGCCGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGgc                                                                                        <  2:1838780/151‑1 (MQ=255)
                                                              cGACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCa                                                                                     >  1:9508051/1‑151 (MQ=255)
                                                               gACATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAg                                                                                    >  1:6303033/1‑151 (MQ=255)
                                                                 cATCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGcc                                                                                  <  1:5114079/151‑1 (MQ=255)
                                                                  aTCGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCg                                                                                 <  1:9046301/151‑1 (MQ=255)
                                                                    cGGAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGt                                                                               >  1:4495132/1‑151 (MQ=255)
                                                                     ggAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTa                                                                              >  1:3731660/1‑151 (MQ=255)
                                                                      gAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTAc                                                                             <  2:4131549/151‑1 (MQ=255)
                                                                      gAGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTAc                                                                             <  1:5689382/151‑1 (MQ=255)
                                                                       aGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTccc                                                                            >  2:10852306/1‑151 (MQ=255)
                                                                       aGTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTAcc                                                                            >  1:8791150/1‑151 (MQ=255)
                                                                        gTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTAccc                                                                           >  2:3036826/1‑151 (MQ=255)
                                                                        gTACTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTAccc                                                                           >  1:2308653/1‑151 (MQ=255)
                                                                          aCTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCcat                                                                         <  1:3901695/151‑1 (MQ=255)
                                                                           cTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCcatc                                                                        <  2:7162224/151‑1 (MQ=255)
                                                                           cTCGCCCAGGGACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCcatc                                                                        >  1:6482762/1‑151 (MQ=255)
                                                                              gcccagggACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCatcatc                                                                     <  1:1876743/146‑1 (MQ=255)
                                                                                 cagggACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGcc                                                                  <  1:7748119/149‑1 (MQ=255)
                                                                                  agggACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCa                                                                 <  2:1378666/150‑1 (MQ=255)
                                                                                  agggACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCa                                                                 <  2:5323310/150‑1 (MQ=255)
                                                                                    ggACAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATg                                                               <  1:575380/151‑1 (MQ=255)
                                                                                      aCAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGcc                                                             <  1:11023988/151‑1 (MQ=255)
                                                                                       cAGGTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGccc                                                            >  1:9220521/1‑151 (MQ=255)
                                                                                         cgTCCGGACCCGCGGTGATGTTCACCAGGAAGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCCGCCGGTACCCATCATCGCCATGCCCAt                                                          >  2:5749738/2‑151 (MQ=255)
                                                                                          gTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCAtt                                                         >  1:2167859/1‑151 (MQ=255)
                                                                                          gTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCAtt                                                         >  1:4845705/1‑151 (MQ=255)
                                                                                          gTCCGGACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCGTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCAtt                                                         <  1:10074432/151‑1 (MQ=255)
                                                                                              ggACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGc                                                     >  1:595774/1‑151 (MQ=255)
                                                                                              ggACCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGc                                                     <  1:3020862/151‑1 (MQ=255)
                                                                                                aCCCGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTc                                                   <  1:3550200/151‑1 (MQ=255)
                                                                                                  ccGCGGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCAt                                                 >  1:5060227/1‑151 (MQ=255)
                                                                                                    gcgGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGa                                               <  1:5011445/151‑1 (MQ=255)
                                                                                                     cgGTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGAc                                              <  2:7999224/151‑1 (MQ=255)
                                                                                                       gTGATGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACgg                                            >  1:2470578/1‑151 (MQ=255)
                                                                                                          aTGTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCt                                         >  1:7519256/1‑151 (MQ=255)
                                                                                                            gTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTc                                       >  1:10178604/1‑151 (MQ=255)
                                                                                                            gTTCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTc                                       <  2:2126954/151‑1 (MQ=255)
                                                                                                             ttCAACAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCa                                      >  2:4594434/1‑151 (MQ=255)
                                                                                                              tCACCAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCAc                                     <  1:8169903/151‑1 (MQ=255)
                                                                                                                 ccAGGATGCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTc                                  <  2:10717721/151‑1 (MQ=255)
                                                                                                                  caggatgCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCg                                 >  2:8665668/4‑151 (MQ=255)
                                                                                                                   aggatgCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATGTCGCTCATGACGGTCTTCACGTCgg                                >  2:266642/3‑151 (MQ=255)
                                                                                                                      atgCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGa                             >  1:1720634/1‑151 (MQ=255)
                                                                                                                       tgCCGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGaa                            <  1:7253609/151‑1 (MQ=255)
                                                                                                                         ccGCGCGCGCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGt                          <  2:5022677/151‑1 (MQ=255)
                                                                                                                          cgcgcgcgCCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTc                         <  1:4918482/151‑1 (MQ=255)
                                                                                                                          cgcgcgcgCCCTGCAGGTTGACGTCGTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTc                         <  1:5260827/151‑1 (MQ=255)
                                                                                                                            cgcgcgcCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGa                       <  1:6842799/151‑1 (MQ=255)
                                                                                                                               gcgcCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGt                    <  2:2479876/151‑1 (MQ=255)
                                                                                                                               gcgcCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGt                    <  1:9709810/151‑1 (MQ=255)
                                                                                                                                cgcCCTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGtt                   <  2:7733188/151‑1 (MQ=255)
                                                                                                                                 gcCCTGCAGGTTGACGTCGTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTg                  <  1:8242768/151‑1 (MQ=255)
                                                                                                                                 gcCCTGCAGGTGGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTg                  <  2:1197881/151‑1 (MQ=255)
                                                                                                                                  cccTGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGa                 <  1:10385795/151‑1 (MQ=255)
                                                                                                                                     tGCAGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCa              >  2:2899283/1‑151 (MQ=255)
                                                                                                                                        aGGTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCATGc           >  2:325966/1‑151 (MQ=255)
                                                                                                                                          gTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCATGccc         <  1:2083087/151‑1 (MQ=255)
                                                                                                                                          gTTGACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCATGccc         >  1:3370610/1‑151 (MQ=255)
                                                                                                                                             gACGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCGCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGTTCATGCCCGGa      >  2:9063026/1‑151 (MQ=255)
                                                                                                                                              aCGTCTTCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCATGCCCGGAc     >  2:9368305/1‑151 (MQ=255)
                                                                                                                                                 tcttCCAGCAGCGGGTTGCGGATTGCCGCCTCGGTGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCCTTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCATGCCCGGACGCt  >  1:8474154/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCGATCACGGTGCCCACCTTCACAGTGGCGTGCTCGGAAGCGAACTGTTCGATGATGTTGCCGACATCGGAGTACTCGCCAGGGGACAGGTCCGGACCCGCGGTGATGTTCACCATGATGCCGCGCGCGCCCTGCAGGTTGACGTCTCCCAGCAGCGGGGTGCGGATTGCCGCCTCGGGGGCCTCGCGGGCACGGTTCGGACCGCTGGCGCAGCCGGTACCCATCATCGCCATGCCCATTTCGCTCATGACGGTCTTCACGTCGGCGAAGTCGACGTTGATCATGCCCGGACGCT  >  NC_008463/5103143‑5103438

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