Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,759,803 A→G intergenic (‑221/‑7) yibF ← / → rhsA putative S‑transferase/RhsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,759,8030AG85.0% 683.4 / 100.7 267intergenic (‑221/‑7)yibF/rhsAputative S‑transferase/RhsA
Reads supporting (aligned to +/- strand):  ref base A (23/17);  new base G (124/103);  total (147/120)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.66e-01

Only 100 of 272 total aligned reads displayed.

TTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGG  >  NC_000913/3759711‑3759895
                                                                                            |                                                                                            
ttattTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATATCCCGATAAAACAAGGa                                                                                              <  1:5344264/93‑1 (MQ=255)
   tttCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTga                                                                                           <  2:3290478/93‑1 (MQ=255)
   tttCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATga                                                                                           <  2:6522915/93‑1 (MQ=255)
     tCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTgaga                                                                                         >  2:1401172/1‑93 (MQ=255)
        ctctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAAt                                                                                      >  2:887496/1‑93 (MQ=255)
        ctctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAAt                                                                                      >  1:5107207/1‑93 (MQ=255)
        ctctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAAt                                                                                      >  2:2318879/1‑93 (MQ=255)
        ctctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAAt                                                                                      >  2:3609624/1‑93 (MQ=255)
        ctctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAAt                                                                                      >  1:1725058/1‑93 (MQ=255)
         tctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAAGGTGAGAAAta                                                                                     >  1:3080348/1‑92 (MQ=255)
         tatGATGGATATCTCTCTTCAGGCTTTCTTATCAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATg                                                                                     <  2:2203254/91‑1 (MQ=14)
          ctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGa                                                                                    >  1:4294072/1‑93 (MQ=255)
          ctGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGa                                                                                    <  1:3886475/93‑1 (MQ=255)
            gATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGc                                                                                  >  2:1466529/1‑93 (MQ=255)
             aTGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCg                                                                                 <  2:5917977/93‑1 (MQ=255)
             aTGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCg                                                                                 >  2:1917004/1‑93 (MQ=255)
             aTGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCg                                                                                 <  1:2329371/93‑1 (MQ=255)
               ggATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGa                                                                               <  2:3523446/93‑1 (MQ=255)
               ggATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGa                                                                               <  2:6107973/93‑1 (MQ=255)
                gATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGaa                                                                              <  2:1821830/93‑1 (MQ=255)
                 atatCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGaaa                                                                             <  1:1794792/93‑1 (MQ=255)
                 atatCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGaaa                                                                             >  2:6595745/1‑93 (MQ=255)
                 atatCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGGGAGAAATGAGCGGaaa                                                                             >  2:1544084/1‑93 (MQ=255)
                  tatCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGaaaa                                                                            <  1:5668685/93‑1 (MQ=255)
                   atCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAAc                                                                           >  1:3303994/1‑93 (MQ=255)
                    tctcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAAcc                                                                          >  2:1631169/1‑93 (MQ=255)
                    tctcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAAcc                                                                          >  2:5282628/1‑93 (MQ=255)
                     ctcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCg                                                                         <  1:2165035/93‑1 (MQ=255)
                       cACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATACCCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGc                                                                       <  1:1089285/93‑1 (MQ=255)
                       cACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGc                                                                       <  2:4619866/93‑1 (MQ=255)
                         cTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAg                                                                     >  1:2795376/1‑93 (MQ=255)
                         cTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAg                                                                     >  2:2397827/1‑93 (MQ=255)
                           tCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAgcg                                                                   >  1:2471725/1‑93 (MQ=255)
                            cAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAgcgc                                                                  >  1:2104623/1‑93 (MQ=255)
                             aGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGcgcg                                                                 >  2:4336186/1‑93 (MQ=255)
                             aGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGcgcg                                                                 <  2:5011604/93‑1 (MQ=255)
                              ggCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGt                                                                <  1:1202167/93‑1 (MQ=255)
                               gCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTc                                                               <  1:4112367/93‑1 (MQ=255)
                               gCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTc                                                               <  1:5963708/93‑1 (MQ=255)
                               gCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTc                                                               >  2:5372257/1‑93 (MQ=255)
                                cTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCa                                                              <  1:4035440/93‑1 (MQ=255)
                                cTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCa                                                              >  2:3028155/1‑93 (MQ=255)
                                cTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCa                                                              >  2:111434/1‑93 (MQ=255)
                                 tttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAg                                                             >  2:4269831/1‑93 (MQ=255)
                                  ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAgg                                                            <  2:6240859/93‑1 (MQ=255)
                                  ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAgg                                                            <  1:4626105/93‑1 (MQ=255)
                                   tCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAggg                                                           <  1:1516551/93‑1 (MQ=255)
                                      tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGa                                                        <  1:4997581/93‑1 (MQ=255)
                                        taAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACa                                                      <  2:3976153/93‑1 (MQ=255)
                                           agcTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGa                                                   <  1:2896571/91‑1 (MQ=255)
                                            tCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGAc                                                  >  2:6045833/1‑93 (MQ=255)
                                             cTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACg                                                 <  1:45232/93‑1 (MQ=255)
                                               gTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCa                                               >  2:1050015/1‑93 (MQ=255)
                                                tAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAg                                              >  1:1369414/1‑93 (MQ=255)
                                                gaGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAg                                              <  1:867360/92‑1 (MQ=255)
                                                 aGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGt                                             <  2:2963783/93‑1 (MQ=255)
                                                   ggTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAt                                           <  2:4452577/93‑1 (MQ=255)
                                                     tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATgg                                         <  1:4887701/93‑1 (MQ=255)
                                                     tttCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATgg                                         <  2:705291/93‑1 (MQ=255)
                                                       tCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCg                                       >  1:6616946/1‑93 (MQ=255)
                                                         gCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCggg                                     >  2:335947/1‑92 (MQ=255)
                                                         gCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGt                                     <  1:1825194/93‑1 (MQ=255)
                                                         gCCTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGt                                     >  1:5288453/1‑93 (MQ=255)
                                                         gCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGt                                     >  1:2573154/1‑93 (MQ=255)
                                                          ccTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTa                                    >  1:4339049/1‑93 (MQ=255)
                                                           cTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAg                                   <  2:4810146/93‑1 (MQ=255)
                                                           cTGTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAg                                   <  1:2758721/93‑1 (MQ=255)
                                                           cTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAg                                   >  2:2677344/1‑93 (MQ=255)
                                                            tGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGc                                  <  1:5871176/93‑1 (MQ=255)
                                                             gTCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCa                                 <  2:4552918/93‑1 (MQ=255)
                                                              tCAGCAGACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAt                                <  1:1277491/93‑1 (MQ=255)
                                                               cagcagACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAtt                               >  2:4198027/1‑93 (MQ=255)
                                                                agcagACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTg                              >  1:708958/1‑93 (MQ=255)
                                                                agcagACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTg                              >  2:4982071/1‑93 (MQ=255)
                                                                   agACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTc                           >  2:5766408/1‑93 (MQ=255)
                                                                   agACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTc                           <  2:5090281/93‑1 (MQ=255)
                                                                    gACAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCa                          <  1:6515536/93‑1 (MQ=255)
                                                                     aCAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAg                         <  2:2589711/93‑1 (MQ=255)
                                                                     aCAAATAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAg                         >  1:5165857/1‑93 (MQ=255)
                                                                     aCAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAg                         >  1:4173474/1‑93 (MQ=255)
                                                                      cAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAgg                        >  1:55185/1‑93 (MQ=255)
                                                                       aaaTAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAggg                       >  1:1063340/1‑93 (MQ=255)
                                                                        aaTAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGt                      <  1:3713178/93‑1 (MQ=255)
                                                                        aaTAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGt                      >  1:4813885/1‑93 (MQ=255)
                                                                         aTAACCCGATAAAACAAGGGCGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGtt                     >  2:2635856/1‑93 (MQ=255)
                                                                          tAACCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTc                    <  2:5257443/93‑1 (MQ=255)
                                                                            aCCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAg                  >  1:320378/1‑93 (MQ=255)
                                                                            aCCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAg                  <  1:526836/93‑1 (MQ=255)
                                                                            aCCCGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAg                  <  1:2552675/93‑1 (MQ=255)
                                                                               cGATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCg               >  2:2431931/1‑93 (MQ=255)
                                                                                gATAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCgg              >  1:195177/1‑93 (MQ=255)
                                                                                 aTAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCggg             <  2:365042/93‑1 (MQ=255)
                                                                                 aTAAAACAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCggg             >  1:6257686/1‑93 (MQ=255)
                                                                                   aaaaCAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGt           >  2:68860/1‑93 (MQ=255)
                                                                                   aaaaCAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGt           >  2:314647/1‑93 (MQ=255)
                                                                                     aaCAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTgc         >  2:5616812/1‑93 (MQ=255)
                                                                                     aaCAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTgc         >  2:377408/1‑93 (MQ=255)
                                                                                     aaCAAGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTgc         >  2:78724/1‑93 (MQ=255)
                                                                                        aaGGGTGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCa      <  1:2516807/93‑1 (MQ=255)
                                                                                         aGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCAt     >  1:4932040/1‑93 (MQ=255)
                                                                                            |                                                                                            
TTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGG  >  NC_000913/3759711‑3759895

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