Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,614,444 A→C 6.7% intergenic (+184/+236) yjjV → / ← yjjW putative DNase/putative pyruvate formate lyase activating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,614,4440AC 6.7% 540.2 / 2.7 197intergenic (+184/+236)yjjV/yjjWputative DNase/putative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base A (86/85);  new base C (2/16);  total (89/108)
Fisher's exact test for biased strand distribution p-value = 2.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 197 total aligned reads displayed.

GCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCC  >  NC_000913/4614296‑4614590
                                                                                                                                                    |                                                                                                                                                  
gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAg                                                                                                                                                   <  1:855627/150‑1 (MQ=255)
gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAg                                                                                                                                                   <  1:195152/150‑1 (MQ=255)
gcgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAg                                                                                                                                                   >  1:1941052/1‑150 (MQ=255)
 cgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGt                                                                                                                                                  <  1:653552/150‑1 (MQ=255)
 cgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGt                                                                                                                                                  <  1:1952375/150‑1 (MQ=255)
 cgcACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGt                                                                                                                                                  <  1:1926895/150‑1 (MQ=255)
   cacaATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTGCACGCCGCATCCGGCAGTGGGCGCACAATGCCTGATGCGACGCTTGTCGCGGCGTCTCATGCCTACACGTCt                                                                                                                                                >  2:738228/1‑150 (MQ=14)
     caATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgt                                                                                                                                              <  1:232954/150‑1 (MQ=255)
     caATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgt                                                                                                                                              <  2:785309/150‑1 (MQ=255)
     caATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgt                                                                                                                                              <  2:785173/150‑1 (MQ=255)
      aaTGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCtgtg                                                                                                                                             <  1:117263/150‑1 (MQ=255)
           cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGaa                                                                                                                                        >  1:762868/1‑150 (MQ=255)
           cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGaa                                                                                                                                        >  2:1444336/1‑150 (MQ=255)
           cTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTACAAGTCTGTGCCGaa                                                                                                                                        >  2:854202/1‑150 (MQ=37)
                 cGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAg                                                                                                                                  >  1:1270593/1‑150 (MQ=37)
                 cGACGCTGATCGGGTCTTTTCATGGCTTAAGTTGTTTGCCGTCCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAg                                                                                                                                  <  1:1929647/150‑1 (MQ=11)
                   aCGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGc                                                                                                                                <  2:2034650/150‑1 (MQ=255)
                    cGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGcc                                                                                                                               <  1:1179319/150‑1 (MQ=255)
                       ttAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGa                                                                                                                            <  1:1790945/150‑1 (MQ=255)
                         aaCGCGTCTTCTCATGCCTACAGGTTTGTGCCGACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATa                                                                                                                          <  2:1461079/150‑1 (MQ=255)
                             cgTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGc                                                                                                                      <  2:2762200/150‑1 (MQ=255)
                             cgTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGc                                                                                                                      <  2:1945069/150‑1 (MQ=255)
                              gTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCg                                                                                                                     <  2:305331/150‑1 (MQ=255)
                              gTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCg                                                                                                                     <  1:2126217/150‑1 (MQ=255)
                              gTCTTATCATGCCTACAGGTTTGTGCCGAAACGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGcc                                                                                                                                           >  2:352531/1‑128 (MQ=38)
                               tCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATCCGGCGt                                                                                                                    >  2:2907907/1‑150 (MQ=25)
                                cTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGtt                                                                                                                   <  2:2425003/150‑1 (MQ=255)
                                 ttATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTc                                                                                                                  >  2:314369/1‑150 (MQ=255)
                                  tATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACACGTCTGTGCCGAACCGTAGGCCGTATCCGGCGtggc                                                                                                                 >  2:641980/1‑147 (MQ=11)
                                     cATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                              <  2:2081200/150‑1 (MQ=255)
                                     cATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                              <  1:2261175/150‑1 (MQ=255)
                                     cATGCCTACAGGTGTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcgc                                                                                                              <  1:175882/150‑1 (MQ=255)
                                            aCAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                       <  1:2402158/150‑1 (MQ=255)
                                            aCAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCAGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATcc                                                                                                       <  1:849799/150‑1 (MQ=255)
                                                gtgtgtgCCGACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                   <  2:2828848/147‑1 (MQ=255)
                                                       ccGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACATTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATCAGGCGTTCACGCCGCATAGGGCAGTCgggt                                                                                            >  2:1836886/1‑148 (MQ=255)
                                                        cGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGgcgc                                                                                           >  1:1540857/1‑150 (MQ=255)
                                                             cGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAAt                                                                                      >  2:1905450/1‑150 (MQ=255)
                                                                aGGCCGGATAAGGCGGTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACACGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGcc                                                                                   <  1:2269645/150‑1 (MQ=255)
                                                                 tgCCGGATGCGGCGTGAACGCCTTATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGc                                                                                                                      <  1:2226057/113‑1 (MQ=14)
                                                                  gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCAAAATGCCTg                                                                                 >  1:2448937/1‑150 (MQ=255)
                                                                    cGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGAt                                                                               >  2:2428109/1‑150 (MQ=255)
                                                                     ggATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATg                                                                              <  2:2729352/150‑1 (MQ=255)
                                                                      gATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGc                                                                             >  2:1559951/1‑150 (MQ=255)
                                                                           ggCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGc                                                                        <  2:1158314/150‑1 (MQ=255)
                                                                           ggCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTCACGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGc                                                                        <  2:765558/150‑1 (MQ=255)
                                                                            gCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCt                                                                       >  2:473553/1‑150 (MQ=21)
                                                                            gCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCt                                                                       >  1:2789775/1‑150 (MQ=255)
                                                                                tCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCg                                                                                                                                          >  2:2395606/1‑79 (MQ=255)
                                                                                 cACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTcg                                                                  >  1:151102/1‑150 (MQ=255)
                                                                                  aCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTAGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTcgc                                                                 >  1:2721939/1‑150 (MQ=255)
                                                                                  aCGCCGCATCCGGCAGGTGGCGCACAATGCCTGCTGCGACGCTTGTCGCGTCTTATCATGCCTACACGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTcgc                                                                 <  2:1717038/150‑1 (MQ=255)
                                                                                   cgccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTcgcg                                                                <  1:248541/150‑1 (MQ=255)
                                                                                     ccgcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTc                                                              >  2:1491596/1‑150 (MQ=255)
                                                                                       gcATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCtt                                                            <  2:2727464/150‑1 (MQ=255)
                                                                                        ccgccGGCAGGTGGCGCACACGGCCGGCGGCGACGCGTGCCGCGTCGTAGCAGGCCTACCCGTCTGTGCCGAACCGTAGGGCGGATAAGGCGTTCACGCCGCAGCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTa                                                           <  2:1799701/147‑1 (MQ=14)
                                                                                          tCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATc                                                         <  2:762868/150‑1 (MQ=255)
                                                                                           ccGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                        >  2:2221331/1‑150 (MQ=255)
                                                                                           ccGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCa                                                        >  2:836483/1‑150 (MQ=255)
                                                                                             ggCTGTTGGAGCAGAATGACTGATGCGAGGCTTGTCGCGTCTTATCATGCCTACACGTCTGTGCCGACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATg                                                      <  2:2585719/150‑1 (MQ=14)
                                                                                             ggCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGACCAATACGGCAGGAGGCGCATACTGCCTCACGCGACGCTAGTCGCGGCTTAtaaac                                                      >  2:2550017/1‑146 (MQ=14)
                                                                                                  ttGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGtt                                                                                                                                                 <  2:2106167/54‑2 (MQ=9)
                                                                                                  ttGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGtt                                                                                                                                                 >  1:2106167/1‑53 (MQ=9)
                                                                                                   tGGCGCACAATGCCTGATGCGACGCTTGTAGGGTCTTATCATGCCTACACGTGTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACa                                                <  2:1696657/150‑1 (MQ=255)
                                                                                                    ggCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAg                                               <  1:314369/150‑1 (MQ=255)
                                                                                                     gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGTTGCGACGCTTGTCGCGTCTTAGCATGCCTACAgg                                              >  1:2389238/1‑150 (MQ=255)
                                                                                                     gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAgg                                              >  1:539681/1‑150 (MQ=255)
                                                                                                     gcgcACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAgg                                              >  2:454681/1‑150 (MQ=255)
                                                                                                        cacaCTGCCTGATGCGCCGCTTGGCGCGTCTTATCATGCCTACACGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGttt                                           <  2:2418084/150‑1 (MQ=255)
                                                                                                        cacaATGCCTGCTGCGACGCTTGTCGCGTCTTATCATGCCTTCACGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGttt                                           <  2:1016428/150‑1 (MQ=255)
                                                                                                        cacaATGCCTGATGCGACGCTTTTCGCGTCTTATCATGCCTACATGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTAGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGttt                                           <  1:716146/150‑1 (MQ=255)
                                                                                                            ctGCCTGCTGCGACGCTTGTCGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGc                                       <  1:2256556/149‑1 (MQ=255)
                                                                                                                cTGTTGCGAGGATTCTGGGGTGTTTTCATGCGTACATGTGTGTGCCGTTCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGaa                                   <  2:373161/150‑1 (MQ=14)
                                                                                                                cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGaa                                   >  2:988863/1‑150 (MQ=255)
                                                                                                                 tGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAAc                                  >  2:209719/1‑150 (MQ=255)
                                                                                                                     gCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTa                              >  2:2222595/1‑150 (MQ=255)
                                                                                                                         cGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGGCGGCGCATAATGCCTGATCCGACGCTTGTCGCGTCTTAACATCCCAACAGGTTAGTGCCGAACCGTAGGcc                          >  2:1209704/1‑150 (MQ=17)
                                                                                                                            ttGTCGCGTCTTATCATGCCTACATGTCTGTGCCGAACCGTAGGCCGGATTAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGa                       <  1:2194035/150‑1 (MQ=255)
                                                                                                                            ttGTCGCGTCTTATCAAGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCCGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGa                       >  1:1263427/1‑150 (MQ=255)
                                                                                                                                cgcgTCTTATCATGCCTACACGTCTGTGCCGAACCGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAg                   <  2:1245619/150‑1 (MQ=255)
                                                                                                                                 gcgTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAgg                  >  1:1709323/1‑150 (MQ=255)
                                                                                                                                  cgTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGc                 >  1:1665712/1‑150 (MQ=255)
                                                                                                                                    tattATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGt               <  2:539681/148‑1 (MQ=255)
                                                                                                                                      ttATCATGCCTACAAGTCTGTGCCGACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTc             <  1:1129634/150‑1 (MQ=255)
                                                                                                                                       tATCATGCCCACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTcg            <  1:45893/150‑1 (MQ=255)
                                                                                                                                         tCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTcgcg          >  2:1307864/1‑150 (MQ=255)
                                                                                                                                         tCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTcgcg          >  1:1993383/1‑150 (MQ=255)
                                                                                                                                         tCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTcgcg          >  1:1993232/1‑150 (MQ=255)
                                                                                                                                           ctGCCTACACGTCTGTGCCGATCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCCcgcc        <  1:2689651/149‑1 (MQ=255)
                                                                                                                                           aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc        >  2:1187037/1‑150 (MQ=255)
                                                                                                                                           aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc        >  1:392443/1‑150 (MQ=255)
                                                                                                                                            tGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgccg       >  1:2697008/1‑150 (MQ=255)
                                                                                                                                               cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCAt    >  2:1957622/1‑150 (MQ=255)
                                                                                                                                               cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCAt    >  2:2102419/1‑150 (MQ=255)
                                                                                                                                               cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCAt    >  2:2102458/1‑150 (MQ=255)
                                                                                                                                               cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCAt    >  1:1204955/1‑150 (MQ=255)
                                                                                                                                               cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCAt    >  1:1196920/1‑150 (MQ=255)
                                                                                                                                               cTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGAAGCGACGCTTGTCGCGTCTAATCATCCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCAt    >  1:1205076/1‑150 (MQ=255)
                                                                                                                                                tACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATc   >  2:2240407/1‑150 (MQ=255)
                                                                                                                                                 aCAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATcc  <  2:2789776/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCC  >  NC_000913/4614296‑4614590

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