Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,599,156 A→G intergenic (‑70/+37) dnaT ← / ← yjjB primosomal protein I/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,599,1560AG96.9% 559.4 / 5.1 192intergenic (‑70/+37)dnaT/yjjBprimosomal protein I/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base G (91/95);  total (97/95)
Fisher's exact test for biased strand distribution p-value = 2.89e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.31e-01

Only 100 of 209 total aligned reads displayed.

GTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGA  >  NC_000913/4599067‑4599245
                                                                                         |                                                                                         
gTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGagggggt                                                                                        >  1:3799335/1‑92 (MQ=255)
gTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGc                                                                                        >  2:2521397/1‑93 (MQ=255)
gTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGc                                                                                        >  2:1497873/1‑93 (MQ=255)
gTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGc                                                                                        <  1:2794687/93‑1 (MQ=255)
   aaaaCTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGt                                                                                     >  1:694842/1‑93 (MQ=255)
    aaaCTCTGGAAGACATGCTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGtc                                                                                    >  1:2446140/1‑93 (MQ=255)
    aaaCTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGtc                                                                                    >  1:1174742/1‑93 (MQ=255)
    aaaCTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGtc                                                                                    <  1:768930/93‑1 (MQ=255)
    aaaCTCTGAAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGtc                                                                                    <  2:1028204/93‑1 (MQ=255)
     aaCTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGtct                                                                                   >  2:4679476/1‑93 (MQ=255)
        tctGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGtctctc                                                                                <  2:512003/93‑1 (MQ=255)
         ctGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCa                                                                               <  1:357509/93‑1 (MQ=255)
         ctGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCa                                                                               <  2:4431803/93‑1 (MQ=255)
          tGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCAt                                                                              >  2:3426822/1‑93 (MQ=255)
           ggAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATc                                                                             >  2:593019/1‑93 (MQ=255)
           ggAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATc                                                                             <  1:1467738/93‑1 (MQ=255)
           ggAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATc                                                                             <  1:2924609/93‑1 (MQ=255)
           ggAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGTCTCTCATc                                                                             >  2:846785/1‑93 (MQ=255)
                aCATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCgtgt                                                                        <  1:567226/93‑1 (MQ=255)
                aCATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCgtgt                                                                        <  1:330554/93‑1 (MQ=255)
                 cATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGtgtg                                                                       >  2:3629503/1‑93 (MQ=255)
                  aTACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGc                                                                      >  1:2155088/1‑93 (MQ=255)
                   tACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCt                                                                     >  2:161195/1‑93 (MQ=255)
                   tACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCt                                                                     >  2:1110499/1‑93 (MQ=255)
                   tACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGTCTCTCATCCGTGTGCt                                                                     >  2:3132957/1‑93 (MQ=255)
                     cTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTa                                                                   >  1:3767545/1‑93 (MQ=255)
                      tGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTaa                                                                  <  2:1381029/93‑1 (MQ=255)
                      tGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTaa                                                                  >  2:2205019/1‑93 (MQ=255)
                       gTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTGCCTCATCCCAGTGCTTAAg                                                                 >  2:1782289/1‑93 (MQ=255)
                        tttCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGc                                                                >  1:232017/1‑93 (MQ=255)
                         ttatcAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCg                                                               <  1:1857038/90‑1 (MQ=255)
                         ttCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCg                                                               <  1:789644/93‑1 (MQ=255)
                          tctcAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGa                                                              <  1:3392630/93‑1 (MQ=255)
                          tctcAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGa                                                              >  2:3277381/1‑93 (MQ=255)
                           ctcAAATATGGAACGATGAACGAGAAGGGAGCTAACTTTTCTAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAg                                                             >  1:4036230/1‑93 (MQ=255)
                           ctcAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAg                                                             <  2:571539/93‑1 (MQ=255)
                            tcAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAgg                                                            <  1:897123/93‑1 (MQ=255)
                              aaaTATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCa                                                          >  2:219685/1‑93 (MQ=255)
                              aaaTATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCa                                                          >  2:2415920/1‑93 (MQ=255)
                               aaTATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCaa                                                         >  2:4324794/1‑93 (MQ=255)
                               aaTATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCaa                                                         <  2:3907167/93‑1 (MQ=255)
                                atatGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAAt                                                        <  1:2069696/93‑1 (MQ=255)
                                   tGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAAtttt                                                     >  2:478426/1‑93 (MQ=255)
                                      aaCGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTtata                                                  <  1:3519905/93‑1 (MQ=255)
                                       aCGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATAc                                                 >  1:4338918/1‑93 (MQ=255)
                                       aCGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATAc                                                 >  2:3083045/1‑93 (MQ=255)
                                        cGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATAcg                                                <  1:453500/93‑1 (MQ=255)
                                         gATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATAcgc                                               <  1:1036285/93‑1 (MQ=255)
                                         gATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATAcgc                                               >  2:2856687/1‑93 (MQ=255)
                                          aTGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATAcgcg                                              >  1:3700506/1‑93 (MQ=255)
                                           tGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGa                                             >  2:2941781/1‑93 (MQ=255)
                                           tGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGa                                             <  1:255399/93‑1 (MQ=255)
                                            gAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAg                                            <  2:4169407/93‑1 (MQ=255)
                                             aaCGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAgg                                           <  2:676900/93‑1 (MQ=255)
                                             aaCGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAgg                                           <  1:3373911/93‑1 (MQ=255)
                                              aCGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAggg                                          <  2:2105123/93‑1 (MQ=255)
                                              aCGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAggg                                          <  1:3605017/93‑1 (MQ=255)
                                               cGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGgc                                         <  1:4022863/93‑1 (MQ=255)
                                                gagaAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGgcg                                        >  1:4666174/1‑93 (MQ=255)
                                                 agaAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGgcgc                                       >  1:2694004/1‑93 (MQ=255)
                                                   aaGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCtt                                     <  1:1989828/93‑1 (MQ=255)
                                                    aGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTg                                    >  1:977787/1‑93 (MQ=255)
                                                    aGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTg                                    >  1:1742117/1‑93 (MQ=255)
                                                     gggAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGc                                   >  1:606852/1‑93 (MQ=255)
                                                     gggAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGc                                   >  1:2834429/1‑93 (MQ=255)
                                                      ggAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCg                                  >  1:880725/1‑93 (MQ=255)
                                                      ggAGATAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCg                                  >  2:2359810/1‑93 (MQ=255)
                                                       gagaTAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCgg                                 >  2:1182177/1‑93 (MQ=255)
                                                       gagaTAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCgg                                 <  1:1330200/93‑1 (MQ=255)
                                                          aTAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTAc                              <  2:3694273/93‑1 (MQ=255)
                                                           tAACTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACa                             >  1:127677/1‑93 (MQ=255)
                                                             aCTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAAc                           <  2:193920/93‑1 (MQ=255)
                                                              cTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAAcc                          >  2:2247279/1‑93 (MQ=255)
                                                              cTTTTATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAAGGCGCTTGCGGTACAAcc                          >  1:676726/1‑93 (MQ=255)
                                                                tttATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCAt                        >  1:2386474/1‑93 (MQ=255)
                                                                 ttATAGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATa                       >  1:1392959/1‑93 (MQ=255)
                                                                 ttATAGAAACACAAGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATa                       >  2:4502360/1‑93 (MQ=255)
                                                                   ataGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAAt                     <  2:648254/93‑1 (MQ=255)
                                                                   ataGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAAt                     >  2:2169868/1‑93 (MQ=255)
                                                                   ataGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAAt                     >  1:2277221/1‑93 (MQ=255)
                                                                   ataGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAAt                     >  2:1490595/1‑93 (MQ=255)
                                                                   ataGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAAt                     <  2:1331103/93‑1 (MQ=255)
                                                                    taGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATc                    >  2:2645893/1‑93 (MQ=255)
                                                                    taGTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATc                    <  2:3506008/93‑1 (MQ=255)
                                                                      gTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGGTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCa                  <  2:173857/93‑1 (MQ=255)
                                                                      gTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCa                  >  2:4382103/1‑93 (MQ=255)
                                                                      gTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCa                  >  1:3054083/1‑93 (MQ=255)
                                                                      gTAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCa                  <  2:204792/93‑1 (MQ=255)
                                                                       tAACACACGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAg                 <  2:2559374/93‑1 (MQ=255)
                                                                          cacacGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGaa              <  2:743343/93‑1 (MQ=255)
                                                                           acacGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAAt             >  1:2857767/1‑93 (MQ=255)
                                                                           acacGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAAt             >  2:2817195/1‑93 (MQ=255)
                                                                              cGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCTCGAATGGa          <  2:4610429/93‑1 (MQ=255)
                                                                              cGGAGGAGAGGGGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGa          >  1:4481048/1‑93 (MQ=255)
                                                                                  ggagaggggGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGa      >  2:377366/1‑93 (MQ=255)
                                                                                   gagaggggGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGAc     <  1:537041/93‑1 (MQ=255)
                                                                                    agaggggGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGAcc    >  2:3443331/1‑93 (MQ=255)
                                                                                     gaggggGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCg   >  1:1994330/1‑93 (MQ=255)
                                                                                      agggggCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGa  <  2:3737961/93‑1 (MQ=255)
                                                                                      agggggCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGa  <  2:915046/93‑1 (MQ=255)
                                                                                         |                                                                                         
GTCAAAACTCTGGAAGACATACTGTTTCTCAAATATGGAACGATGAACGAGAAGGGAGATAACTTTTATAGTAACACACGGAGGAGAGGAGGCTGTCTCTCATCCGTGTGCTTAAGCGAGGCAATTTTATACGCGAGGGCGCTTGCGGTACAACCATAATCCAGGAATGGAAAGACCGA  >  NC_000913/4599067‑4599245

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