Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,815,106 C→A 8.2% intergenic (+16/+21) slmA → / ← pyrE nucleoid occlusion factor SlmA/orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,815,1060CA 8.2% 461.5 / 16.9 141intergenic (+16/+21)slmA/pyrEnucleoid occlusion factor SlmA/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (60/69);  new base A (5/7);  total (65/76)
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 142 total aligned reads displayed.

AAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTA  >  NC_000913/3814957‑3815249
                                                                                                                                                     |                                                                                                                                               
aaaCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGAc                                                                                                                                                 <  1:561129/150‑1 (MQ=255)
 aaCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACt                                                                                                                                                >  1:981215/1‑150 (MQ=255)
 aaCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACt                                                                                                                                                >  1:111295/1‑150 (MQ=255)
  aCCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACtt                                                                                                                                               >  2:521038/1‑150 (MQ=255)
    ggatctgGCTTGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGCCtttt                                                                                                                                             >  2:711687/5‑150 (MQ=255)
    ccTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACtttt                                                                                                                                             >  1:1373210/1‑150 (MQ=255)
    ccTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACtttt                                                                                                                                             >  1:1054953/1‑150 (MQ=255)
     ctgctgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTc                                                                                                                                            >  2:928200/1‑150 (MQ=255)
     ctcctgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTc                                                                                                                                            >  1:217883/4‑150 (MQ=255)
      tgctgGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCa                                                                                                                                           <  1:485327/150‑1 (MQ=255)
        ctgGCAAGCCTGATCCTGGCCTTCTGTGAAGTTATGCTGTCAAGTTTTGTCCGCAGGGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGGCGGATGTCTTTTCATc                                                                                                                                         <  2:1086723/150‑1 (MQ=255)
           gCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCgg                                                                                                                                      <  2:75525/150‑1 (MQ=255)
            cAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGAAAAGTCATCCGGc                                                                                                                                     <  1:1622007/150‑1 (MQ=255)
             aaGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTTGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCg                                                                                                                                    <  2:1219170/150‑1 (MQ=255)
             aaGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGAAAAGTCATCCGGCg                                                                                                                                    <  2:1525534/150‑1 (MQ=255)
             aaGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGAAAAGTCATCCGGCg                                                                                                                                    <  1:1509939/150‑1 (MQ=255)
             aaGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCACAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGAAAAGTCATCCGGCg                                                                                                                                    <  2:513449/150‑1 (MQ=255)
              aGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGa                                                                                                                                   >  2:972272/1‑150 (MQ=255)
               gCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCAGGTTTTGTCGGCAGCGAATTTAAATACAGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAg                                                                                                                                  >  2:430511/1‑150 (MQ=255)
               gCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGAAAAGTCATCCGGCgtc                                                                                                                                  <  1:274859/150‑3 (MQ=255)
                ccAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGAAAAGTCATCCGGCgtca                                                                                                                                 <  1:16733/150‑4 (MQ=255)
                 cAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGtt                                                                                                                                <  1:475845/150‑1 (MQ=255)
                    aTCCTGGCCTTCTGTGAAGGTATGCTGTCACTTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTTCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCt                                                                                                                             <  2:443702/150‑1 (MQ=255)
                     tCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCtt                                                                                                                            >  1:780120/1‑150 (MQ=255)
                          gCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAAc                                                                                                                       <  1:1498703/150‑1 (MQ=255)
                             ttCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGcc                                                                                                                    <  1:1118714/150‑1 (MQ=255)
                               cTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCaa                                                                                                                  >  2:305114/1‑150 (MQ=255)
                                tgtgAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCaaa                                                                                                                 >  2:1374155/1‑150 (MQ=255)
                                 gtgAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAAc                                                                                                                <  2:538443/150‑1 (MQ=255)
                                 gtgAAGGGATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCCCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAAc                                                                                                                <  2:1458784/150‑1 (MQ=255)
                                        tATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTcgc                                                                                                         <  1:65044/150‑1 (MQ=255)
                                        tATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTcgc                                                                                                         >  2:1316750/1‑150 (MQ=255)
                                         aTGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTcgcg                                                                                                        <  1:1146348/150‑1 (MQ=255)
                                           gCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGAt                                                                                                      <  1:451164/150‑1 (MQ=255)
                                             tGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAg                                                                                                    <  2:217883/150‑1 (MQ=255)
                                              gTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAgg                                                                                                   >  2:1010084/1‑150 (MQ=255)
                                                cACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGcc                                                                                                 >  2:800595/1‑150 (MQ=255)
                                                 aCGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCt                                                                                                <  1:1374155/150‑1 (MQ=255)
                                                  cGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCtt                                                                                               <  2:1458420/150‑1 (MQ=255)
                                                   gTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTa                                                                                              <  2:350620/150‑1 (MQ=255)
                                                     tttGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAAc                                                                                            <  1:447803/150‑1 (MQ=255)
                                                     tttGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAAc                                                                                            <  2:553302/150‑1 (MQ=255)
                                                           cGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACcgccgc                                                                                      >  1:1450015/1‑150 (MQ=255)
                                                              aGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAg                                                                                   <  1:249338/150‑1 (MQ=255)
                                                                cGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGAt                                                                                 >  2:1109166/1‑150 (MQ=255)
                                                                 gAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCAACCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATg                                                                                <  1:1386354/150‑1 (MQ=255)
                                                                  aaTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGt                                                                               >  1:848451/1‑150 (MQ=255)
                                                                   aTTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGtt                                                                              >  1:910359/1‑150 (MQ=255)
                                                                      tAAATACCGCCCGACGGATGAGTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTGCTTTAAACGCCAAACTCTGCGCGATAGGCCTTAACCGCCGCCAGATCTTCCg                                                                           >  1:1097472/1‑150 (MQ=255)
                                                                        aaTACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGcc                                                                         <  1:37702/150‑1 (MQ=255)
                                                                        aaTACCGCCCGACGGATGATTTTGACGCCCGCGGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGcc                                                                         <  1:1454202/150‑1 (MQ=255)
                                                                              gaccGACGGATGATTTTGAAGAACGCTGCCCCCTAAATTAGCCACAGTTGCAGTAATATGACGCCAGATGACTTTTCATCCCGACAGATTCTTTAAACGACAAACTCTTCGCGATAGGCATTACCATCCGCAACATGTTACCCCATTTcc                                                                   >  2:158173/3‑150 (MQ=255)
                                                                              gCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTcc                                                                   <  2:111295/150‑1 (MQ=255)
                                                                              gCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTcc                                                                   <  1:1066145/150‑1 (MQ=255)
                                                                              gCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTcc                                                                   >  2:513534/1‑150 (MQ=255)
                                                                               cccGCCGGATGATTTTGACGCCCGCGGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCg                                                                  <  2:542741/150‑1 (MQ=255)
                                                                                cggACGGATGATTTTGTCGCCCGCTGGCCGGTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCgg                                                                 <  2:1543368/148‑1 (MQ=255)
                                                                                    cGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTtc                                                             <  1:5137/150‑1 (MQ=255)
                                                                                    cGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTtc                                                             <  1:642743/150‑1 (MQ=255)
                                                                                     ggATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTtct                                                            >  2:1233902/1‑150 (MQ=255)
                                                                                      gatgatTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTtctc                                                           >  1:450483/1‑150 (MQ=255)
                                                                                          atTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTcc                                                       >  1:522058/1‑150 (MQ=255)
                                                                                           ttttGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCa                                                      <  1:1259895/150‑1 (MQ=255)
                                                                                           ttttGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCa                                                      <  1:28268/150‑1 (MQ=255)
                                                                                           ttttGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCa                                                      <  2:1450015/150‑1 (MQ=255)
                                                                                            tttGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAg                                                     <  1:1276713/150‑1 (MQ=255)
                                                                                             ttGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAgg                                                    <  1:521038/150‑1 (MQ=255)
                                                                                                 cGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAg                                                <  2:828499/150‑1 (MQ=255)
                                                                                                 cGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAg                                                <  2:981215/150‑1 (MQ=255)
                                                                                                  gCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGc                                               <  2:904566/150‑1 (MQ=255)
                                                                                                  gCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGc                                               <  2:1407143/150‑1 (MQ=255)
                                                                                                  gCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGc                                               <  2:953589/150‑1 (MQ=255)
                                                                                                    aactgtgCCGCAATTAGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCTGCCAGATGTTCCGCAATTTACGGCTTCTCTTACAGGTAAGCaa                                             >  1:693901/7‑150 (MQ=255)
                                                                                                     cGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAAt                                            <  2:912356/150‑1 (MQ=255)
                                                                                                      gCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATc                                           >  1:752157/1‑150 (MQ=255)
                                                                                                      ctgtgCCGCAATTAGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATc                                           >  1:831551/5‑150 (MQ=255)
                                                                                                      ctgtgCCGCAATTAGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATc                                           >  2:570572/5‑150 (MQ=255)
                                                                                                       cTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCa                                          >  2:657962/1‑150 (MQ=255)
                                                                                                         ggCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAgg                                        <  1:1564344/150‑1 (MQ=255)
                                                                                                         cgCCGCTAATAGACGCCCAGTTGCAGTAATATGACGCCGGATGAATTTTCATCCTGCGAGTGTCTTTAAACGCCAAACTCTTAGAGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTGCCAGGTAAGCAATCAgg                                        >  2:1133658/2‑150 (MQ=255)
                                                                                                           ccGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTc                                      >  2:1364060/1‑150 (MQ=255)
                                                                                                               tAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTc                                  >  2:126163/1‑150 (MQ=255)
                                                                                                                   tGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAggg                              <  1:1316750/150‑1 (MQ=255)
                                                                                                                    gCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGt                             <  1:972272/150‑1 (MQ=255)
                                                                                                                    gCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGACCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGt                             >  1:1613949/1‑150 (MQ=255)
                                                                                                                      ggCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGtga                           <  2:910359/150‑1 (MQ=255)
                                                                                                                        cacaGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAggg                              >  2:641314/1‑145 (MQ=255)
                                                                                                                         acaGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGtgatga                        <  1:32027/150‑1 (MQ=255)
                                                                                                                           aGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATa                      <  1:559648/150‑1 (MQ=255)
                                                                                                                                cAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGAt                 >  1:1179479/1‑150 (MQ=255)
                                                                                                                                     atatGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACtt            >  2:283327/1‑150 (MQ=255)
                                                                                                                                        tGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGc         <  1:1233902/150‑1 (MQ=255)
                                                                                                                                         gACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCa        >  2:746366/1‑150 (MQ=255)
                                                                                                                                         ctcGCCGGATGAAAAGTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTGCAGGTAAGGAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCa        >  1:1015194/3‑150 (MQ=255)
                                                                                                                                         ctcGCCGGATGAAAAGTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCa        >  2:1058098/3‑150 (MQ=255)
                                                                                                                                         ctcGCCGGATGAAAAGTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCa        >  2:131434/3‑150 (MQ=255)
                                                                                                                                          tcGCCGGATGAAAAGTCATCCGGCGAGTTTCTTTAAACGCCAAACTATTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCGCTTCCAGGTAAGCAATCAGGTCTTTCAGCGTGATGATAGAGATCACTTTGCAg       >  1:1600922/2‑150 (MQ=255)
                                                                                                                                          aCGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAg       <  1:283327/150‑1 (MQ=255)
                                                                                                                                           cGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGt      >  2:1128212/1‑150 (MQ=255)
                                                                                                                                               ggATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTa  <  1:928200/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AAACCCTGCTGGCAAGCCAGATCCTGGCCTTCTGTGAAGGTATGCTGTCACGTTTTGTCCGCAGCGAATTTAAATACCGCCCGACGGATGATTTTGACGCCCGCTGGCCGCTAATTGCGGCACAGTTGCAGTAATATGACGCCGGATGACTTTTCATCCGGCGAGTTTCTTTAAACGCCAAACTCTTCGCGATAGGCCTTAACCGCCGCCAGATGTTCCGCCATTTCCGGCTTCTCTTCCAGGTAAGCAATCAGGTCTTTCAGGGTGATGATAGAGATCACTTTGCAGTTGTA  >  NC_000913/3814957‑3815249

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