Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,960,578 2 bp→CG coding (1266‑1267/1485 nt) gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,960,5780GC99.8% 1886.0 / NA 490L423V (CTC→GTC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (209/280);  minor base T (0/1);  total (209/281)
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
*NC_0009133,960,5790CG99.8% 1900.8 / NA 489R422R (CGG→CGC) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (207/281);  minor base A (0/1);  total (207/282)
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 525 total aligned reads displayed.

CCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTG  >  NC_000913/3960488‑3960668
                                                                                          ||                                                                                         
ccTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                       <  2:5083816/96‑1 (MQ=255)
ccTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                       <  1:4041746/96‑1 (MQ=255)
ccTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAcgcgcg                                                                                       >  1:6256657/1‑96 (MQ=255)
 cTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcgg                                                                                      >  1:13955695/1‑96 (MQ=255)
  tGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                                                     >  1:1003142/1‑96 (MQ=255)
  tGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                                                     <  1:769693/96‑1 (MQ=255)
  tGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCgcggc                                                                                     <  2:5898903/96‑1 (MQ=255)
     aaTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTGTGCGACGCGCGGCGGt                                                                                  >  1:10320484/1‑96 (MQ=255)
     aaTGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTGTGCGACGCGCGGCGGt                                                                                  >  1:12810696/1‑96 (MQ=255)
       tGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTAc                                                                                <  2:4682423/96‑1 (MQ=255)
        gTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACg                                                                               >  1:1859025/1‑96 (MQ=255)
          cATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                                             <  1:7058821/96‑1 (MQ=255)
          cATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGc                                                                             <  2:4277950/96‑1 (MQ=255)
           aTCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                                            >  1:5079814/1‑96 (MQ=255)
           aTCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCa                                                                            >  2:13840129/1‑96 (MQ=255)
            tctcTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCTACGCGCGGCGGTACGGCAt                                                                           >  1:6538553/1‑96 (MQ=255)
             ctctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTGAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                          <  1:8976568/96‑1 (MQ=255)
             ctctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                          >  1:12095865/1‑96 (MQ=255)
             ctctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtt                                                                          <  1:9854487/96‑1 (MQ=255)
              tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAttt                                                                         <  2:14349126/96‑1 (MQ=255)
              tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAttt                                                                         >  2:5478307/1‑96 (MQ=255)
              tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAttt                                                                         >  2:6141027/1‑96 (MQ=255)
              tctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAAGGATGGCCAGGCGTAGTAAACGGCAGAGCTGTGCTGCGACGCGCGGCGGTACGGCAttt                                                                         >  1:2176709/1‑96 (MQ=255)
               ctGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCAtttt                                                                        >  2:2185012/1‑96 (MQ=255)
                tGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTg                                                                       >  1:8528634/1‑96 (MQ=255)
                tGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTg                                                                       >  1:2890155/1‑96 (MQ=255)
                 ggCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGc                                                                      <  2:2943973/96‑1 (MQ=255)
                    aCGAGATCGTCACGGCTACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                   >  2:9176353/1‑96 (MQ=255)
                    aCGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGa                                                                   <  2:14198101/96‑1 (MQ=255)
                     cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                                                  >  2:12688234/1‑96 (MQ=255)
                     cGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAt                                                                  <  1:13080606/96‑1 (MQ=255)
                      gagaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGGACGGCATTTTGCTGATg                                                                 >  2:9623194/1‑96 (MQ=255)
                       agaTCGTCACAGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                                <  2:13217448/96‑1 (MQ=255)
                       agaTCGGCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGc                                                                <  2:11885078/96‑1 (MQ=255)
                        gaTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCa                                                               >  2:6970740/1‑96 (MQ=255)
                         aTCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAg                                                              <  2:6285250/96‑1 (MQ=255)
                            gTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                           >  1:4731600/1‑96 (MQ=255)
                            gTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                           <  2:4905397/96‑1 (MQ=255)
                            gTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGa                                                           <  2:11737863/96‑1 (MQ=255)
                             tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGAt                                                          <  2:7822702/96‑1 (MQ=255)
                             tCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGAt                                                          <  1:2672229/96‑1 (MQ=255)
                               aCGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATga                                                        <  1:1228202/96‑1 (MQ=255)
                                 ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGGAGCGATgaga                                                      <  1:9245190/96‑1 (MQ=255)
                                 ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                      <  2:3825629/96‑1 (MQ=255)
                                 ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgaga                                                      >  2:9868670/1‑96 (MQ=255)
                                 ggCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGATGTTCTGCGACGCGCGGCGGAACGGCATTTTGCTGATGCAGCGATgaga                                                      >  1:9496444/1‑96 (MQ=255)
                                   cGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATgagaga                                                    <  2:9386163/96‑1 (MQ=255)
                                    gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGAt                                                   <  1:1186450/96‑1 (MQ=255)
                                    gACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGAGGCAGCGATGAGAGAt                                                   >  1:2530196/1‑96 (MQ=255)
                                      gagCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCg                                                 >  1:11686727/3‑96 (MQ=255)
                                       ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                                                <  1:4516240/96‑1 (MQ=255)
                                       ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGa                                                >  2:7265436/1‑96 (MQ=255)
                                       ggCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGCTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGGTGAGAGATCGa                                                >  1:10623800/1‑96 (MQ=255)
                                         cTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGAcc                                              <  1:14346247/96‑1 (MQ=255)
                                           ggCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCgg                                            >  2:13396464/1‑96 (MQ=255)
                                           ggCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGCTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCgg                                            <  1:12704516/96‑1 (MQ=255)
                                             cAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGAt                                          <  1:9969618/96‑1 (MQ=255)
                                                aaaTGATGGCCAGGCGTAGTAAATGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                       <  1:12586704/96‑1 (MQ=255)
                                                aaaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAg                                       >  2:8119567/1‑96 (MQ=255)
                                                 aaTGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGt                                      >  2:12263494/1‑96 (MQ=255)
                                                  atgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTc                                     <  2:8543701/96‑1 (MQ=255)
                                                   tgatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCt                                    >  1:10456360/1‑96 (MQ=255)
                                                    gatgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTg                                   <  2:5607260/96‑1 (MQ=255)
                                                     atgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg                                  <  2:10334111/96‑1 (MQ=255)
                                                     atgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTgg                                  >  2:12381721/1‑96 (MQ=255)
                                                      tgGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGt                                 >  1:12446152/1‑96 (MQ=255)
                                                       ggCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGtt                                >  2:6683072/1‑96 (MQ=255)
                                                        gCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                               <  2:1137565/96‑1 (MQ=255)
                                                        gCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTg                               >  2:1814868/1‑96 (MQ=255)
                                                           aGGCGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACGGGATTAGTCTGGTTGAGc                            >  1:8694744/1‑96 (MQ=255)
                                                              cGTAGTAAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGt                         >  2:11101219/1‑96 (MQ=255)
                                                               gtagtaTATGGCAGAGGTGTTCTGCTACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTa                        <  1:1575806/96‑1 (MQ=255)
                                                               gtagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAgga                        >  1:7498698/1‑94 (MQ=255)
                                                                tagtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAg                       <  2:8772277/96‑1 (MQ=255)
                                                                  gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                     >  2:2670066/1‑96 (MQ=255)
                                                                  gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                     <  1:677624/96‑1 (MQ=255)
                                                                  gtaAACGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCg                     >  1:13370549/1‑96 (MQ=255)
                                                                    aaaGGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                   >  2:12574275/1‑96 (MQ=255)
                                                                    aaaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                   <  2:1470768/96‑1 (MQ=255)
                                                                    aaaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTc                   <  1:5630931/96‑1 (MQ=255)
                                                                     taCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCg                  >  1:11269980/2‑96 (MQ=255)
                                                                     aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGGGCAGTAGCGTCg                  >  1:3616105/1‑96 (MQ=255)
                                                                     aaCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCg                  >  2:10125264/1‑96 (MQ=255)
                                                                      aCGGCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGGTTAGTCTGGTTGAGCAGTAGCGTCGc                 >  1:10932283/1‑96 (MQ=255)
                                                                         gCAGAGTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAg              >  1:5425896/1‑96 (MQ=255)
                                                                          cAGAGTTGTTCTGCGACGCGCGGCGGCACGGCAATTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGc             >  2:14050021/1‑96 (MQ=255)
                                                                             agTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCGAGCAAt          >  2:12937575/1‑96 (MQ=255)
                                                                             agTTGTTCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGGCGCCAGCAAt          >  2:3936170/1‑96 (MQ=255)
                                                                              gttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtt         <  1:8160269/96‑1 (MQ=255)
                                                                              gttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtt         >  1:13542085/1‑96 (MQ=255)
                                                                               ttgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAttt        <  1:5292619/96‑1 (MQ=255)
                                                                                tgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttt       <  2:3398007/96‑1 (MQ=255)
                                                                                tgttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttt       >  2:3061120/1‑96 (MQ=255)
                                                                                  ttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttttt     <  1:3113675/96‑1 (MQ=255)
                                                                                  ttCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAAtttttt     >  2:11558670/1‑96 (MQ=255)
                                                                                   tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTc    >  2:7536782/1‑96 (MQ=255)
                                                                                   tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTc    <  1:4361144/96‑1 (MQ=255)
                                                                                   tCTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTc    >  2:10786512/1‑96 (MQ=255)
                                                                                    cTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCt   <  1:3228740/96‑1 (MQ=255)
                                                                                    cTGCGACGCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCt   <  1:8064595/96‑1 (MQ=255)
                                                                                          ||                                                                                         
CCTGTAATGTCATCTCTGGCACGAGATCGTCACGGCGACGGCTGGCAAAAATGATGGCCAGGCGTAGTAAACGGCAGAGTTGTTCTGCGAGCCGCGGCGGTACGGCATTTTGCTGATGCAGCGATGAGAGATCGACCGGATTAGTCTGGTTGAGCAGTAGCGTCGCCAGCAATTTTTTCTG  >  NC_000913/3960488‑3960668

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑