Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,235,689 +A intergenic (+282/‑217) pgi → / → yjbE glucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,235,6851.A97.7% 296.5 / ‑0.2 87intergenic (+278/‑221)pgi/yjbEglucosephosphate isomerase/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base A (41/44);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

GCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATT  >  NC_000913/4235588‑4235779
                                                                                                  |                                                                                               
gCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaa                                                                                               >  1:1162644/1‑100 (MQ=255)
gCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaa                                                                                               <  1:1199500/100‑1 (MQ=255)
  aTTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaaaa                                                                                            >  2:1573297/1‑100 (MQ=255)
   tttttttATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaaaaa                                                                                           <  1:1244887/100‑2 (MQ=255)
    ttttttATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaaaaat                                                                                          <  1:402723/100‑3 (MQ=255)
       tttATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATCCCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaaaaatacg                                                                                       <  2:1004298/100‑6 (MQ=255)
       tttATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGaaaaatacg                                                                                       >  1:636163/1‑95 (MQ=255)
        ttATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGc                                                                                       <  1:1296919/100‑1 (MQ=255)
           tCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGa                                                                                    <  1:1152955/100‑1 (MQ=255)
            cAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGAt                                                                                   >  2:460772/1‑100 (MQ=255)
             aGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGAtt                                                                                  >  1:136890/1‑100 (MQ=255)
              gATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTa                                                                                 >  1:889228/1‑100 (MQ=255)
               aTAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTaa                                                                                <  2:563615/100‑1 (MQ=255)
                tAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAAt                                                                               <  2:1366405/100‑1 (MQ=255)
                  aaTCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTg                                                                             >  1:233132/1‑100 (MQ=255)
                      cccTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCttt                                                                         >  1:1416804/1‑100 (MQ=255)
                      cccTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCttt                                                                         <  2:749453/100‑1 (MQ=255)
                        cTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTg                                                                       <  1:87626/100‑1 (MQ=255)
                        cTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTg                                                                       >  1:1286206/1‑100 (MQ=255)
                        cTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACACAGATTAATTGCTTTTg                                                                       <  2:1037781/100‑1 (MQ=255)
                            tCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGtttt                                                                   <  1:554812/100‑1 (MQ=255)
                             cTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGttttt                                                                  >  2:1247887/1‑100 (MQ=255)
                             cTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGttttt                                                                  >  1:1268535/1‑100 (MQ=255)
                             cTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGttttt                                                                  <  2:867208/100‑1 (MQ=255)
                             cTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGttttt                                                                  >  1:480177/1‑100 (MQ=255)
                             cTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGttttt                                                                  <  1:695473/100‑1 (MQ=255)
                              tGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTtttta                                                                 <  2:860524/100‑1 (MQ=255)
                              tGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTtttta                                                                 >  1:719929/1‑100 (MQ=255)
                               gTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTttttat                                                                >  1:1047591/1‑100 (MQ=255)
                                tAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTttttatt                                                               <  2:1523724/100‑1 (MQ=255)
                                 aaTTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTttttattt                                                              <  2:1451822/100‑1 (MQ=255)
                                   tttAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTttttatttta                                                            >  1:309084/1‑100 (MQ=255)
                                   tttAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTttttatttta                                                            >  1:429622/1‑100 (MQ=255)
                                   tttAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTttttatttta                                                            <  2:269342/100‑1 (MQ=255)
                                    ttAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTaa                                                           >  1:612268/1‑100 (MQ=255)
                                     tAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAg                                                          <  2:1446982/100‑1 (MQ=255)
                                     tAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAg                                                          <  2:1559608/100‑1 (MQ=255)
                                        cGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGttt                                                       >  2:977564/1‑100 (MQ=255)
                                          gAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTAt                                                     <  1:1149794/100‑1 (MQ=255)
                                           aaaTCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATg                                                    <  2:198702/100‑1 (MQ=255)
                                                aTACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGAtttt                                               <  2:989244/100‑1 (MQ=255)
                                                aTACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGAtttt                                               <  1:974903/100‑1 (MQ=255)
                                                aTACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGAtttt                                               <  1:1219644/100‑1 (MQ=255)
                                                  aCCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTa                                             >  2:772855/1‑100 (MQ=255)
                                                  aCCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTa                                             <  2:680480/100‑1 (MQ=255)
                                                   ccGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTAt                                            >  2:1280695/1‑100 (MQ=255)
                                                   ccGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTAt                                            <  1:433085/100‑1 (MQ=255)
                                                   ccGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTAt                                            >  2:1534637/1‑100 (MQ=255)
                                                     gTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTg                                          >  1:138779/1‑100 (MQ=255)
                                                     gTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTg                                          <  2:1286197/100‑1 (MQ=255)
                                                      tGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGt                                         <  1:1573307/100‑1 (MQ=255)
                                                        aGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTa                                       <  2:558658/100‑1 (MQ=255)
                                                         ggTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTAt                                      >  2:1245132/1‑100 (MQ=255)
                                                         ggTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTAtt                                     >  2:1461244/1‑100 (MQ=255)
                                                          gTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTAtt                                     >  1:1175720/1‑100 (MQ=255)
                                                            tAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTa                                   <  1:1247139/100‑1 (MQ=255)
                                                            tAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTa                                   <  1:455971/100‑1 (MQ=255)
                                                              aTCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTaaa                                 <  1:1279785/100‑1 (MQ=255)
                                                               tCCTAAAATAGAGTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAat                                <  2:53227/100‑1 (MQ=255)
                                                                ccTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAata                               >  1:572726/1‑100 (MQ=255)
                                                                 cTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAatat                              >  1:641418/1‑100 (MQ=255)
                                                                 cTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAatat                              >  1:640255/1‑100 (MQ=255)
                                                                 cTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAatat                              >  2:1078058/1‑100 (MQ=255)
                                                                 cTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAatat                              >  2:1078075/1‑100 (MQ=255)
                                                                   aaaaTAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATaa                            <  1:795507/100‑1 (MQ=255)
                                                                    aaaTAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTTATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAg                           <  2:363947/100‑1 (MQ=255)
                                                                     aaTAGATTTTTAATCGTTGTTTATTTCGGAAAAATATTCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAg                           <  2:936201/99‑1 (MQ=255)
                                                                     aaTAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGt                          <  2:1420762/100‑1 (MQ=255)
                                                                     aaTAGATTTTTAATCGTTGTTTATTTAGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGt                          <  2:1365504/100‑1 (MQ=255)
                                                                      aTAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGtt                         >  1:759861/1‑100 (MQ=255)
                                                                       tAGATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTg                        <  2:429620/100‑1 (MQ=255)
                                                                         gATTTTTAATCGTTGTTTATTTCGGAAAAATACGTAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGaa                      <  2:829831/100‑1 (MQ=255)
                                                                         gATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGGTTTTTATTGTTATTTAAATATAAGtt                         >  2:1541800/1‑97 (MQ=255)
                                                                         gATTTTTAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGaa                      >  1:456076/1‑100 (MQ=255)
                                                                           tttttAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAAc                    >  1:80744/1‑100 (MQ=255)
                                                                            ttttAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACt                   >  2:1373951/1‑100 (MQ=255)
                                                                            ttttAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACt                   >  2:1373733/1‑100 (MQ=255)
                                                                            ttttAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACt                   >  2:804970/1‑100 (MQ=255)
                                                                            ttttAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACt                   >  2:723334/1‑100 (MQ=255)
                                                                               tAATCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTtat                >  2:281721/1‑100 (MQ=255)
                                                                                aaTCGTTGTTTAGTTCGGAAAAAGACGCAGGTTAATTGATTTTGTGTTTAATTTCAATATATGATTTTTATTGTTATTTGAATATAAGTGGaaatt                   >  2:1004526/1‑94 (MQ=255)
                                                                                 aTCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTatat              >  2:177391/1‑100 (MQ=255)
                                                                                 aTCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTatat              >  2:855247/1‑100 (MQ=255)
                                                                                  tCGTTGTTTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATAtt             >  2:1000286/1‑100 (MQ=255)
                                                                                    gttgttTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTg           <  2:480174/100‑1 (MQ=255)
                                                                                    attgttTATTTCGTAACAATAAGCAGATCAAATACTTTTGATTTTATTTTAAGTTTATTACTTTTATTGTTATTTAAATATAAGTTGAAACTtata               >  2:247838/2‑96 (MQ=255)
                                                                                     ttgttTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGa          <  2:641415/100‑1 (MQ=255)
                                                                                      tgttTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGat         <  2:1556496/100‑1 (MQ=255)
                                                                                      tgttTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGa          <  2:640252/99‑1 (MQ=255)
                                                                                       gttTATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGata        <  2:1126869/100‑1 (MQ=255)
                                                                                         ttATTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATatt      <  2:87107/100‑1 (MQ=255)
                                                                                           aTTTCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATattca    >  1:208871/1‑100 (MQ=255)
                                                                                             ttCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATattcatt  >  1:129055/1‑100 (MQ=255)
                                                                                             ttCGGAAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATattcatt  <  2:1510832/100‑1 (MQ=255)
                                                                                                  |                                                                                               
GCATTTTTTTATCAGATAAATCCCCTTGTCTGTAATTTAACGGAAATCATACCGTGAGGTTAATCCTAAAATAGATTTTTAATCGTTGTTTATTTCGGAAAATACGCAGATTAATTGCTTTTGTTTTTATTTTAAGTTTATGATTTTTATTGTTATTTAAATATAAGTTGAAACTTATATTTGATATTCATT  >  NC_000913/4235588‑4235779

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