Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,947,713 (T)7→8 intergenic (+30/+44) metV → / ← amiC tRNA‑Met/N‑acetylmuramoyl‑L‑alanine amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,7061.T95.9% 337.5 / 6.1 98intergenic (+23/+51)metV/amiCtRNA‑Met/N‑acetylmuramoyl‑L‑alanine amidase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base T (50/44);  total (54/44)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

CGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAA  >  NC_000913/2947609‑2947803
                                                                                                   |                                                                                                   
cgGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tt                                                                                                  >  1:978096/1‑100 (MQ=255)
 ggggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑ttt                                                                                                 >  2:2044103/1‑100 (MQ=255)
  gggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt                                                                                                <  1:405819/100‑1 (MQ=255)
  gggTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt                                                                                                >  1:1481556/1‑100 (MQ=255)
  gcgTGGAGCAGCCTGGTAGCTCGCCGGGCTCATACCCCGCAGTTCGTCGATTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt                                                                                                <  2:1471033/98‑1 (MQ=18)
   ggTGGAGCAGTCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAATCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑t                                                                                              <  2:2239486/100‑1 (MQ=255)
     tGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑ttt                                                                                            <  1:255645/100‑1 (MQ=255)
      ggAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCTCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑tttt                                                                                           <  2:141659/100‑2 (MQ=255)
      ggAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑tttt                                                                                           >  1:2358313/1‑99 (MQ=255)
      ggAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑tttt                                                                                           <  1:732313/100‑2 (MQ=255)
       gAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑ttttg                                                                                          <  2:285748/100‑3 (MQ=255)
       gAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑ttttg                                                                                          <  1:3380978/100‑3 (MQ=255)
         gcagcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑tttt‑ttttgtt                                                                                        <  2:2122632/100‑5 (MQ=255)
           agcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑TTTTCTTTGATTc                                                                                       >  2:2326492/1‑100 (MQ=255)
            gcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCt                                                                                      >  1:1340502/1‑100 (MQ=255)
            gcCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCt                                                                                      <  1:399290/100‑1 (MQ=255)
            ctcTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCt                                                                                      <  2:805468/98‑1 (MQ=255)
             ccTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTg                                                                                     <  1:1393498/100‑1 (MQ=255)
              cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCtag                                                                                    >  1:805468/1‑98 (MQ=255)
               tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGt                                                                                   <  1:1806185/100‑1 (MQ=255)
                ggTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGtc                                                                                  <  1:436694/100‑1 (MQ=255)
                ggTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGtc                                                                                  >  2:3120074/1‑100 (MQ=255)
                   accTCGTCCGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGtc                                                                                  <  2:3318009/95‑1 (MQ=21)
                   aGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTcc                                                                               >  2:328938/1‑100 (MQ=255)
                     cTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCa                                                                             <  2:2358306/100‑1 (MQ=255)
                       cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATa                                                                           >  1:3108201/1‑100 (MQ=255)
                       cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATa                                                                           <  2:2031330/100‑1 (MQ=255)
                       cgtcgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATa                                                                           >  1:1034753/1‑100 (MQ=255)
                          cgGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATaaaa                                                                        >  2:2568568/1‑100 (MQ=255)
                           gggCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATaaaaa                                                                       >  1:1448271/1‑100 (MQ=255)
                            ggCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATaaaaaa                                                                      >  1:10359/1‑100 (MQ=255)
                             gCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCACCGCCACCGATTATTGAACCCCATAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATaaaaaa                                                                      <  2:280550/99‑1 (MQ=25)
                               tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgcg                                                                   >  2:2567710/1‑100 (MQ=255)
                               tCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgcg                                                                   <  2:1034752/100‑1 (MQ=255)
                                cATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgcgc                                                                  <  2:1719723/100‑1 (MQ=255)
                                cATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAgcgc                                                                  <  1:1951471/100‑1 (MQ=255)
                                cATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGcc                                                                 >  2:769512/1‑100 (MQ=255)
                                 aTAACCCGAAGGTAGTCGGTTCAAATCCGGCGCCCGCAACCAATCATTTAACACCCTAACGGGC‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATaa                                                                          >  2:154671/1‑91 (MQ=21)
                                  tAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCa                                                                <  1:876134/100‑1 (MQ=255)
                                  tAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCa                                                                >  2:2140963/1‑100 (MQ=255)
                                    aCCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCCACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCAtt                                                              >  1:1839285/1‑100 (MQ=255)
                                     cccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTTGTCTCCCATAAAAAAGCGCCAt                                                               >  2:2060172/1‑98 (MQ=255)
                                     cccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTc                                                             >  2:1529021/1‑100 (MQ=255)
                                      ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCa                                                            <  2:3108189/100‑1 (MQ=255)
                                      ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCa                                                            >  2:516336/1‑100 (MQ=255)
                                          aGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAgcg                                                         <  1:1242259/99‑1 (MQ=255)
                                           ggTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGcc                                                       >  2:1115606/1‑100 (MQ=255)
                                            gtcgtcGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCt                                                      >  2:943824/1‑100 (MQ=255)
                                              ttccgatctaaaTCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCttt                                                    <  1:2080337/91‑1 (MQ=255)
                                                 cgGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTa                                                 >  1:1109849/1‑100 (MQ=255)
                                                 cgGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTa                                                 <  1:516337/100‑1 (MQ=255)
                                                 cgGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTa                                                 >  1:2130079/1‑100 (MQ=255)
                                                  ggTTCAAATCCGGCCCCCGCAACCAATTATTGAACACGCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTat                                                >  2:1825655/1‑100 (MQ=255)
                                                  ggTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GTTTTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTa                                                 >  1:3147052/1‑100 (MQ=255)
                                                  ggTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTat                                                >  1:1755828/1‑100 (MQ=255)
                                                   gctCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTatc                                               >  1:2309905/3‑100 (MQ=255)
                                                   gTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTatc                                               <  1:2060176/100‑1 (MQ=255)
                                                   gTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTat                                                >  1:532195/1‑99 (MQ=255)
                                                     tCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTatcat                                             <  1:1115608/100‑1 (MQ=255)
                                                      cAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTatcatc                                            <  1:2567718/100‑1 (MQ=255)
                                                      cAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATcc                                           >  1:276039/1‑100 (MQ=255)
                                                       aaaTCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATcc                                           >  2:2497149/1‑100 (MQ=255)
                                                       aaaTCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATcc                                           <  1:3214097/100‑1 (MQ=255)
                                                        aaTCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATccc                                          >  1:2349675/1‑100 (MQ=255)
                                                         aTCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATcccc                                         <  1:1529024/100‑1 (MQ=255)
                                                         aTCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGG‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATcccc                                         >  1:1653810/1‑100 (MQ=255)
                                                          tCCGGCCCCCGCAACCGATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCt                                        <  2:978095/100‑1 (MQ=255)
                                                          tCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCt                                        <  2:1544680/100‑1 (MQ=255)
                                                          tCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCt                                        <  2:1448268/100‑1 (MQ=255)
                                                             ggCCCCCGCAACCGATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTtct                                     >  1:2287306/1‑100 (MQ=255)
                                                             ggCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGATTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTtct                                     >  2:1893011/1‑100 (MQ=255)
                                                              gCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTtctc                                    <  2:1769494/100‑1 (MQ=255)
                                                              gCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTtctc                                    <  1:1465307/100‑1 (MQ=255)
                                                               cccccGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCg                                   <  2:1002521/100‑1 (MQ=255)
                                                                ccccGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGc                                  <  2:1481553/100‑1 (MQ=255)
                                                                  ccGCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCa                                >  1:3248885/1‑100 (MQ=255)
                                                                    gCAACCAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCACCATTCAGCGCCTTATTATCATCCCCTTCTCGCCAGc                              >  2:2289500/1‑100 (MQ=255)
                                                                        ccAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCg                          >  2:381018/1‑100 (MQ=255)
                                                                        ccAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTAAGCGCCTTTTTATCATCACCTTCTCGCCAGCGTCg                          >  2:381113/1‑100 (MQ=255)
                                                                         cAATTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGc                         >  1:2631353/1‑100 (MQ=255)
                                                                          aaTTATTGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGcc                        <  1:1893015/100‑1 (MQ=255)
                                                                           attattGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGccc                       <  2:1656289/100‑1 (MQ=255)
                                                                           attattGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGccc                       <  2:1714262/100‑1 (MQ=255)
                                                                             tattGAACACCCTAACGGGTGG‑‑TTTT‑TTTGTGTCTGTTCACCCATAAAAAAGCGCCATTCAACGCCTTTTTATCATCCCCTTCTCGCCAGCGTCTCCCCa                     >  2:1841062/1‑100 (MQ=255)
                                                                              attGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCAt                    >  1:2656236/1‑100 (MQ=255)
                                                                                tGAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCg                  <  2:532194/100‑1 (MQ=255)
                                                                                 gAACACCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCgg                 >  2:1137480/1‑100 (MQ=255)
                                                                                    cacCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCaa              <  1:193538/100‑1 (MQ=255)
                                                                                    cacCCTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCaa              >  1:3312169/1‑100 (MQ=255)
                                                                                      cccTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATACAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCaaaa            >  2:1830165/1‑100 (MQ=255)
                                                                                      cccTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCaaaa            <  2:2309899/100‑1 (MQ=255)
                                                                                       ccTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAAt           <  1:2126981/100‑1 (MQ=255)
                                                                                       ccTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAAt           >  2:2809288/1‑100 (MQ=255)
                                                                                       ccTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAAt           <  2:1380440/100‑1 (MQ=255)
                                                                                       ccTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAAt           <  1:1137483/100‑1 (MQ=255)
                                                                                        cTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATa          >  1:556807/1‑100 (MQ=255)
                                                                                        cTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATa          >  2:679676/1‑100 (MQ=255)
                                                                                        cTAACGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATa          >  1:1835974/1‑100 (MQ=255)
                                                                                           aCGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGc       <  1:980753/100‑1 (MQ=255)
                                                                                            cGGGT‑GT‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCt      <  1:227311/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTATTGAACACCCTAACGGGT‑G‑‑TTTT‑TTTGTTTCTGGTCTCCCATAAAAAAGCGCCATTCAGCGCCTTTTTATCATCCCCTTCTCGCCAGCGTCGCCCCATCGGCAAAATACGCTTTAA  >  NC_000913/2947609‑2947803

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