Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,899,288 G→A R108W (CGG→TGG)  yoaE ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,899,2880GA99.7% 1081.4 / ‑5.5 303R108W (CGG→TGG) yoaEputative transport protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (131/171);  total (131/172)
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.82e-01

Only 100 of 306 total aligned reads displayed.

CGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGA  >  NC_000913/1899198‑1899377
                                                                                          |                                                                                         
cGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCAtt                                                                                         <  1:8781780/93‑1 (MQ=255)
   cGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCAt                                                                                      <  2:7156086/93‑1 (MQ=255)
   cGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCAt                                                                                      >  2:2702461/1‑93 (MQ=255)
     tCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGc                                                                                    >  1:7156786/1‑93 (MQ=255)
     tCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGc                                                                                    <  2:6013332/93‑1 (MQ=255)
      cAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCa                                                                                   <  1:6434887/93‑1 (MQ=255)
      cAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCAAGCCATTCATGCa                                                                                   >  1:2512484/1‑93 (MQ=255)
        aGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGt                                                                                 >  1:1050065/1‑93 (MQ=255)
        aGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGt                                                                                 >  1:6297332/1‑93 (MQ=255)
        aGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGt                                                                                 <  2:2489972/93‑1 (MQ=255)
           aTGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCGGCCATTCATGCAGTTCg                                                                              <  2:5344130/93‑1 (MQ=255)
           aTGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCg                                                                              >  1:3582401/1‑93 (MQ=255)
             gACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGgt                                                                            >  1:1906345/1‑93 (MQ=255)
             gACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGgt                                                                            <  1:7990872/93‑1 (MQ=255)
                aaTCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGgttgt                                                                         >  2:569903/1‑93 (MQ=255)
                 aTCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGgttgtt                                                                        <  1:4974936/93‑1 (MQ=255)
                 aTCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGgttgtt                                                                        >  1:633945/1‑93 (MQ=255)
                  tCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGttgttg                                                                       <  2:4443440/93‑1 (MQ=255)
                   cTGTGTGACGACCACCCAGAACCTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGc                                                                      <  2:4068293/93‑1 (MQ=255)
                    tgtgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCt                                                                     <  2:1780893/93‑1 (MQ=255)
                    tgtgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCt                                                                     <  2:7346932/93‑1 (MQ=255)
                    tgtgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCt                                                                     >  1:115613/1‑93 (MQ=255)
                      tgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCttt                                                                   <  2:5891032/93‑1 (MQ=255)
                      tgtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCttt                                                                   <  2:2236588/93‑1 (MQ=255)
                       gtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCCTGCAGTTCGGTTGTTGCTTTg                                                                  >  2:6066322/1‑93 (MQ=255)
                       gtgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTg                                                                  >  1:6741070/1‑93 (MQ=255)
                       gtgACGACAACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTg                                                                  <  1:5303183/93‑1 (MQ=255)
                        tgACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGa                                                                 <  1:5584443/93‑1 (MQ=255)
                          acgacCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAAc                                                               >  2:8632715/1‑93 (MQ=255)
                           cgacCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCTAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACa                                                              <  2:5496662/93‑1 (MQ=255)
                           cgacCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACa                                                              <  2:5558623/93‑1 (MQ=255)
                            gacCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAg                                                             >  2:1606487/1‑93 (MQ=255)
                            gacCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAg                                                             <  2:1315605/93‑1 (MQ=255)
                            gacCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAg                                                             >  1:4329858/1‑93 (MQ=255)
                             accaccCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGc                                                            <  1:5874901/93‑1 (MQ=255)
                               caccCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTCTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCaa                                                          >  1:7283186/1‑93 (MQ=255)
                                accCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAg                                                         >  1:8592675/1‑93 (MQ=255)
                                accCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAg                                                         <  2:6267570/93‑1 (MQ=255)
                                accCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAg                                                         >  1:8638231/1‑93 (MQ=255)
                                accCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAg                                                         >  1:8122496/1‑93 (MQ=255)
                                accCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAg                                                         >  2:2523677/1‑93 (MQ=255)
                                 cccAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGa                                                        >  2:8536219/1‑93 (MQ=255)
                                 cccAGAAACTGGCGTAGCCTTGACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGa                                                        <  2:3924454/93‑1 (MQ=255)
                                  ccAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGaa                                                       >  1:2606088/1‑93 (MQ=255)
                                    aGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAata                                                     >  1:7453260/1‑93 (MQ=255)
                                     gAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAatat                                                    <  1:3651125/93‑1 (MQ=255)
                                     gAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAatat                                                    <  2:148366/93‑1 (MQ=255)
                                      aaaCTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATc                                                   >  1:5493317/1‑93 (MQ=255)
                                       aaCTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcc                                                  <  1:237103/93‑1 (MQ=255)
                                       aaCTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcc                                                  >  1:1371867/1‑93 (MQ=255)
                                        aCTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATccc                                                 >  1:2161317/1‑93 (MQ=255)
                                         cTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATcccc                                                <  1:6083358/93‑1 (MQ=255)
                                               tAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGaac                                          >  1:303148/1‑93 (MQ=255)
                                               tAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGaac                                          <  2:1189557/93‑1 (MQ=255)
                                                aGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGaaca                                         <  1:3796855/93‑1 (MQ=255)
                                                 gCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGaacaa                                        <  1:2720895/93‑1 (MQ=255)
                                                  ccTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGaacaac                                       <  1:5069896/93‑1 (MQ=255)
                                                  ccTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGaacaac                                       >  1:7814060/1‑93 (MQ=255)
                                                   cTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAacaaca                                      <  1:7603670/93‑1 (MQ=255)
                                                    tttACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACAt                                     >  2:6562286/1‑93 (MQ=255)
                                                    tttACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACAt                                     >  2:717190/1‑93 (MQ=255)
                                                      tACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATaa                                   <  1:1713511/93‑1 (MQ=255)
                                                       aCCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAAt                                  <  1:3277491/93‑1 (MQ=255)
                                                        ccGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATc                                 <  2:3522702/93‑1 (MQ=255)
                                                        ccGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATc                                 <  1:5775550/93‑1 (MQ=255)
                                                        ccGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATc                                 >  1:1925749/1‑93 (MQ=255)
                                                         cGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCa                                <  2:3401287/93‑1 (MQ=255)
                                                         cGTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCa                                >  1:656601/1‑93 (MQ=255)
                                                          gTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAg                               <  2:2161317/93‑1 (MQ=255)
                                                          gTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAg                               >  2:1896309/1‑93 (MQ=255)
                                                          gTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAg                               <  2:7448821/93‑1 (MQ=255)
                                                          gTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAg                               >  1:3415488/1‑93 (MQ=255)
                                                          gTGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAg                               <  2:4985398/93‑1 (MQ=255)
                                                           tGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAgg                              <  1:3110738/93‑1 (MQ=255)
                                                           tGGCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAgg                              <  1:6162325/93‑1 (MQ=255)
                                                            ggCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGt                             <  1:388028/93‑1 (MQ=255)
                                                            ggCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGt                             <  2:8232216/93‑1 (MQ=255)
                                                             gCCGGAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTc                            <  1:8465031/93‑1 (MQ=255)
                                                                ggAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGCACAACATAATCAGGTcgcg                         >  2:1099431/1‑93 (MQ=255)
                                                                ggAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTcgcg                         <  1:4883269/93‑1 (MQ=255)
                                                                 gAATCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTcgcgc                        >  1:421823/1‑92 (MQ=255)
                                                                  aaTCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTc                       >  1:3955510/1‑93 (MQ=255)
                                                                   aTCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACa                                                              >  2:6864911/1‑53 (MQ=255)
                                                                    tCATGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCg                     <  1:1278704/93‑1 (MQ=255)
                                                                      aTGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCgga                   >  2:4505768/1‑93 (MQ=255)
                                                                      aTGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCgga                   >  1:8265928/1‑93 (MQ=255)
                                                                      aTGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCgga                   <  2:7938995/93‑1 (MQ=255)
                                                                      aTGATCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCgga                   >  2:6919071/1‑93 (MQ=255)
                                                                          tCGCGGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCggagaag               >  2:3716346/1‑93 (MQ=255)
                                                                            gcgGTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCggagaagga             <  1:468910/93‑1 (MQ=255)
                                                                               gTTTTCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCggagaaggagaa          >  1:1011459/1‑93 (MQ=255)
                                                                                ttttCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGaaa         >  2:5252588/1‑93 (MQ=255)
                                                                                 tttCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAAt        <  2:384105/93‑1 (MQ=255)
                                                                                   tCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATcc      >  1:4548139/1‑93 (MQ=255)
                                                                                   tCCAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATcc      >  2:8387450/1‑93 (MQ=255)
                                                                                    ccAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCa     >  1:5705185/1‑93 (MQ=255)
                                                                                     cAGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCAt    <  1:8477022/93‑1 (MQ=255)
                                                                                      aGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATGCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATg   >  2:2958530/1‑93 (MQ=255)
                                                                                      aGCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATg   <  2:321992/93‑1 (MQ=255)
                                                                                       gCCATTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCTTCCGGAGAAGGAGAAATCCATGa  >  2:4760287/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGGCGTCAAGGATGACAATCTGTGTGACGACCACCCAGAAACTGGCGTAGCCTTTACCGTGGCCGGAATCATGATCGCGGTTTTCCAGCCGTTCATGCAGTTCGGTTGTTGCTTTGAACAGCAAGAATATCCCCCCGAACAACATAATCAGGTCGCGTCCGGAGAAGGAGAAATCCATGA  >  NC_000913/1899198‑1899377

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: