Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,406,129 T→C 100% intergenic (+144/‑61) purA → / → nsrR adenylosuccinate synthetase/nitric oxide‑sensitive repressor for NO regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,406,1290TC97.5% 1005.3 / 6.9 273intergenic (+144/‑61)purA/nsrRadenylosuccinate synthetase/nitric oxide‑sensitive repressor for NO regulon
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (124/142);  total (129/144)
Fisher's exact test for biased strand distribution p-value = 2.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 273 total aligned reads displayed.

GTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTTC  >  NC_000913/4405982‑4406275
                                                                                                                                                   |                                                                                                                                                  
gTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCtg                                                                                                                                                  >  1:160748/1‑150 (MQ=255)
  aaTTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCtgtg                                                                                                                                                >  2:3678628/1‑150 (MQ=255)
   aTTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACTGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTATATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTgg                                                                                                                                               <  1:2327716/150‑1 (MQ=255)
    ttCTGTTACGCCTGTCATATATTTTGCCTGCCGGGCGAACTGTGTGATACATTGCTGTGTCGGGTAAGCAATTACGCTATCCGACAAAGTTTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGc                                                                                                                                              <  1:1197191/150‑1 (MQ=255)
       tGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTgg                                                                                                                                           <  2:2561818/150‑1 (MQ=255)
       tGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTgg                                                                                                                                           <  2:1286907/150‑1 (MQ=255)
          tACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACTTTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGttt                                                                                                                                        <  1:2624008/150‑1 (MQ=255)
             gCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatc                                                                                                                                     <  1:2324506/150‑1 (MQ=255)
               cTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCGTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatcat                                                                                                                                   <  2:3131634/150‑1 (MQ=255)
                  gCAGATATTTTGCCTGCCGGGCGAACTGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCaa                                                                                                                                <  1:1534705/150‑1 (MQ=255)
                      atatTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAAtata                                                                                                                            <  1:2493892/150‑1 (MQ=255)
                      atatTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCGGGTTTATCATCAAtata                                                                                                                            >  2:1016327/1‑150 (MQ=255)
                         ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                         <  1:3003879/150‑1 (MQ=255)
                         ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                         <  1:2067907/150‑1 (MQ=255)
                         ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                         <  1:3763055/150‑1 (MQ=255)
                         ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                         <  1:2251589/150‑1 (MQ=255)
                          tttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATg                                                                                                                        >  2:3479831/1‑150 (MQ=255)
                           ttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGt                                                                                                                       >  1:846568/1‑150 (MQ=255)
                             gcctgccGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTAt                                                                                                                     >  1:4118825/1‑150 (MQ=255)
                               ctgccGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTAttt                                                                                                                   >  2:2636010/1‑150 (MQ=255)
                                 gccGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTAttttt                                                                                                                 >  2:191694/1‑150 (MQ=255)
                                     ggCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACTCAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTccc                                                                                                             <  2:2816041/150‑1 (MQ=255)
                                     ggCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTccc                                                                                                             <  2:1971626/150‑1 (MQ=255)
                                     ggCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTccc                                                                                                             <  1:219168/150‑1 (MQ=255)
                                      gCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCg                                                                                                            >  1:2107769/1‑150 (MQ=255)
                                      gCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCg                                                                                                            >  2:179564/1‑150 (MQ=255)
                                          aCAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGAttt                                                                                                        >  2:818295/1‑150 (MQ=255)
                                            aGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTcc                                                                                                      >  1:4050981/1‑150 (MQ=255)
                                               gtgATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCtt                                                                                                   <  1:3236462/150‑1 (MQ=255)
                                                tgATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCttt                                                                                                  >  2:3013033/1‑150 (MQ=255)
                                                tgATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCttt                                                                                                  <  1:1369671/150‑1 (MQ=255)
                                                   tACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGa                                                                                               <  2:1294430/150‑1 (MQ=255)
                                                   tACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGa                                                                                               <  1:1766758/150‑1 (MQ=255)
                                                    aCATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAg                                                                                              <  2:387219/150‑1 (MQ=255)
                                                     cATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAgg                                                                                             <  2:2655777/150‑1 (MQ=255)
                                                          cTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGAt                                                                                        >  2:2870616/1‑150 (MQ=255)
                                                          cTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGAt                                                                                        <  1:3392289/150‑1 (MQ=255)
                                                          cTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGAt                                                                                        >  1:591842/1‑150 (MQ=255)
                                                          cTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTGGAGGTTGAt                                                                                        >  1:730551/1‑150 (MQ=255)
                                                           tgtgTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATCTAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGAtg                                                                                       >  1:3114566/1‑150 (MQ=255)
                                                              gtCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGc                                                                                    <  1:986067/150‑1 (MQ=255)
                                                               tCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                   >  1:3530112/1‑150 (MQ=255)
                                                               tCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                   <  1:3445500/150‑1 (MQ=255)
                                                               tCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                   <  1:3179630/150‑1 (MQ=255)
                                                               tCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                   <  2:2161968/150‑1 (MQ=255)
                                                               tCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                   <  2:3827179/150‑1 (MQ=255)
                                                               tCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                   >  2:4057548/1‑150 (MQ=255)
                                                                cGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAg                                                                                  <  1:3920641/150‑1 (MQ=255)
                                                                 gggTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGt                                                                                 >  2:1607253/1‑150 (MQ=255)
                                                                 gggTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTGGAGGTTGATGTGCAGt                                                                                 >  1:1563175/1‑150 (MQ=255)
                                                                   gTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTa                                                                               >  1:1455867/1‑150 (MQ=255)
                                                                    tAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTaa                                                                              <  1:1236975/150‑1 (MQ=255)
                                                                     aaGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAAc                                                                             >  2:614549/1‑150 (MQ=255)
                                                                       gCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGa                                                                           <  2:958701/150‑1 (MQ=255)
                                                                          aTTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGtt                                                                        <  2:1368680/150‑1 (MQ=255)
                                                                          aTTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGtt                                                                        >  2:2059745/1‑150 (MQ=255)
                                                                          aTTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGtt                                                                        <  1:2713255/150‑1 (MQ=255)
                                                                          aTTACGCTATCCGACACAGTGTTAAAGCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCATTATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGtt                                                                        <  2:1444246/150‑1 (MQ=255)
                                                                             aCGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCa                                                                     <  2:1857605/150‑1 (MQ=255)
                                                                             aCGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCa                                                                     <  2:2933221/150‑1 (MQ=255)
                                                                             aCGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCa                                                                     <  2:3896620/150‑1 (MQ=255)
                                                                              cGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCAc                                                                    <  1:3309418/150‑1 (MQ=255)
                                                                               gCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACt                                                                   >  2:157608/1‑150 (MQ=255)
                                                                               gCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACt                                                                   <  2:327157/150‑1 (MQ=255)
                                                                                     cGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattac                                                             <  2:3086198/150‑1 (MQ=255)
                                                                                     cGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTGGATGTGCAGTTAACGAGTTTCACTgattac                                                             >  1:896353/1‑150 (MQ=255)
                                                                                      gACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacg                                                            >  2:900616/1‑150 (MQ=255)
                                                                                      gACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGGTGATGTGCAGTTAACGAGTTTCACTgattacg                                                            >  1:1089324/1‑150 (MQ=255)
                                                                                       acacAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacgg                                                           <  2:2999086/150‑1 (MQ=255)
                                                                                           aGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacggatta                                                       >  1:400072/1‑150 (MQ=255)
                                                                                            gtgtTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacggattac                                                      <  2:3708288/150‑1 (MQ=255)
                                                                                                tAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTgc                                                  <  2:4050981/150‑1 (MQ=255)
                                                                                                tAAATCCTCGCTTTTTTCCGTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTgc                                                  <  2:3752422/150‑1 (MQ=255)
                                                                                                  aaTCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTgcgc                                                <  1:1742304/150‑1 (MQ=255)
                                                                                                       tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATc                                           <  1:101047/150‑1 (MQ=255)
                                                                                                       tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATc                                           <  2:2295569/150‑1 (MQ=255)
                                                                                                        cGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCt                                          >  2:1639733/1‑150 (MQ=255)
                                                                                                        cGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCt                                          <  2:3096613/150‑1 (MQ=255)
                                                                                                        cGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCt                                          >  2:95375/1‑150 (MQ=255)
                                                                                                          cTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTAc                                        <  1:783914/150‑1 (MQ=255)
                                                                                                          cTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTAc                                        <  1:1890528/150‑1 (MQ=255)
                                                                                                               ttccttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGc                                   <  2:506465/150‑1 (MQ=255)
                                                                                                                   ttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTct                               >  1:3645976/1‑149 (MQ=255)
                                                                                                                   ttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCa                               >  1:3014909/1‑150 (MQ=255)
                                                                                                                   ttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCa                               >  2:3150681/1‑150 (MQ=255)
                                                                                                                     ccccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCAtt                             >  2:1270323/1‑150 (MQ=255)
                                                                                                                       ccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGc                           >  2:3648237/1‑150 (MQ=255)
                                                                                                                       ccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGc                           >  2:2569770/1‑150 (MQ=255)
                                                                                                                        cGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGcc                          <  2:2334601/150‑1 (MQ=255)
                                                                                                                        cGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGcc                          <  2:1381591/150‑1 (MQ=255)
                                                                                                                         gAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCa                         <  1:535924/150‑1 (MQ=255)
                                                                                                                           aCTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGa                       <  1:1055917/150‑1 (MQ=255)
                                                                                                                               aaataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAggg                   >  1:1989062/1‑150 (MQ=255)
                                                                                                                               aaataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAggg                   >  1:2445444/1‑150 (MQ=255)
                                                                                                                               aaataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAggg                   >  2:3273017/1‑150 (MQ=255)
                                                                                                                                aataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGc                  <  2:730551/150‑1 (MQ=255)
                                                                                                                                     atTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATg             >  2:3216234/1‑150 (MQ=255)
                                                                                                                                            acacAGCCTGTGTCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTg                                                                                                <  1:2807795/60‑1 (MQ=255)
                                                                                                                                            acacAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTg                                                                                                >  2:2807795/1‑60 (MQ=255)
                                                                                                                                                aGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTTc  >  2:2221327/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTTC  >  NC_000913/4405982‑4406275

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