Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,533,408 (A)6→5 100% intergenic (+28/‑356) cysK → / → ptsH cysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,533,4030A.83.3% 220.3 / 49.1 78intergenic (+23/‑361)cysK/ptsHcysteine synthase A, O‑acetylserine sulfhydrolase A subunit/phosphohistidinoprotein‑hexose phosphotransferase component of PTS system (Hpr)
Reads supporting (aligned to +/- strand):  ref base A (5/7);  new base . (36/29);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 5.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.77e-01

ACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTA  >  NC_000913/2533305‑2533501
                                                                                                  |                                                                                                  
aCCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTaa                                                                                                   >  2:1146615/1‑100 (MQ=255)
     catcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTaaaaagc                                                                                              <  1:504746/100‑3 (MQ=255)
      atcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTaaaaagca                                                                                             >  2:1499249/1‑97 (MQ=255)
      atcGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGc                                                                                             >  1:2498912/1‑100 (MQ=255)
        cGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCAc                                                                                           <  2:2437241/100‑1 (MQ=255)
           gTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCttt                                                                                       >  2:2551951/1‑100 (MQ=255)
             gAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTCATACCCGCTTGTTAAAAATGAGTAAAAACCACCttttt                                                                                     >  1:2194886/1‑100 (MQ=255)
              aGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTa                                                                                    >  2:1735149/1‑100 (MQ=255)
               gCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAg                                                                                   >  2:2274763/1‑100 (MQ=255)
               gCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATACGTAAAAAGCACCTTTTTAg                                                                                   <  1:618589/100‑1 (MQ=255)
                cGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCGCCTTTTTAgg                                                                                  <  1:2122694/100‑1 (MQ=255)
                cGTGATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTCCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAgg                                                                                  <  1:994628/100‑1 (MQ=255)
                 gTTATTTAAGCACCGCATGGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAgg                                                                                  <  1:2450754/99‑1 (MQ=255)
                     tttAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCtt                                                                             <  1:2421975/100‑1 (MQ=255)
                       tAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCtttt                                                                           >  1:1783783/1‑100 (MQ=255)
                       tAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCttt                                                                            >  1:1786438/1‑99 (MQ=255)
                       tAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTATTAGGTGCtttt                                                                           >  1:867541/1‑100 (MQ=255)
                        aaGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCttttt                                                                          <  1:2374713/100‑1 (MQ=255)
                          gCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTtg                                                                        >  1:1896779/1‑100 (MQ=255)
                           cACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTtg                                                                        <  1:2046714/100‑1 (MQ=255)
                           cACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTAAAAAGGTGCtttttacg                                                                       <  1:2098815/100‑4 (MQ=38)
                            aCCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTtgtg                                                                      >  1:1275561/1‑100 (MQ=255)
                            aCCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTAAAAAGGTGCtttttacgc                                                                      <  2:1189459/100‑5 (MQ=37)
                             ccGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAACACCTTTTTAGGTGCTTTTTTGTgg                                                                     >  1:1043205/1‑100 (MQ=255)
                              cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGc                                                                    <  2:2438027/100‑1 (MQ=255)
                              cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTAAAAAGGTGCtttttacgcat                                                                    <  2:1104465/100‑7 (MQ=37)
                              cGCATTGTTTGCCGATCTCTTCACTGAGAAAGAAGTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTAAAAAGGTGCttttttacg                                                                      <  2:79581/99‑4 (MQ=37)
                                  ttagttGCCTAGCTCTACACTTAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTAAAAAGGTGCtttttac                                                                        <  1:2386383/88‑3 (MQ=19)
                                     tttGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTc                                                             >  2:85441/1‑100 (MQ=255)
                                      ttGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCGAGCTTGTTAAAAATGCGTAAAAAGCACCTTTGAAGGTTCGTTTTTGTTGCCGGCTTCa                                                            >  1:2585381/1‑100 (MQ=255)
                                        gCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCaaa                                                          <  1:1923495/100‑1 (MQ=255)
                                        gCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCaaa                                                          >  2:816885/1‑100 (MQ=255)
                                         ccGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAAc                                                         >  2:1722559/1‑100 (MQ=255)
                                         ccGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAAc                                                         <  2:1034813/100‑1 (MQ=255)
                                          cGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACt                                                        <  2:1932498/100‑1 (MQ=255)
                                           gATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACt                                                        <  2:998425/100‑1 (MQ=255)
                                             tctcTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTc                                                     <  2:1208512/100‑1 (MQ=255)
                                              ctctTCACTGAGAAAGAATTTCAACAGTAATGCTAGCTAGTTTAAATTGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCg                                                    <  1:481647/100‑1 (MQ=255)
                                              ctctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCg                                                    <  1:726469/100‑1 (MQ=255)
                                               tctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCg                                                    >  2:563779/1‑99 (MQ=255)
                                               tctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGc                                                   >  1:395941/1‑100 (MQ=255)
                                               tctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGc                                                   >  2:1251834/1‑100 (MQ=255)
                                                ctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGcc                                                  >  2:691053/1‑100 (MQ=255)
                                                ctTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGcc                                                  <  1:1257103/100‑1 (MQ=255)
                                                   cACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGcccc                                                >  2:2109369/1‑100 (MQ=255)
                                                   cACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCcct                                               >  2:1248179/1‑100 (MQ=255)
                                                   cACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCcct                                               >  1:1856727/1‑100 (MQ=255)
                                                     cTGTGAAAGAATTGCAACAGTAATCCCAGCTTATTAAAAATGCGTAAAAAGCACCTTCTTAGGTGCTTGTTTGTGGCCTGCTTCAAACTTTCGCCcctc                                              >  1:16197/1‑99 (MQ=255)
                                                     cTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCcctcc                                             >  2:13369/1‑100 (MQ=255)
                                                     cTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCcctcc                                             >  2:1524168/1‑100 (MQ=255)
                                                       gagaAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTg                                           >  2:14728/1‑100 (MQ=255)
                                                       gagaAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTg                                           <  1:1825564/100‑1 (MQ=255)
                                                           aagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTACAAAGCACCTTATTAGGTGCTCTTGTGTGGCCTGCTTCGAACTTTCGCCCCTCCTTGCAt                                       >  1:1941354/1‑100 (MQ=255)
                                                           aagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAACGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCAt                                       >  1:801338/1‑100 (MQ=255)
                                                           aagaaTTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCAt                                       >  1:915708/1‑100 (MQ=255)
                                                                 tgaaaCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTGTTGTGGTCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTc                                 <  1:237805/97‑1 (MQ=255)
                                                                  gCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCAttactgtt                                >  2:137429/1‑94 (MQ=255)
                                                                   cAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCCCCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAg                               <  2:2321883/100‑1 (MQ=255)
                                                                   cAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGTTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAg                               >  1:1789374/1‑100 (MQ=255)
                                                                     aCAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGcc                             <  2:264870/100‑1 (MQ=255)
                                                                      cAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCt                            >  2:1350241/1‑100 (MQ=255)
                                                                        gTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTg                           <  1:131533/100‑1 (MQ=255)
                                                                        gTAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGt                          >  2:2121847/1‑100 (MQ=255)
                                                                         tAATGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTc                         <  2:1837047/100‑1 (MQ=255)
                                                                          aaTGCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCg                        >  2:563228/1‑100 (MQ=255)
                                                                          aaTGCCAGCTTGTTAAAAATGCGGAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGcc                             <  1:2057248/95‑1 (MQ=255)
                                                                           aTGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTACTGGCATTGATTCAGCCTGTCg                        >  2:2356214/1‑100 (MQ=255)
                                                                             gCCAGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGaa                     >  2:699582/1‑100 (MQ=255)
                                                                               cAGCTTGTTAAAAATGCGTAAAACGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACt                   <  2:2383557/100‑1 (MQ=255)
                                                                                atcTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACt                   <  2:1554303/98‑1 (MQ=255)
                                                                                aGCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTg                  >  2:690668/1‑100 (MQ=255)
                                                                                 gCTTGTTAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACt                   >  1:2412917/1‑98 (MQ=255)
                                                                                    tGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGAGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCAGGCATTGATTCAGCCTGTCGGAACtagat               >  1:1554294/1‑96 (MQ=255)
                                                                                      ttAAAAATGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTAttt            <  2:348944/100‑1 (MQ=255)
                                                                                         aaaaTGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCTTGGCATTGATTCAGCCTGTCGGAACTGGTCTTTAAc         >  1:2201034/1‑100 (MQ=255)
                                                                                         aaaaTGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAAc         <  2:306703/100‑1 (MQ=255)
                                                                                          aaaTGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAAcc        <  1:2551932/100‑1 (MQ=255)
                                                                                           aaTGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAAcc        <  1:1116975/100‑1 (MQ=255)
                                                                                           aaTGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCa       <  1:1532181/100‑1 (MQ=255)
                                                                                             tGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGa     >  2:1968972/1‑100 (MQ=255)
                                                                                             tGCGTAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGa     <  2:521093/100‑1 (MQ=255)
                                                                                             gccacaaaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGGCTGTCGGAACTGGTATTTAACCAg      >  2:251413/6‑100 (MQ=37)
                                                                                                acaaaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACt   >  1:762320/3‑100 (MQ=38)
                                                                                                 caaaaaaGCACCTAAAAAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTa  >  1:225797/2‑100 (MQ=38)
                                                                                                  |                                                                                                  
ACCATCATCGGGTGAGCGTTATTTAAGCACCGCATTGTTTGCCGATCTCTTCACTGAGAAAGAATTGCAACAGTAATGCCAGCTTGTTAAAAATGCGTAAAAAAGCACCTTTTTAGGTGCTTTTTTGTGGCCTGCTTCAAACTTTCGCCCCTCCTGGCATTGATTCAGCCTGTCGGAACTGGTATTTAACCAGACTA  >  NC_000913/2533305‑2533501

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