Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,614,580 G→A 6.1% intergenic (+320/+100) 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,5800GA 6.1% 991.9 / 34.0 304intergenic (+320/+100)yjjV/yjjWputative DNase/putative pyruvate formate lyase activating enzyme
Reads supporting (aligned to +/- strand):  ref base G (197/79);  new base A (9/9);  total (207/97)
Fisher's exact test for biased strand distribution p-value = 6.53e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 309 total aligned reads displayed.

ATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTA  >  NC_000913/4614432‑4614729
                                                                                                                                                    |                                                                                                                                                     
   aTGCCTACAAGTCTGTGCCGACCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGGGTCTTATCATGCCTACAGGTTTGTGCCGACCCGTAGGCCGGATAAGGCGTTCGcgcc                                                                                                                                                   <  2:338544/150‑1 (MQ=255)
   aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc                                                                                                                                                   <  1:3575590/150‑1 (MQ=255)
   aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc                                                                                                                                                   >  2:868500/1‑150 (MQ=255)
   aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc                                                                                                                                                   >  1:3088775/1‑150 (MQ=255)
   aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc                                                                                                                                                   >  1:1919830/1‑150 (MQ=255)
   aTGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGATTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgcc                                                                                                                                                   <  1:478288/150‑1 (MQ=255)
    tGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgccg                                                                                                                                                  >  2:2064487/1‑150 (MQ=255)
    tGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGcgccg                                                                                                                                                  >  2:3164783/1‑150 (MQ=255)
        tACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATc                                                                                                                                              >  2:3679235/1‑150 (MQ=255)
         aCAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATcc                                                                                                                                             <  1:1012243/150‑1 (MQ=255)
          cAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGACTTATCATGCCTACAGCTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCg                                                                                                                                            >  1:824394/1‑150 (MQ=255)
           aaGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCgg                                                                                                                                           >  1:971551/1‑150 (MQ=255)
           aaGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCgg                                                                                                                                           >  1:2712766/1‑150 (MQ=255)
             gTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCa                                                                                                                                         >  2:3535985/1‑150 (MQ=255)
              tCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATCAGGCGTTCGCGCCGCATCCGGccg                                                                                                                                        >  1:1756426/1‑148 (MQ=255)
              tCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAg                                                                                                                                        >  1:3898396/1‑150 (MQ=255)
                tgtgCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGtt                                                                                                                                      >  1:3164215/1‑150 (MQ=255)
                  tgCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTgg                                                                                                                                    >  1:3785295/1‑150 (MQ=255)
                    ccGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTATGTGCCGAACCGTAGGCCGGATAAGGCGTTCCCCCCGCATACGGCAGTTggag                                                                                                                                  >  2:216218/1‑148 (MQ=255)
                      gAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGca                                                                                                                                >  1:3932597/1‑150 (MQ=255)
                        aCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGcaca                                                                                                                              >  1:3749636/1‑150 (MQ=255)
                         ccGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACaa                                                                                                                             >  1:1868623/1‑150 (MQ=255)
                               gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTg                                                                                                                       >  1:2940820/1‑150 (MQ=255)
                                ccGGATGCGGCGTGAACGCCTTATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGa                                                                                                                      >  2:2387867/1‑150 (MQ=14)
                                ccGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGa                                                                                                                      <  1:2033754/150‑1 (MQ=255)
                                 cGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCAGAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGAt                                                                                                                     <  1:3027422/150‑1 (MQ=255)
                                   gATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATCAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                                                                                   >  2:3614974/1‑150 (MQ=255)
                                   gATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGc                                                                                                                   >  1:584965/1‑150 (MQ=255)
                                      aaGGCGTTCACGCCGCATCCGGCAGTCGGCGCACAATGCCTGATGCGACGCTTATCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGAc                                                                                                                <  1:1382873/150‑1 (MQ=17)
                                         gCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGtt                                                                                                                                      >  2:1625578/1‑125 (MQ=255)
                                         gCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGtt                                                                                                                                      <  1:1625578/125‑1 (MQ=255)
                                         gCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCt                                                                                                             >  1:1097811/1‑150 (MQ=255)
                                          cGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCtt                                                                                                            <  2:2777693/150‑1 (MQ=255)
                                           gTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTGGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCCCAATGCCTGAGGCGCCGCGTg                                                                                                           >  2:1497521/1‑150 (MQ=14)
                                            ttCACGCCGCCTCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGCCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGa                                                                                                          <  1:674818/150‑1 (MQ=255)
                                                cgccgcATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGAcgcg                                                                                                      >  2:1838521/1‑150 (MQ=255)
                                                 gccgcATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGt                                                                                                     >  2:1113033/1‑150 (MQ=255)
                                                     cATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCGTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATCAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTa                                                                                                 >  1:1845218/1‑150 (MQ=255)
                                                         cGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                             >  1:3218651/1‑150 (MQ=255)
                                                         cGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAg                                                                                             >  2:2580575/1‑150 (MQ=255)
                                                               tCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTAc                                                                                       >  1:506163/1‑150 (MQ=255)
                                                               tCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCGTATCAGGCCTAc                                                                                       >  2:849984/1‑150 (MQ=255)
                                                                 ggCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACaa                                                                                     >  1:1140050/1‑150 (MQ=255)
                                                                 ggCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACaa                                                                                     >  2:3220731/1‑150 (MQ=255)
                                                                   cgcAGAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGt                                                                                   <  1:260521/150‑1 (MQ=255)
                                                                     cacaaTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCt                                                                                 <  1:1406199/147‑1 (MQ=255)
                                                                        aaTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCtgtg                                                                              >  2:554592/1‑150 (MQ=255)
                                                                         aTGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGc                                                                             >  2:885023/1‑150 (MQ=255)
                                                                         aTGCCTGATGCGACGCTTCTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGc                                                                             <  1:3468347/150‑1 (MQ=255)
                                                                             cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATCAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGac                                                                         >  2:1015927/1‑149 (MQ=255)
                                                                             cTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGaa                                                                         >  2:488603/1‑150 (MQ=255)
                                                                                  gCGACGCTTGTCGCGTCTTATCATGCCTACAGGTGTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTa                                                                    <  2:2010724/150‑1 (MQ=255)
                                                                                   cGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAg                                                                   >  2:895403/1‑150 (MQ=255)
                                                                                   cGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGAGAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAg                                                                   >  2:2417133/1‑150 (MQ=255)
                                                                                    gACGCTTGTCGCGTCTTATCATGCCTACAGGTGTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAgg                                                                  <  2:1251394/150‑1 (MQ=255)
                                                                                       gCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCg                                                               >  1:3529496/1‑150 (MQ=255)
                                                                                       gCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCg                                                               >  1:3840160/1‑150 (MQ=255)
                                                                                            tagcgACTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTAGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATcc                                                          >  2:2654760/3‑150 (MQ=255)
                                                                                            tCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATcc                                                          >  2:3407660/1‑150 (MQ=255)
                                                                                            tCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATcc                                                          >  2:3715048/1‑150 (MQ=255)
                                                                                             cgcgTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCg                                                         <  2:2906663/150‑1 (MQ=255)
                                                                                                gTCTTATCCTGCCTACAGGTTTGTGCCGTCCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCa                                                      <  2:1043880/150‑1 (MQ=255)
                                                                                                gTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCa                                                      <  2:2254709/150‑1 (MQ=255)
                                                                                                      tCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTcac                                                <  1:777611/150‑1 (MQ=255)
                                                                                                        aTGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACaa                                              >  2:767059/1‑150 (MQ=255)
                                                                                                         tGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACaaa                                             >  1:72602/1‑150 (MQ=255)
                                                                                                          gCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAAt                                            >  2:3652154/1‑150 (MQ=255)
                                                                                                            cTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATag                                          >  2:1850090/1‑150 (MQ=255)
                                                                                                             tACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGa                                                                                                          <  1:3145752/85‑1 (MQ=35)
                                                                                                                  gTTTGTGCCGAACCGTAGGCCGTATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGcc                                    <  1:868500/150‑1 (MQ=255)
                                                                                                                      gtgCCGAACCGTAGGCCGGATAAGGCGTTCTCGCCGCATCCGGCAGTGGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTACATGTGTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGaa                                <  2:1672265/150‑1 (MQ=17)
                                                                                                                       tgCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGaaa                               >  1:2288449/1‑150 (MQ=255)
                                                                                                                         ccGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACCAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAata                             >  2:526088/1‑150 (MQ=255)
                                                                                                                         ccGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAata                             >  1:2533641/1‑150 (MQ=255)
                                                                                                                         ccGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAata                             >  1:3866083/1‑150 (MQ=255)
                                                                                                                          cGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCt                                                                                 <  1:94026/97‑1 (MQ=18)
                                                                                                                           gAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAatat                            >  1:1317972/1‑149 (MQ=255)
                                                                                                                           gAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAatat                            <  2:1317972/149‑1 (MQ=255)
                                                                                                                             aCCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCaa                         >  1:3670558/1‑150 (MQ=255)
                                                                                                                              ccGTAGGCCGGATACGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAAc                        <  1:1432307/150‑1 (MQ=255)
                                                                                                                              ccGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATg                                                                                                                           <  1:537361/51‑1 (MQ=35)
                                                                                                                                   ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGc                                                                             >  2:1397912/1‑92 (MQ=18)
                                                                                                                                     ccGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGc                                                                                           <  2:2804228/76‑1 (MQ=255)
                                                                                                                                     ccGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCa                 >  1:2644071/1‑150 (MQ=255)
                                                                                                                                     ccGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCa                 >  1:235398/1‑150 (MQ=255)
                                                                                                                                      cGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCAc                <  1:914306/150‑1 (MQ=255)
                                                                                                                                        gATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCAccc              >  2:3299944/1‑150 (MQ=255)
                                                                                                                                         aTAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAgg                                                                                            <  1:3348796/71‑1 (MQ=17)
                                                                                                                                          tAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCcg            >  2:599781/1‑150 (MQ=255)
                                                                                                                                           aaGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCcgc           >  2:1960127/1‑150 (MQ=255)
                                                                                                                                            aGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGc                                                                 <  2:2984389/95‑1 (MQ=18)
                                                                                                                                              gCGTTCGCGCCGCATCCGGCATTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCa        <  1:2049764/150‑1 (MQ=255)
                                                                                                                                              gCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCa        >  2:453791/1‑150 (MQ=255)
                                                                                                                                                gTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGa                                                                                                                      <  2:1781722/38‑1 (MQ=40)
                                                                                                                                                gTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTACAAGTCTGTGc                                                                             <  2:3856431/79‑1 (MQ=21)
                                                                                                                                                gTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCAcc      >  2:2850047/1‑150 (MQ=255)
                                                                                                                                                gTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCtgt                                                                               >  2:1152318/1‑77 (MQ=21)
                                                                                                                                                 ttCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCt     <  1:3439803/150‑1 (MQ=255)
                                                                                                                                                  tCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCtt    <  2:3938311/150‑1 (MQ=255)
                                                                                                                                                   cgcgCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCttt   <  1:853721/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATCATGCCTACAAGTCTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGGCGCATAATGCCTGATGCGACGCTTGTCGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCGCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACAAGTCTGTGCCGAACCGTAGGCCGTATCCGGCATGTCACAAATAGAGCGCCGGAAATATCAACCGGCTCACCCCGCGCACCTTTA  >  NC_000913/4614432‑4614729

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: