Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,978 A→T 6.6% intergenic (‑34/+29) ampC ← / ← frdD beta‑lactamase/fumarate reductase membrane protein FrdD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,9780AT 6.6% 819.2 / 30.5 255intergenic (‑34/+29)ampC/frdDbeta‑lactamase/fumarate reductase membrane protein FrdD
Reads supporting (aligned to +/- strand):  ref base A (130/107);  new base T (8/9);  total (139/116)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

Only 100 of 257 total aligned reads displayed.

TGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGT  >  NC_000913/4378829‑4379127
                                                                                                                                                     |                                                                                                                                                     
tgttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATa                                                                                                                                                       <  2:1002154/150‑1 (MQ=255)
    gCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCg                                                                                                                                                   >  1:1993079/1‑150 (MQ=255)
     ctctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                                                                                                                                                  <  2:559811/150‑1 (MQ=255)
     ctctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                                                                                                                                                  >  1:1940448/1‑150 (MQ=255)
              caacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                                                                                                                                         >  1:1405022/8‑150 (MQ=255)
              cGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCgttttg                                                                                                                                         <  2:1973674/150‑6 (MQ=255)
              cGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCgttttg                                                                                                                                         <  1:1906169/150‑6 (MQ=255)
               ggggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                                                                                                                                        <  2:771299/150‑1 (MQ=255)
               aacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                                                                                                                                        >  1:560141/7‑150 (MQ=255)
                 cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa                                                                                                                                      >  2:483897/5‑150 (MQ=255)
                  gTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAt                                                                                                                                     >  1:1750673/1‑150 (MQ=255)
                  gTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAt                                                                                                                                     >  2:1809167/1‑150 (MQ=255)
                      ttGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGc                                                                                                                                 <  1:1288041/150‑1 (MQ=255)
                       tgtgCGCTGCACATTTTCGTTTATTTGTTCAGGGGCAGCATTTGTGGAGCATTAGGCGGTAATTAATTAGGCGCAGAGCGTCGTTTTGTACATTGGGTCTGGTTTCCATACAAAACGGCCGGCGTTTGGCGGTCAGGATTTAGATTGGCg                                                                                                                                <  2:1816093/150‑1 (MQ=255)
                         tgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGAt                                                                                                                              >  1:57389/1‑150 (MQ=255)
                          gCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATg                                                                                                                             <  2:270442/150‑1 (MQ=255)
                               gCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTa                                                                                                                        >  1:667560/1‑150 (MQ=255)
                                cacaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAg                                                                                                                       <  2:1405022/150‑1 (MQ=255)
                                cacaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAg                                                                                                                       >  2:946117/1‑150 (MQ=255)
                                cacaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAg                                                                                                                       <  2:456259/150‑1 (MQ=255)
                                 acaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGa                                                                                                                      <  1:1415513/150‑1 (MQ=255)
                                  caATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCgg                                                                                                                                             <  2:122104/126‑1 (MQ=255)
                                  caATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAt                                                                                                                     >  1:383736/1‑150 (MQ=255)
                                  caATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGAt                                                                                                                     <  2:1422602/150‑1 (MQ=255)
                                      atCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                                                                                                                                              >  2:166311/1‑121 (MQ=255)
                                        cGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAAc                                                                                                               <  2:1473574/150‑1 (MQ=255)
                                         gTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACg                                                                                                              <  1:631955/150‑1 (MQ=255)
                                           tGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGac                                                                                                            >  1:1526430/1‑150 (MQ=255)
                                              tttGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACAcc                                                                                                         <  1:1573511/150‑1 (MQ=255)
                                                   tgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCa                                                                                                    >  1:102579/1‑150 (MQ=255)
                                                   tgAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCa                                                                                                    >  2:40353/1‑150 (MQ=255)
                                                   tgAGGGGCAGCAAACGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACAcc                                                                                                         <  1:321375/145‑1 (MQ=255)
                                                      ggggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCg                                                                                                 <  1:530611/150‑1 (MQ=255)
                                                       gggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGt                                                                                                >  2:1830954/1‑150 (MQ=255)
                                                        ggCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTg                                                                                               <  2:123899/150‑1 (MQ=255)
                                                          cagcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGAc                                                                                             >  1:283880/1‑150 (MQ=255)
                                                          cagcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGAc                                                                                             <  2:1677506/150‑1 (MQ=255)
                                                            gcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACaa                                                                                           >  1:32390/1‑150 (MQ=255)
                                                            gcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACaa                                                                                           >  1:421259/1‑150 (MQ=255)
                                                               aaTGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTg                                                                                        >  1:933571/1‑150 (MQ=255)
                                                               aaTGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTg                                                                                        >  1:1260311/1‑150 (MQ=255)
                                                                aTGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGt                                                                                       >  1:12348/1‑150 (MQ=255)
                                                                aTGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAATGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGt                                                                                       >  1:818830/1‑150 (MQ=255)
                                                                 tgtgGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTc                                                                                      <  1:1809167/150‑1 (MQ=255)
                                                                    ggAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAgg                                                                                   >  2:1154928/1‑150 (MQ=255)
                                                                    ggAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCACACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAgg                                                                                   >  1:3741/1‑150 (MQ=255)
                                                                      aGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGAt                                                                                 <  1:1027723/150‑1 (MQ=255)
                                                                       gCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATa                                                                                >  2:1311751/1‑150 (MQ=255)
                                                                       gCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATa                                                                                >  1:1319556/1‑150 (MQ=255)
                                                                             ggCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGcc                                                                          >  2:72072/1‑150 (MQ=255)
                                                                             ggCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGcc                                                                          <  2:1064364/150‑1 (MQ=255)
                                                                                ggTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGa                                                                       >  1:1815023/1‑150 (MQ=255)
                                                                                ggTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGa                                                                       >  2:696856/1‑150 (MQ=255)
                                                                                   aattaatAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCg                                                                    <  2:1993079/150‑1 (MQ=255)
                                                                                      taataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCTTAg                                                                 <  2:571827/150‑1 (MQ=255)
                                                                                       aataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGa                                                                >  2:547825/1‑150 (MQ=255)
                                                                                        ataaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaa                                                               >  1:792839/1‑150 (MQ=255)
                                                                                          aaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaaaa                                                             <  2:223582/150‑1 (MQ=255)
                                                                                          aaGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGaaaa                                                             <  1:1817978/150‑1 (MQ=255)
                                                                                            ggCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAcc                                                           <  1:1372231/150‑1 (MQ=255)
                                                                                            ggCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAcc                                                           <  2:1615938/150‑1 (MQ=255)
                                                                                            ggCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAAcc                                                           <  2:1598845/150‑1 (MQ=255)
                                                                                                  gagCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTg                                                     <  1:1544931/150‑1 (MQ=255)
                                                                                                    gCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAGAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGcc                                                   <  1:1457760/150‑1 (MQ=255)
                                                                                                     cgtcgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGccc                                                  >  1:1160241/1‑150 (MQ=255)
                                                                                                      gtcgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCg                                                 >  1:1619207/1‑150 (MQ=255)
                                                                                                       tcgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGc                                                <  2:1533232/150‑1 (MQ=255)
                                                                                                       tcgtTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGc                                                >  2:1117856/1‑150 (MQ=255)
                                                                                                          ttttGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAgg                                             <  1:1154928/150‑1 (MQ=255)
                                                                                                            ttGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTa                                           >  2:1904663/1‑150 (MQ=255)
                                                                                                              gAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACg                                         <  1:1867358/150‑1 (MQ=255)
                                                                                                                   tAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGcc                                                   >  1:993152/1‑135 (MQ=255)
                                                                                                                    aGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAtt                                   <  2:1330069/150‑1 (MQ=255)
                                                                                                                     gggTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAttt                                  >  1:1484789/1‑150 (MQ=255)
                                                                                                                     gggTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGAttt                                  >  1:1430046/1‑150 (MQ=255)
                                                                                                                        tCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCa                               <  2:603774/150‑1 (MQ=255)
                                                                                                                          tGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGa                             >  2:1084071/1‑150 (MQ=255)
                                                                                                                           ggTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGAt                            >  1:177514/1‑150 (MQ=255)
                                                                                                                             tttCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCg                          >  1:1093994/1‑150 (MQ=255)
                                                                                                                             tttCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCg                          >  2:1267545/1‑150 (MQ=255)
                                                                                                                              ttCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGt                         >  2:1102006/1‑150 (MQ=255)
                                                                                                                               tCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCCTg                        <  1:1710195/150‑1 (MQ=255)
                                                                                                                                   tACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATc                    >  2:835056/1‑150 (MQ=255)
                                                                                                                                   tACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATc                    >  1:551713/1‑150 (MQ=255)
                                                                                                                                    taaatccGGCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcg                   >  2:1711990/7‑150 (MQ=255)
                                                                                                                                    aCAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcg                   >  1:68912/1‑150 (MQ=255)
                                                                                                                                     cAAAACGGCCCGCCATGGGCGGGCGGGATTTCCATTGGCGATGGGTTAGATGGTAACGACCCCAAGCGGCGTGACACCTGTCCGGATAGCAGCCAGACCGTTGAAAACCCATTTGCCCGCAGGAGCGTGGATTTTAAGATCTTGctgcgc                  >  2:707137/1‑145 (MQ=255)
                                                                                                                                     cAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgc                  >  2:1878719/1‑150 (MQ=255)
                                                                                                                                     cAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgc                  >  1:852310/1‑150 (MQ=255)
                                                                                                                                       aaaCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGc                       >  1:473711/1‑143 (MQ=255)
                                                                                                                                       aaaCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGc                       <  2:473711/143‑1 (MQ=255)
                                                                                                                                         aCGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtgg              <  1:1068325/150‑1 (MQ=255)
                                                                                                                                          cGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt             <  2:1940448/150‑1 (MQ=255)
                                                                                                                                           ggCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg            >  2:362668/1‑150 (MQ=255)
                                                                                                                                           ggCCCGCCTATGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg            >  1:459182/1‑150 (MQ=255)
                                                                                                                                           ggCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg            >  1:1220602/1‑150 (MQ=255)
                                                                                                                                           ggCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg            <  2:802192/150‑1 (MQ=255)
                                                                                                                                           ggCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggtg            <  2:302393/150‑1 (MQ=255)
                                                                                                                                             cccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCa          <  2:340396/150‑1 (MQ=255)
                                                                                                                                              ccGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCAt         >  1:727360/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACGGTGT  >  NC_000913/4378829‑4379127

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: