Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,799,168 C→A 5.3% intergenic (‑226/+58) pheS ← / ← pheM phenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,1680CA 5.3% 586.8 / 13.0 170intergenic (‑226/+58)pheS/pheMphenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (79/82);  new base A (7/2);  total (86/84)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 174 total aligned reads displayed.

ATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAG  >  NC_000913/1799019‑1799317
                                                                                                                                                     |                                                                                                                                                     
aTTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCAc                                                                                                                                                       <  2:1044375/150‑1 (MQ=255)
 ttttGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTACGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACt                                                                                                                                                      >  2:1641981/1‑150 (MQ=255)
  tttGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTg                                                                                                                                                     <  2:463901/150‑1 (MQ=255)
      gTCCCGGGCTGTTCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAg                                                                                                                                                 <  1:13859/150‑1 (MQ=255)
      gTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAg                                                                                                                                                 >  2:2166139/1‑150 (MQ=255)
      gTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAg                                                                                                                                                 >  1:1473600/1‑150 (MQ=255)
       tCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAgg                                                                                                                                                <  2:2200763/150‑1 (MQ=255)
         ccGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCt                                                                                                                                              >  1:600812/1‑150 (MQ=255)
            ggCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTc                                                                                                                                           >  1:22802/1‑150 (MQ=255)
             gCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGATACAGGCTTTTTATTTTTCAAAACGCGCATAGAAAAAAAGTCTCCCAGTGGAGGCttttt                                                                                                                                          <  2:187061/150‑3 (MQ=37)
               tgtgCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAgg                                                                                                                                        <  1:449500/150‑1 (MQ=255)
                gtgCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCCAGTGGAGGCtttttttg                                                                                                                                       <  1:1499/150‑6 (MQ=38)
                 tgCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGgcg                                                                                                                                      <  2:1659517/150‑1 (MQ=255)
                  gCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGgcgc                                                                                                                                     >  2:366636/1‑150 (MQ=255)
                  gCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGgcgc                                                                                                                                     >  2:714158/1‑150 (MQ=255)
                   ctagAAATTCTCACGTACATAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGACTACACTGGGAGGCTTTAAGGCGCt                                                                                                                                    >  2:908055/4‑150 (MQ=39)
                   cTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCt                                                                                                                                    <  1:2366204/150‑1 (MQ=255)
                    tCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTg                                                                                                                                   <  1:337265/150‑1 (MQ=255)
                     cGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGt                                                                                                                                  <  2:379456/150‑1 (MQ=255)
                            ctcACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGt                                                                                                                           >  1:2015041/1‑150 (MQ=255)
                            ctcACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGt                                                                                                                           <  1:1164418/150‑1 (MQ=255)
                            ctcACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGt                                                                                                                           >  1:229948/1‑150 (MQ=255)
                                  tACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctct                                                                                                                     <  2:1003832/150‑1 (MQ=255)
                                    cTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttc                                                                                                                   >  1:1915248/1‑150 (MQ=255)
                                      tAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttctc                                                                                                                 >  1:553776/1‑150 (MQ=255)
                                      tAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTttctcttctc                                                                                                                 >  2:308038/1‑150 (MQ=255)
                                       aaaTACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCa                                                                                                                >  1:1264166/1‑150 (MQ=255)
                                        aaTACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCAc                                                                                                               <  2:127627/150‑1 (MQ=255)
                                            cGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACgcgc                                                                                                           >  2:469817/1‑150 (MQ=255)
                                            cGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACgcgc                                                                                                           >  2:1527658/1‑150 (MQ=255)
                                            cGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACgcgc                                                                                                           >  1:2092019/1‑150 (MQ=255)
                                            cGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACgcgc                                                                                                           >  1:1762098/1‑150 (MQ=255)
                                              cTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTa                                                                                                         >  2:1429567/1‑150 (MQ=255)
                                               acaggTTTATCCGCGCTGGCAGTTTCCATTCTGGATGCGACAATTACACATGATGATACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAg                                                                                                        >  2:2099560/4‑150 (MQ=255)
                                                ccGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGc                                                                                                       <  1:366636/150‑1 (MQ=255)
                                                       ctcCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTgg                                                                                                >  1:1439313/1‑150 (MQ=255)
                                                         ccGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTATTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGAt                                                                                              >  2:1580747/1‑150 (MQ=255)
                                                              cTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAgg                                                                                         <  2:605492/150‑1 (MQ=255)
                                                              cTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAgg                                                                                         <  2:333518/150‑1 (MQ=255)
                                                                ggCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTg                                                                                       >  1:1858890/1‑150 (MQ=255)
                                                                ggCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTg                                                                                       >  1:96406/1‑150 (MQ=255)
                                                                  ccGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCt                                                                                     >  2:1229925/1‑150 (MQ=255)
                                                                  ccGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCt                                                                                     >  2:977608/1‑150 (MQ=255)
                                                                  ccGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCt                                                                                     >  2:898984/1‑150 (MQ=255)
                                                                   cGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTa                                                                                    <  2:1324651/150‑1 (MQ=255)
                                                                     tgtCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTaaa                                                                                  >  2:1816855/1‑150 (MQ=255)
                                                                        ccAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTa                                                                               >  1:1196521/1‑150 (MQ=255)
                                                                         cAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTaa                                                                              <  2:1890424/150‑1 (MQ=255)
                                                                           gTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTaaaa                                                                            <  1:85448/150‑1 (MQ=255)
                                                                                 cTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAaagaag                                                                     >  1:1788344/1‑150 (MQ=255)
                                                                                  tGCGACAATTACGCCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGAGTCAGGTGCTAAAGTAAAAAaagaag                                                                     <  1:943463/150‑1 (MQ=255)
                                                                                   gCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGc                                                                    <  1:23279/150‑1 (MQ=255)
                                                                                      aCAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGa                                                                 >  2:1381707/1‑150 (MQ=255)
                                                                                       cAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGaa                                                                <  1:1895618/150‑1 (MQ=255)
                                                                                         aTTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGaaaa                                                              >  1:1544744/1‑150 (MQ=255)
                                                                                         aTTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGaaaa                                                              <  2:572161/150‑1 (MQ=255)
                                                                                          ttACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAAt                                                             <  2:1858890/150‑1 (MQ=255)
                                                                                             cacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATagc                                                          <  2:45546/150‑1 (MQ=255)
                                                                                                     gagaCAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGcattca                                                  >  2:1624554/1‑150 (MQ=255)
                                                                                                      agaCAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGcattcat                                                 >  1:1766587/1‑150 (MQ=255)
                                                                                                       gaCAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGcattcatt                                                <  1:1565264/150‑1 (MQ=255)
                                                                                                       gaCAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGcattcatt                                                <  2:1196521/150‑1 (MQ=255)
                                                                                                          aGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgct                                             <  2:1473600/150‑1 (MQ=255)
                                                                                                           ggCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgctt                                            >  2:1961738/1‑150 (MQ=255)
                                                                                                           ggCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCAttgctt                                            >  2:1219593/1‑150 (MQ=255)
                                                                                                              tttttATTTTTCAAAACGCGCATACAAAAACAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCg                                         <  2:424785/150‑1 (MQ=255)
                                                                                                              tttttATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCg                                         >  1:374841/1‑150 (MQ=255)
                                                                                                                   aTTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTAcc                                    >  1:378572/1‑150 (MQ=255)
                                                                                                                       ttCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt                                >  1:2148999/1‑150 (MQ=255)
                                                                                                                       ttCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt                                <  1:665853/150‑1 (MQ=255)
                                                                                                                       ttCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt                                <  2:150850/150‑1 (MQ=255)
                                                                                                                       ttCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt                                <  1:324212/150‑1 (MQ=255)
                                                                                                                       ttCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCtttt                                <  1:972148/150‑1 (MQ=255)
                                                                                                                        tCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTg                               >  2:2085244/1‑150 (MQ=255)
                                                                                                                         cAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTgg                              <  1:1607195/150‑1 (MQ=255)
                                                                                                                            aaCGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTAc                           <  2:1721180/150‑1 (MQ=255)
                                                                                                                              cgcgCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTActc                         <  1:1429567/150‑1 (MQ=255)
                                                                                                                                 gcATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTc                      <  2:553776/150‑1 (MQ=255)
                                                                                                                                   aTACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCaa                    <  2:374841/150‑1 (MQ=255)
                                                                                                                                     aCAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCATCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATACCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCaaaa                  <  2:433139/150‑1 (MQ=255)
                                                                                                                                      cAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAg                 >  1:246193/1‑150 (MQ=255)
                                                                                                                                       cctgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGa                >  1:1065487/5‑150 (MQ=255)
                                                                                                                                       cctgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGa                >  1:255450/5‑150 (MQ=255)
                                                                                                                                        ctgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAc               >  1:1219493/4‑150 (MQ=255)
                                                                                                                                        aaaaaaGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAc               >  1:1727753/1‑150 (MQ=255)
                                                                                                                                        aaaaaaGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAc               >  1:1814973/1‑150 (MQ=255)
                                                                                                                                         tgaaaGCCTCCCAGTGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAcc              >  1:2221715/3‑150 (MQ=255)
                                                                                                                                         aaaaaGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAcc              <  2:328316/150‑1 (MQ=255)
                                                                                                                                         aaaaaGCCTCCACTGGGAGGCTTCCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGAcc              <  1:1046117/150‑1 (MQ=255)
                                                                                                                                            aaGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTTAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCttt           <  1:624430/150‑1 (MQ=255)
                                                                                                                                              gCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAt         >  2:993413/1‑150 (MQ=255)
                                                                                                                                              gCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAt         >  2:2035014/1‑150 (MQ=255)
                                                                                                                                               ccTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAtt        <  2:509666/150‑1 (MQ=255)
                                                                                                                                               ccTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTAtt        <  2:1168259/150‑1 (MQ=255)
                                                                                                                                                cTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTg       >  1:90005/1‑150 (MQ=255)
                                                                                                                                                  ccACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaa     >  1:717160/1‑150 (MQ=255)
                                                                                                                                                  ccACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaa     >  2:2081295/1‑150 (MQ=255)
                                                                                                                                                    aCTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGaaaa   <  2:869989/150‑1 (MQ=255)
                                                                                                                                                     cTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAg  <  2:2323604/150‑1 (MQ=255)
                                                                                                                                                     cTGGGAGGCTTTCAGGCGCTGTTTTCCGTATCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAg  <  1:774127/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTTTGGTCCCGGGCTGTGCTCGAAATTCTCACGTACTTAAATACGCTCCGGTTTCTCCGCGCTGGCCGTGTCCAGTCTGGCTGCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCGGAAAATAGCAGCATTCATTGCTTGCGTTACCTTTTGGTACTCTTCAAAAGACCTTTATTGAAAAG  >  NC_000913/1799019‑1799317

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: