Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,161,237 T→G 6.3% intergenic (‑91/+20) ygiQ ← / ← ftsP radical SAM superfamily protein YgiQ/cell division protein required during stress conditions

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,161,2370TG 6.3% 471.4 / 20.0 159intergenic (‑91/+20)ygiQ/ftsPradical SAM superfamily protein YgiQ/cell division protein required during stress conditions
Reads supporting (aligned to +/- strand):  ref base T (74/75);  new base G (3/7);  total (77/82)
Fisher's exact test for biased strand distribution p-value = 3.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 170 total aligned reads displayed.

GTACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGT  >  NC_000913/3161090‑3161386
                                                                                                                                                   |                                                                                                                                                     
gTACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCcctcc                                                                                                                                                     <  1:145002/150‑1 (MQ=255)
     gCGGCCAGTAGAAAAGGTCGCGATCCGGTTGTATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATTATAAAAAAATAATCAAAGGCCGGGGATTAAACGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCggg                                                                                                                                                <  2:2296465/150‑1 (MQ=255)
     gCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCggg                                                                                                                                                <  1:1768272/150‑1 (MQ=255)
       ggCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATCATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCCAAAGGTATTCCATACTCGCCCTCCTCGGGCg                                                                                                                                              >  1:1400225/1‑150 (MQ=255)
        gCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGa                                                                                                                                             >  1:214154/1‑150 (MQ=255)
         ccAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAg                                                                                                                                            <  2:718905/150‑1 (MQ=255)
          cAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGt                                                                                                                                           >  1:370959/1‑150 (MQ=255)
            ggAGAACAGGTCGCGATCCGGTTGGATGAGGGAGATAGAGCTGATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGGCGGAACGAAAGAGAAATCGATAGGTATTCCATTGTCGCCCGTGGAGGGCGAGTAt                                                                                                                                         <  2:494479/150‑1 (MQ=255)
            ggAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTAt                                                                                                                                         <  2:1988350/150‑1 (MQ=255)
            ggAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCGAGGAGGGCGAGTAt                                                                                                                                         <  1:2441953/150‑1 (MQ=255)
            ggAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCACTCCTCGGGCGAGTAt                                                                                                                                         >  2:789586/1‑150 (MQ=255)
             gagaACAGGGCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCGAGGAGGGCGAGTATg                                                                                                                                        <  1:2259012/150‑1 (MQ=255)
               gaACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGaa                                                                                                                                      >  1:998791/1‑150 (MQ=255)
                aaCAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAACGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAg                                                                                                                                     <  1:853674/150‑1 (MQ=255)
                 aCAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGa                                                                                                                                    >  1:1110424/1‑150 (MQ=255)
                       cgcgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACgg                                                                                                                              <  2:1661820/150‑1 (MQ=255)
                         cgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTa                                                                                                                            >  1:2028979/1‑150 (MQ=255)
                         cgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTa                                                                                                                            >  1:1848419/1‑150 (MQ=255)
                         cgATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTa                                                                                                                            >  2:1885770/1‑150 (MQ=255)
                             ccGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCgg                                                                                                                        >  2:1805963/1‑150 (MQ=255)
                               ggTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGAt                                                                                                                      >  1:2180172/1‑150 (MQ=255)
                                gTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGAtt                                                                                                                     >  1:989184/1‑150 (MQ=255)
                                 ttGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTg                                                                                                                    <  2:672001/150‑1 (MQ=255)
                                   ggATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGAc                                                                                                                  >  1:347441/1‑150 (MQ=255)
                                   ggATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGAc                                                                                                                  >  1:2139633/1‑150 (MQ=255)
                                    gATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGAcc                                                                                                                 <  1:1512275/150‑1 (MQ=255)
                                    gATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGAcc                                                                                                                 >  1:1597615/1‑150 (MQ=255)
                                       cAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGCCCAAc                                                                                                              >  1:480361/1‑150 (MQ=255)
                                         gggAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAg                                                                                                            >  1:1484406/1‑150 (MQ=255)
                                         gggAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAg                                                                                                            >  1:1335143/1‑150 (MQ=255)
                                                  agCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAAt                                                                                                   >  1:123449/1‑150 (MQ=255)
                                                     tCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGa                                                                                                >  1:166425/1‑150 (MQ=255)
                                                      cATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAg                                                                                               >  2:772855/1‑150 (MQ=255)
                                                      cATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAg                                                                                               <  1:138533/150‑1 (MQ=255)
                                                      cATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAg                                                                                               >  1:1377614/1‑150 (MQ=255)
                                                       aTTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGc                                                                                              >  2:32267/1‑150 (MQ=255)
                                                             cTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGt                                                                                        >  2:1193865/1‑150 (MQ=255)
                                                                ccAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCg                                                                                     <  2:1684214/150‑1 (MQ=255)
                                                                ccAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCg                                                                                     <  1:2291683/150‑1 (MQ=255)
                                                                ccAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCg                                                                                     >  2:108631/1‑150 (MQ=255)
                                                                 cAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGc                                                                                    >  2:2070681/1‑150 (MQ=255)
                                                                 cAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCACTCCTCGGGCGAGTATGAAGATTACGGTACAGGATTGACAAACAGTTGCCCAATCGAGCCACGGTCAGc                                                                                    >  1:154949/1‑150 (MQ=255)
                                                                  aGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGcc                                                                                   <  2:214154/150‑1 (MQ=255)
                                                                           aaaaaaaTAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGc                                                                          <  2:1267713/150‑1 (MQ=255)
                                                                           aaaaaaaTAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGc                                                                          <  1:1926349/150‑1 (MQ=255)
                                                                            aaaaaaTAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGGCATTTCCAGCg                                                                         >  1:891485/1‑150 (MQ=255)
                                                                            aaaaaaTAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCg                                                                         >  2:2330231/1‑150 (MQ=255)
                                                                            aaaaaaTAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCg                                                                         >  2:1842371/1‑150 (MQ=255)
                                                                                aataatCAAAGGCCGGGGATTATAAGCCCGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTg                                                                     <  1:1005362/150‑1 (MQ=255)
                                                                                  taatCAAAGGCCGGGGATTATAAGGCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGAc                                                                   <  1:630336/150‑1 (MQ=255)
                                                                                  taatCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGAc                                                                   <  2:2360049/150‑1 (MQ=255)
                                                                                    atCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTg                                                                 >  1:503094/1‑150 (MQ=255)
                                                                                     tCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGt                                                                >  2:1436058/1‑150 (MQ=255)
                                                                                      cAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGtt                                                               <  2:847260/150‑1 (MQ=255)
                                                                                      cAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGtt                                                               >  1:51638/1‑150 (MQ=255)
                                                                                        aaGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGa                                                             <  1:2004965/150‑1 (MQ=255)
                                                                                          ggCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAg                                                           <  2:1507889/150‑1 (MQ=255)
                                                                                             cGGGGATTATAAGCCGGAACGAAAGAGAAATCTAAAGGTATTCCATACTCGCCCTCCTCGTGCGAGTATGAAGATTACGGTACCGGATTGACCAAGAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAg                                                        <  2:2288850/150‑1 (MQ=255)
                                                                                                 gATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGGGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCTACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACg                                                    <  2:2458204/150‑1 (MQ=255)
                                                                                                   ttATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACggg                                                  <  2:2086588/150‑1 (MQ=255)
                                                                                                    tataAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGa                                                 <  1:1805963/150‑1 (MQ=255)
                                                                                                        aGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTg                                             <  1:2439312/150‑1 (MQ=255)
                                                                                                         gCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTgc                                            >  2:1466613/1‑150 (MQ=255)
                                                                                                          ccGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTgcg                                           >  2:1343777/1‑150 (MQ=255)
                                                                                                                ctaaaGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAgg                                     <  1:1127639/148‑1 (MQ=255)
                                                                                                                    agagAAATCGAAAGGTATTCCATACTCGTCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAgg                                 <  2:1604116/150‑1 (MQ=255)
                                                                                                                     gagaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGc                                >  2:476325/1‑150 (MQ=255)
                                                                                                                     gagaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGc                                <  2:1848419/150‑1 (MQ=255)
                                                                                                                     gagaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGc                                >  1:1919865/1‑150 (MQ=255)
                                                                                                                       gaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATGGAGTGACTGTCCGCCATTTCCAGCGTGTGACTGTTGAAGTATAACGGGAAGTGCGCCCAGGAAGGCTg                              <  2:2130434/150‑1 (MQ=255)
                                                                                                                       gaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTg                              <  1:2380474/150‑1 (MQ=255)
                                                                                                                       gaAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTg                              >  1:815233/1‑150 (MQ=255)
                                                                                                                        aaaTCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGa                             <  2:1932835/150‑1 (MQ=255)
                                                                                                                         aaTCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGAc                            <  2:347441/150‑1 (MQ=255)
                                                                                                                          aTCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGAcc                           <  2:1140751/150‑1 (MQ=255)
                                                                                                                          aTCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGAcc                           >  1:81054/1‑150 (MQ=255)
                                                                                                                              aaaGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGaaa                       <  2:93415/150‑1 (MQ=255)
                                                                                                                              aaaGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGaaa                       <  1:2464280/150‑1 (MQ=255)
                                                                                                                               aaGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAAt                      <  1:2064303/150‑1 (MQ=255)
                                                                                                                                   atcttcATACTCGCCCGAGGAGGGCGAGTATGGAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGCGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAAc                  >  1:2480187/6‑150 (MQ=255)
                                                                                                                                     ttCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACaa                <  1:1283674/150‑1 (MQ=255)
                                                                                                                                     cttcATACTCGCCCGAGGAGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGc                                                                                    <  2:1952670/79‑1 (MQ=25)
                                                                                                                                     cttcATACTCGCCCGAGGAGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGc                                                                                    >  1:1952670/4‑82 (MQ=25)
                                                                                                                                      ttcATACTCGCCCGAGGAGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAag               >  1:2223174/3‑150 (MQ=255)
                                                                                                                                      tCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAag               >  2:1017161/1‑150 (MQ=255)
                                                                                                                                        cATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAagca             <  2:235131/150‑1 (MQ=255)
                                                                                                                                        cATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAagca             <  1:789586/150‑1 (MQ=255)
                                                                                                                                           aCTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCt          <  2:400906/150‑1 (MQ=255)
                                                                                                                                            cTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATAGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTc         >  2:1412151/1‑150 (MQ=255)
                                                                                                                                              cGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCa       >  1:645358/1‑150 (MQ=255)
                                                                                                                                              cGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCa       <  2:1196205/150‑1 (MQ=255)
                                                                                                                                               gCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGCCCGAAATAAACAAGCAGCTCCAc      >  2:2420905/1‑150 (MQ=255)
                                                                                                                                               gCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCAc      <  1:108631/150‑1 (MQ=255)
                                                                                                                                                cccTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGGCCGAAATAAACAAGCAGCTCCACt     >  1:269922/1‑150 (MQ=255)
                                                                                                                                                cccTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACt     <  2:1575404/150‑1 (MQ=255)
                                                                                                                                                  ctcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAag               <  2:104323/138‑1 (MQ=255)
                                                                                                                                                  ctcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAag               >  1:104323/1‑138 (MQ=255)
                                                                                                                                                  ctcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTg   >  2:2237386/1‑150 (MQ=255)
                                                                                                                                                   tcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGCCCGAAATAAACAAGCAGCTCCACTTGt  >  2:62220/1‑150 (MQ=255)
                                                                                                                                                   tcctcGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGt  >  2:899780/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTACTGCGGCCAGGAGAACAGGTCGCGATCCGGTTGGATCAGGGAGATAGAGCTCATTTTGCTTCCAGAAATGATAAAAAAATAATCAAAGGCCGGGGATTATAAGCCGGAACGAAAGAGAAATCGAAAGGTATTCCATACTCGCCCTCCTCGGGCGAGTATGAAGATTACGGTACCGGATTGACCAACAGTTGCCCAATCGAGCCACGGTCCGCCATTTCCAGCGTCTGACTGTTGAAGTAGAACGGGAAGTGCGCCCAGGAAGGCTGACCGAAATAAACAAGCAGCTCCACTTGT  >  NC_000913/3161090‑3161386

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