Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 206,330 T→G F402C (TTC→TGC)  dnaE → DNA polymerase III alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913206,3300TG98.7% 294.5 / ‑3.5 75F402C (TTC→TGC) dnaEDNA polymerase III alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (44/30);  total (45/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

AACGGCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTTCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGC  >  NC_000913/206185‑206476
                                                                                                                                                 |                                                                                                                                                  
aaCGGCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGAc                                                                                                                                                <  1:920574/150‑1 (MQ=255)
 aCGGCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGt                                                                                                                                              <  1:148300/151‑1 (MQ=255)
  cGGCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCGCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGtt                                                                                                                                             <  1:177892/151‑1 (MQ=255)
   ggCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTg                                                                                                                                            >  1:1130668/1‑151 (MQ=255)
   ggCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTg                                                                                                                                            >  2:354610/1‑151 (MQ=255)
   ggCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGGTg                                                                                                                                            >  1:920768/1‑151 (MQ=255)
     cGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGAc                                                                                                                                          >  1:1279197/1‑151 (MQ=255)
      gTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACt                                                                                                                                         >  1:1158089/1‑151 (MQ=255)
      gTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACt                                                                                                                                         <  2:971782/151‑1 (MQ=255)
        aCCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACtt                                                                                                                                        <  1:934353/150‑1 (MQ=255)
        aCCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACtt                                                                                                                                        <  2:934278/150‑1 (MQ=255)
          cGGTAGGGCCAGGCCGTGGCTCCGGTGCGAGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTg                                                                                                                                     <  1:424551/151‑1 (MQ=255)
           ggTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGt                                                                                                                                    <  2:417880/151‑1 (MQ=255)
                ggCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATgg                                                                                                                                >  2:1221865/1‑150 (MQ=255)
                ggCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATgg                                                                                                                                >  1:799421/1‑150 (MQ=255)
                 gCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGgag                                                                                                                              <  1:394354/151‑1 (MQ=255)
                        cGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAAcgcg                                                                                                                       <  2:221973/151‑1 (MQ=255)
                         gTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAAcgcg                                                                                                                       >  1:425228/1‑150 (MQ=255)
                           ggCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATc                                                                                                                    <  2:1237795/151‑1 (MQ=255)
                            gCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCa                                                                                                                   <  1:906737/151‑1 (MQ=255)
                               ccGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAg                                                                                                                  <  1:625658/149‑1 (MQ=255)
                                      gggTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAg                                                                                                         <  1:474280/151‑1 (MQ=255)
                                         ttCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCAc                                                                                                      >  1:85397/1‑151 (MQ=255)
                                              tggtggCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGc                                                                                                 >  2:249047/1‑151 (MQ=255)
                                               ggtggCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGc                                                                                                 <  1:188346/150‑1 (MQ=255)
                                                 tggCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCgg                                                                                               >  2:820094/1‑150 (MQ=255)
                                                    ccTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACAt                                                                                           <  2:46128/151‑1 (MQ=255)
                                                       aCGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATg                                                                                          <  1:81481/149‑1 (MQ=255)
                                                        cgcgCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTAc                                                                                       <  2:796773/151‑1 (MQ=255)
                                                          cgcTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACgg                                                                                     >  1:734534/1‑151 (MQ=255)
                                                          cgcTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACgg                                                                                     <  1:294910/151‑1 (MQ=255)
                                                            cTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTc                                                                                   >  2:87230/1‑151 (MQ=255)
                                                             tGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCg                                                                                  >  2:928625/1‑151 (MQ=255)
                                                               aaaaTCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTTCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTg                                                                                >  1:590419/1‑151 (MQ=255)
                                                                   tCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGGTATCGAGCACGTAGCGGACATGTACGGTCGTGATg                                                                             >  1:808774/1‑150 (MQ=255)
                                                                               aTCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGat                                                                <  1:1208848/151‑1 (MQ=255)
                                                                                 ccGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGatc                                                               >  2:373357/1‑150 (MQ=255)
                                                                                  cGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGatcat                                                             <  1:1057425/151‑1 (MQ=255)
                                                                                    cTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGatcatc                                                            >  2:644340/1‑150 (MQ=255)
                                                                                         aTTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGTGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTc                                                      <  2:993153/151‑1 (MQ=255)
                                                                                         aTTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTc                                                      >  2:762969/1‑151 (MQ=255)
                                                                                         aTTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTc                                                      >  2:712949/1‑151 (MQ=255)
                                                                                                   ctgTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATgg                                            >  2:1053157/1‑151 (MQ=255)
                                                                                                     gTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATggc                                           <  2:126845/150‑1 (MQ=255)
                                                                                                       tCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATggcg                                          >  2:814529/1‑149 (MQ=255)
                                                                                                         gAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATggcgg                                         >  2:821529/1‑148 (MQ=255)
                                                                                                         gAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGgcggcg                                       >  1:533844/1‑150 (MQ=255)
                                                                                                         gAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGa                                      <  2:1137071/151‑1 (MQ=255)
                                                                                                            cGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGaa                                     >  2:951325/1‑149 (MQ=255)
                                                                                                               ttCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCgg                                >  1:166304/1‑151 (MQ=255)
                                                                                                                tCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGGAAGCg                                 >  1:826044/1‑149 (MQ=255)
                                                                                                                tCCTTAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGt                               >  2:415792/1‑151 (MQ=255)
                                                                                                                   ttAACCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGAt                            >  2:549614/1‑151 (MQ=255)
                                                                                                                     aaCCCGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGa                             >  1:579075/1‑148 (MQ=255)
                                                                                                                       cccGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGTTCAGGTTATCGAGCCCGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGGACAATGGCGGCGCAAGCGGTGa                             >  2:436952/1‑146 (MQ=255)
                                                                                                                       cccGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTGGGTACAATGGCGGCGAAAGCGGTGATCcgc                        >  2:706685/1‑151 (MQ=255)
                                                                                                                       cccGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATcc                          >  1:658424/1‑149 (MQ=255)
                                                                                                                        ccGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCcgcg                       >  2:992045/1‑151 (MQ=255)
                                                                                                                         cGGAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCcgcg                       >  2:302813/1‑150 (MQ=255)
                                                                                                                          ggAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGAc                     <  2:804350/151‑1 (MQ=255)
                                                                                                                          ggAACGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGAc                     <  2:1248560/151‑1 (MQ=255)
                                                                                                                              cGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGt                   >  2:1131649/1‑149 (MQ=255)
                                                                                                                              cGTGTCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTa                  <  1:2998/150‑1 (MQ=255)
                                                                                                                                tgtCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGc               >  1:783841/1‑151 (MQ=255)
                                                                                                                                 gtCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGcc              >  2:440111/1‑151 (MQ=255)
                                                                                                                                 gtCTCCATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGATATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGGATCGCAGATCATCACCTTCGGTACCATGGCGGCGAAAGCGGTGATCCGCGCCGCAGGc               >  2:1089947/1‑150 (MQ=255)
                                                                                                                                     ccATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCcgc            <  2:90362/149‑1 (MQ=255)
                                                                                                                                      cATGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTg         >  1:76435/1‑151 (MQ=255)
                                                                                                                                       aTGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCcgcg           >  1:51880/1‑148 (MQ=255)
                                                                                                                                       aTGCCTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGc        >  2:504552/1‑151 (MQ=255)
                                                                                                                                          ccTGACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTgg     >  1:397316/1‑151 (MQ=255)
                                                                                                                                             gACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTgggg   <  2:85398/150‑1 (MQ=255)
                                                                                                                                             gACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTgggg   >  1:923076/1‑150 (MQ=255)
                                                                                                                                             gACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTggg    <  2:379811/149‑1 (MQ=255)
                                                                                                                                             gACTGCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGc  >  2:869903/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AACGGCGTACCGGTAGGGCCAGGCCGTGGCTCCGGTGCGGGTTCACTGGTGGCCTACGCGCTGAAAATCACCGACCTCGATCCGCTGGAATTTGACCTGCTGTTCGAACGTTTCCTTAACCCGGAACGTGTCTCCATGCCTGACTTCGACGTTGACTTCTGTATGGAGAAACGCGATCAGGTTATCGAGCACGTAGCGGACATGTACGGTCGTGATGCGGTATCGCAGATCATCACCTTCGGTACAATGGCGGCGAAAGCGGTGATCCGCGACGTAGGCCGCGTGCTGGGGC  >  NC_000913/206185‑206476

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