Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,906 A→G 5.9% I166V (ATT→GTT)  gsiB → glutathione ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913869,9060AG 5.9% 341.8 / 7.8 102I166V (ATT→GTT) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (50/46);  new base G (3/3);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

Only 100 of 109 total aligned reads displayed.

CGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAAC  >  NC_000913/869769‑870053
                                                                                                                                         |                                                                                                                                                   
cGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAAATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCa                                                                                                                                          <  1:426560/149‑1 (MQ=255)
    tcatcTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAgcgc                                                                                                                                      >  3:989961/1‑149 (MQ=255)
            aaCGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGaaaaa                                                                                                                              <  4:989961/149‑1 (MQ=255)
              cGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAata                                                                                                                            <  1:1307202/149‑1 (MQ=255)
                 tataACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATgg                                                                                                                         >  3:925796/1‑149 (MQ=255)
                          tataAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGTGACCGCGATGATTTCAcc                                                                                                                                              <  4:1426059/119‑1 (MQ=255)
                          tataAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCAcc                                                                                                                                              >  3:1426059/1‑119 (MQ=255)
                             aaGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTgg                                                                                                             >  2:1285730/1‑149 (MQ=255)
                              aGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAatat                                                                                                                           >  1:766287/1‑134 (MQ=255)
                              aGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAatat                                                                                                                           <  2:766287/134‑1 (MQ=255)
                              aGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGt                                                                                                            <  4:925796/149‑1 (MQ=255)
                                aaTATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGttt                                                                                                          >  1:632587/1‑149 (MQ=255)
                                        tAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGAt                                                                                                                >  4:126517/1‑135 (MQ=255)
                                        tAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCgg                                                                                                  >  2:439720/1‑149 (MQ=255)
                                        tAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCgg                                                                                                  >  2:758734/1‑149 (MQ=255)
                                        tAAAACGGAAGCGATCGATCCGACAACGGGAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGAt                                                                                                                <  3:126517/135‑1 (MQ=255)
                                          aaaCGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTg                                                                                                <  2:632587/149‑1 (MQ=255)
                                          aaaCGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTg                                                                                                >  1:814975/1‑149 (MQ=255)
                                           aaCGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTc                                                                                                                                                 <  3:215899/99‑1 (MQ=255)
                                           aaCGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTc                                                                                                                                                 >  4:215899/1‑99 (MQ=255)
                                                 aGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCg                                                                                         <  1:439720/149‑1 (MQ=255)
                                                  gCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAg                                                                                                                     <  3:1183237/120‑1 (MQ=255)
                                                    gatcgatcCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGAc                                                                                      >  1:1458419/1‑149 (MQ=255)
                                                     atcgatcCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCa                                                                                                                                          >  1:1227780/1‑96 (MQ=255)
                                                     atcgatcCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCa                                                                                                                                          <  2:1227780/96‑1 (MQ=255)
                                                     atcgatcCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGAcc                                                                                     >  3:1019019/1‑149 (MQ=255)
                                                     atcgatcCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGAcc                                                                                     >  1:247625/1‑149 (MQ=255)
                                                         atcCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTAt                                                                                 >  3:1070304/1‑149 (MQ=255)
                                                            cGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGaa                                                                              <  2:814975/149‑1 (MQ=255)
                                                            cGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGaa                                                                              >  3:1361531/1‑149 (MQ=255)
                                                               cAACTGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTg                                                                           <  1:1285730/149‑1 (MQ=255)
                                                                    gTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTc                                                                                                                                                 >  4:801560/1‑74 (MQ=255)
                                                                    gTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTc                                                                                                                                                 <  3:801560/74‑1 (MQ=255)
                                                                     tAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTgg                                                                                               <  4:663074/123‑1 (MQ=255)
                                                                     tAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGa                                                                         >  2:106272/1‑145 (MQ=255)
                                                                     tAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGa                                                                         <  1:106272/145‑1 (MQ=255)
                                                                                cTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGACACCTGGAATCAGa                                                           >  4:1118031/1‑148 (MQ=255)
                                                                                cTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGACACCTGGAATCAGa                                                           <  3:1118031/148‑1 (MQ=255)
                                                                                 tCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGAcc                                                         >  1:338550/1‑149 (MQ=255)
                                                                                   aaaCAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGa                                                       <  4:1361531/149‑1 (MQ=255)
                                                                                     aCAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATcc                                                                                                    <  3:1057957/102‑1 (MQ=255)
                                                                                     aCAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATcc                                                                                                    >  4:1057957/1‑102 (MQ=255)
                                                                                     aCAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGAt                                                      >  3:496052/1‑148 (MQ=255)
                                                                                       aGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGAtttt                                                   >  1:953035/1‑149 (MQ=255)
                                                                                        gCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTtg                                                  <  2:338550/149‑1 (MQ=255)
                                                                                           gTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTg                                                                   <  1:1014722/129‑1 (MQ=255)
                                                                                           gTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTg                                                                   >  2:1014722/1‑129 (MQ=255)
                                                                                              ctcAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAgg                                            <  2:247625/149‑1 (MQ=255)
                                                                                                 aGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTg                                                                           >  2:140364/1‑115 (MQ=255)
                                                                                                   cGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGaaa                                       <  2:953035/149‑1 (MQ=255)
                                                                                                    gTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTg                                                                   <  2:714893/120‑1 (MQ=255)
                                                                                                    gTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTg                                                                   >  1:714893/1‑120 (MQ=255)
                                                                                                    gTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGaaaa                                      <  2:1458419/149‑1 (MQ=255)
                                                                                                        attaATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAAtt                                  >  4:1370098/1‑149 (MQ=255)
                                                                                                         ttaATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTc                                 >  4:1171159/1‑149 (MQ=255)
                                                                                                           aaTATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGa                                                       >  2:306489/1‑125 (MQ=255)
                                                                                                           aaTATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGa                                                       <  1:306489/125‑1 (MQ=255)
                                                                                                           aaTATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTcgc                               <  3:1171159/149‑1 (MQ=255)
                                                                                                            atatTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAAc                                                                             >  1:1031291/1‑102 (MQ=255)
                                                                                                            atatTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAAc                                                                             <  2:1031291/102‑1 (MQ=255)
                                                                                                            atatTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTcgcg                              >  4:646624/1‑149 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGGTTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  3:28639/1‑110 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGGTTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  4:778551/110‑1 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGGTTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  4:28639/110‑1 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGGTTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  1:751036/1‑110 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGGTTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  3:778551/1‑110 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTcgc                               >  1:1387163/1‑147 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTcgc                               <  2:1387163/147‑1 (MQ=255)
                                                                                                             tatTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCgg                             <  4:1019019/149‑1 (MQ=255)
                                                                                                                  atGCCCATCCGGCAACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTAGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTCGTGAAGGTGAAAAAATTCGCGAGATAc                        >  4:944878/2‑149 (MQ=255)
                                                                                                                    gCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCgg                             <  2:722054/142‑1 (MQ=255)
                                                                                                                    gCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCgg                             >  1:722054/1‑142 (MQ=255)
                                                                                                                       cATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGttttt                                                                                                        <  4:336395/64‑1 (MQ=255)
                                                                                                                       cATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGttttt                                                                                                        >  3:336395/1‑64 (MQ=255)
                                                                                                                       cATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGaaaa                                      >  3:731141/1‑130 (MQ=255)
                                                                                                                       cATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGaaaa                                      <  4:731141/130‑1 (MQ=255)
                                                                                                                        aTCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  4:1251821/99‑1 (MQ=255)
                                                                                                                        aTCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAg                  >  1:1383355/1‑149 (MQ=255)
                                                                                                                        aTCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAg                  >  3:1404406/1‑149 (MQ=255)
                                                                                                                          ccGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  3:361194/1‑97 (MQ=255)
                                                                                                                          ccGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  4:361194/97‑1 (MQ=255)
                                                                                                                            ggCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAg              <  3:1370098/149‑1 (MQ=255)
                                                                                                                             gCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  1:865888/1‑94 (MQ=255)
                                                                                                                             gCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  2:865888/94‑1 (MQ=255)
                                                                                                                             gCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTcgc                               >  1:1320704/1‑131 (MQ=255)
                                                                                                                               gACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  2:1152721/92‑1 (MQ=255)
                                                                                                                               gACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  1:1152721/1‑92 (MQ=255)
                                                                                                                               gACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  3:878150/92‑1 (MQ=255)
                                                                                                                               gACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGAtt                                                     <  1:1094115/107‑1 (MQ=255)
                                                                                                                               gACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGAtt                                                     >  2:1094115/1‑107 (MQ=255)
                                                                                                                               gACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGACTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  4:878150/1‑92 (MQ=255)
                                                                                                                                  cgcgATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGc        <  4:1404406/149‑1 (MQ=255)
                                                                                                                                    cgATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGccc      >  1:368825/1‑149 (MQ=255)
                                                                                                                                     gatgatTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAggtt           >  3:703619/1‑141 (MQ=255)
                                                                                                                                     gatgatTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGAt           <  4:703619/143‑1 (MQ=255)
                                                                                                                                      atgatTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGAc                                                          >  4:500973/1‑95 (MQ=255)
                                                                                                                                      atgatTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGAc                                                          <  3:500973/95‑1 (MQ=255)
                                                                                                                                        gatTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAAc  >  4:321724/1‑149 (MQ=255)
                                                                                                                                         atTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    >  2:1390054/1‑82 (MQ=255)
                                                                                                                                         atTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCt                                                                    <  1:1390054/82‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                                   
CGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTTATTAATATTCTTGCCCATCCGGCGACCGCGATGATTTCACCGGCAGCGCTGGAAAAATATGGCAAGGAGATTGGTTTTTATCCGGTGGGAACCGGACCGTATGAACTGGATACCTGGAATCAGACCGATTTTGTGAAGGTGAAAAAATTCGCGGGTTACTGGCAGCCAGGATTGCCCAAAC  >  NC_000913/869769‑870053

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