Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,558,365 A→G intergenic (‑104/‑134) yjiH ← / → b4331 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,558,3650AG95.5% 815.4 / 18.3 245intergenic (‑104/‑134)yjiH/b4331hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (6/5);  new base G (129/105);  total (135/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 249 total aligned reads displayed.

TTGTGTGTTGCTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTATGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGTGACAA  >  NC_000913/4558275‑4558454
                                                                                          |                                                                                         
 tgtgtgTTGCTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGc                                                                                        >  1:3550628/1‑93 (MQ=255)
 tgtgtgTTGCTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGCGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTATGc                                                                                        >  1:772073/1‑93 (MQ=255)
       ttGCTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGttttt                                                                                  >  2:4200405/1‑93 (MQ=255)
       ttGCTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGttttt                                                                                  <  1:221805/93‑1 (MQ=255)
        tGCTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGtttttt                                                                                 <  2:5838518/93‑1 (MQ=255)
          cTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGc                                                                               >  2:933352/1‑93 (MQ=255)
          cTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGc                                                                               >  1:4004913/1‑93 (MQ=255)
          cTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTCTGc                                                                               <  1:2089325/93‑1 (MQ=255)
           ttGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCa                                                                              <  1:2392638/93‑1 (MQ=255)
              aGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAg                                                                           <  2:55852/93‑1 (MQ=255)
                ccAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCg                                                                         <  1:4933558/93‑1 (MQ=255)
                 cAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCgg                                                                        >  2:5377791/1‑93 (MQ=255)
                  aGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGa                                                                       >  2:3308412/1‑93 (MQ=255)
                  aGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGa                                                                       >  2:1345460/1‑93 (MQ=255)
                  aGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGa                                                                       >  1:4758056/1‑93 (MQ=255)
                  aGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGa                                                                       >  1:6127531/1‑93 (MQ=255)
                  aGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTATGCGTTTTTTGCATAGCGGa                                                                       >  2:5511893/1‑93 (MQ=255)
                   ggAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGaa                                                                      >  2:4829672/1‑93 (MQ=255)
                         ttCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTa                                                                <  1:2968857/93‑1 (MQ=255)
                         ttCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTa                                                                >  2:4497743/1‑93 (MQ=255)
                          tCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTAc                                                               >  1:2821750/1‑93 (MQ=255)
                          tCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTAc                                                               <  2:1902461/93‑1 (MQ=255)
                          tCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTAc                                                               <  2:5025912/93‑1 (MQ=255)
                          tCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTAc                                                               >  1:1155112/1‑93 (MQ=255)
                            cTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGc                                                             >  2:3213422/1‑93 (MQ=255)
                            cTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGc                                                             >  2:5044874/1‑93 (MQ=255)
                            cTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGc                                                             <  1:3724327/93‑1 (MQ=255)
                            cTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGc                                                             >  2:5690605/1‑93 (MQ=255)
                             tGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCt                                                            >  1:1899528/1‑93 (MQ=255)
                             tGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCt                                                            >  2:4248492/1‑93 (MQ=255)
                              ggAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCtg                                                           >  1:152721/1‑93 (MQ=255)
                               gAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCtgt                                                          <  1:3775805/93‑1 (MQ=255)
                               gAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCtgt                                                          >  2:2851318/1‑93 (MQ=255)
                                aaCTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCtgtg                                                         >  1:18/1‑93 (MQ=255)
                                     aaaaaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGaa                                                    >  2:3888143/1‑93 (MQ=255)
                                      aaaaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGaaa                                                   <  1:1030122/93‑1 (MQ=255)
                                      aaaaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGaaa                                                   >  2:479430/1‑93 (MQ=255)
                                       aaaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGaaaa                                                  <  1:5233755/93‑1 (MQ=255)
                                       aaaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGaaaa                                                  <  2:4040047/93‑1 (MQ=255)
                                        aaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAg                                                 >  2:2581286/1‑93 (MQ=255)
                                        aaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAg                                                 <  2:1875118/93‑1 (MQ=255)
                                        aaCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAg                                                 >  2:3655949/1‑93 (MQ=255)
                                         aCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAgg                                                >  2:4926033/1‑93 (MQ=255)
                                         aCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAgg                                                >  2:3100589/1‑93 (MQ=255)
                                         aCCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAgg                                                >  2:2991064/1‑93 (MQ=255)
                                          ccTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGc                                               <  1:740252/93‑1 (MQ=255)
                                          ccTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGc                                               >  1:3921725/1‑93 (MQ=255)
                                          ccTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGc                                               >  1:3865612/1‑93 (MQ=255)
                                           cTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGct                                              >  1:3957815/1‑93 (MQ=255)
                                            ttATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGctc                                             >  1:5750665/1‑93 (MQ=255)
                                             tATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGctct                                            >  1:3644275/1‑93 (MQ=255)
                                             tATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGctct                                            >  2:5124605/1‑93 (MQ=255)
                                             tATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGctct                                            >  1:5116115/1‑93 (MQ=255)
                                              aTCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTg                                           >  1:5896843/1‑93 (MQ=255)
                                              aTCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTg                                           >  2:52832/1‑93 (MQ=255)
                                                 aGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTgcgc                                        >  1:1083803/1‑93 (MQ=255)
                                                  gCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCa                                       >  2:1772002/1‑93 (MQ=255)
                                                  gCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTATGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCa                                       <  2:2354055/93‑1 (MQ=255)
                                                    cATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAgg                                     >  2:5357336/1‑93 (MQ=255)
                                                      tgtgGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGAt                                   <  2:2184845/93‑1 (MQ=255)
                                                      tgtgGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGAt                                   >  2:2243494/1‑93 (MQ=255)
                                                      tgtgGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGAt                                   >  2:4640059/1‑93 (MQ=255)
                                                       gtgGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGAtt                                  <  2:2510274/93‑1 (MQ=255)
                                                       gtgGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGAtt                                  >  1:5033981/1‑93 (MQ=255)
                                                        tgGGTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTg                                 >  1:1092877/1‑93 (MQ=255)
                                                         gggTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGa                                <  1:373544/93‑1 (MQ=255)
                                                          ggTATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAg                               <  2:3960280/93‑1 (MQ=255)
                                                            tATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGtt                             >  1:586745/1‑93 (MQ=255)
                                                            tATGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGtt                             >  2:2951757/1‑93 (MQ=255)
                                                              tGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGc                           <  2:1195143/93‑1 (MQ=255)
                                                              tGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGc                           >  2:2131543/1‑93 (MQ=255)
                                                              tGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGc                           >  1:6089593/1‑93 (MQ=255)
                                                              tGGGGAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGc                           >  2:6033834/1‑93 (MQ=255)
                                                                gggAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAg                         >  2:5989435/1‑93 (MQ=255)
                                                                 ggAAAATGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGa                        >  2:767831/1‑93 (MQ=255)
                                                                   aaaaTGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGaaa                      <  1:4235488/93‑1 (MQ=255)
                                                                   aaaaTGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGaaa                      <  2:4717736/93‑1 (MQ=255)
                                                                    aaaTGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGaaaa                     <  2:5948003/93‑1 (MQ=255)
                                                                      aTGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAAtt                   >  1:2581237/1‑93 (MQ=255)
                                                                      aTGCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAAtt                   <  2:3495203/93‑1 (MQ=255)
                                                                        gCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCa                 >  2:633461/1‑93 (MQ=255)
                                                                        gCAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCa                 <  2:1063638/93‑1 (MQ=255)
                                                                         cAAGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCAt                <  2:2914339/93‑1 (MQ=255)
                                                                          aaGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATa               >  1:5663129/1‑93 (MQ=255)
                                                                           aGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATaa              >  2:4379911/1‑93 (MQ=255)
                                                                           aGTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATaa              <  1:331718/93‑1 (MQ=255)
                                                                            gTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATaaa             <  2:4563412/93‑1 (MQ=255)
                                                                            gTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATaaa             <  1:688661/93‑1 (MQ=255)
                                                                            gTATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATaaa             <  1:1336181/93‑1 (MQ=255)
                                                                             tATGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAAc            <  2:3447320/93‑1 (MQ=255)
                                                                               tGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCa          >  2:272719/1‑93 (MQ=255)
                                                                               tGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCa          >  1:1973923/1‑93 (MQ=255)
                                                                               tGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCa          >  1:2628641/1‑93 (MQ=255)
                                                                               tGTTTCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCa          <  2:1605721/93‑1 (MQ=255)
                                                                                  ttCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGt       <  1:3270863/93‑1 (MQ=255)
                                                                                  ttCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGt       <  2:1684058/93‑1 (MQ=255)
                                                                                   tCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGTg      >  1:1563775/1‑93 (MQ=255)
                                                                                   tCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGTg      <  1:4831663/93‑1 (MQ=255)
                                                                                   tCTGGCTGTGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGTg      >  1:5394661/1‑93 (MQ=255)
                                                                                     tGGCTATGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGTGAc    <  1:5538385/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTGTGTGTTGCTTGAGCCAGGAAAATTCCTGGAACTGAAAAACCTTATCAGCCATGTGGGTATGGGGAAAATGCAAGTATGTTTCTGGCTATGCGTTTTTTGCATAGCGGAAAAATTACGCTGTGTGGAAAAGGCTCTGCGCAGGATTGAGTTGCAGAAAATTCATAAACCAAGTGACAA  >  NC_000913/4558275‑4558454

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