Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 869,724 A→C 5.7% N105T (AAC→ACC)  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,7240AC 5.7% 249.4 / 6.0 87N105T (AAC→ACC) gsiBglutathione ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (44/38);  new base C (2/3);  total (46/41)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

TCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869579‑869872
                                                                                                                                                 |                                                                                                                                                    
tCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcgc                                                                                                                                                   >  7:420185/1‑149 (MQ=255)
   ttttACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcgccgc                                                                                                                                                >  3:244168/1‑149 (MQ=255)
   ttttACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAAcgccgc                                                                                                                                                >  2:339416/1‑149 (MQ=255)
        ccAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGcggcgg                                                                                                                                           >  4:125471/1‑149 (MQ=255)
         cAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGt                                                                                                                                          >  2:321680/1‑149 (MQ=255)
             ggCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGaaa                                                                                                                                      <  8:175345/149‑1 (MQ=255)
               cTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGGTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGaaatc                                                                                                                                    <  2:437608/149‑3 (MQ=255)
               cTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGc                                                                                                                                    <  1:339416/149‑1 (MQ=255)
                          ggATAAAGAGATGAAACTGAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGaaa                                                                                                                                      <  7:55015/135‑1 (MQ=255)
                          ggATAAAGAGATGAAACTGAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGaaa                                                                                                                                      >  8:55015/1‑135 (MQ=255)
                           gATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGAc                                                                                                                          <  8:420185/147‑1 (MQ=255)
                           gATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGAc                                                                                                                          <  2:104253/147‑1 (MQ=255)
                            aTAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCt                                                                                                                              <  4:276765/142‑1 (MQ=255)
                            aTAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCt                                                                                                                              >  3:276765/1‑142 (MQ=255)
                            aTAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCgg                                                                                                                       >  5:84743/1‑149 (MQ=255)
                            aTAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCgg                                                                                                                       <  7:313131/149‑1 (MQ=255)
                             tAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCggg                                                                                                                      <  7:151310/149‑1 (MQ=255)
                               aaGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATc                                                                                                                               >  2:67781/1‑138 (MQ=255)
                               aaGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATc                                                                                                                               <  1:67781/138‑1 (MQ=255)
                                agagATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCa                                                                                                                   <  4:414241/149‑1 (MQ=255)
                                     tGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGAt                                                                                                              >  2:408106/1‑149 (MQ=255)
                                      gAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATc                                                                                                             >  1:348571/1‑149 (MQ=255)
                                        aaCTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCg                                                                                                           >  6:91740/1‑149 (MQ=255)
                                        aaCTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCg                                                                                                           >  5:344850/1‑149 (MQ=255)
                                         aCTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCgg                                                                                                          <  8:149533/149‑1 (MQ=255)
                                               aaaCGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACTCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAat                                                                                                     >  4:200383/1‑148 (MQ=255)
                                                    tGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCtt                                                                                               <  7:343277/149‑1 (MQ=255)
                                                         gCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTaaa                                                                                            <  5:373949/147‑1 (MQ=255)
                                                          cGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGc                                                                                         <  4:291036/149‑1 (MQ=255)
                                                           ggAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGTCATTAAATTCCAGCATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCAGGCGAATCATCTTAAACGc                                                                                         >  5:155918/1‑148 (MQ=255)
                                                                 ttATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGCTTTCACCGCCGCGGCGTTGAAATCggtg                                                                                                                                >  6:104117/1‑100 (MQ=38)
                                                                 ttATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGCTTTCACCGCCGCGGCGTTGAAATCggtg                                                                                                                                <  5:104117/103‑4 (MQ=38)
                                                                   ataCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCt                                                                                 <  6:84743/148‑1 (MQ=255)
                                                                    taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTg                                                                                >  2:203269/1‑148 (MQ=255)
                                                                    taCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                                                                               <  3:125471/149‑1 (MQ=255)
                                                                    aaCCCTTTCCGATGACGGCATTACTTACACCCTGAAATTACGGCAAAGCATTAAAATCCAGCATGGCACCGATTTCAACGCCGCTGCGGTGAAAGCGAATCTTGACCGGGACAGCGATCCTTCGAATCATCTAAAACTATATAAccac                                                                                >  8:443964/2‑146 (MQ=255)
                                                                     aCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGAACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                                                                               >  5:237012/1‑148 (MQ=255)
                                                                     aCCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGAATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGta                                                                              >  1:187094/1‑149 (MQ=255)
                                                                      ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCggt                                                                                                                                 <  1:127357/97‑3 (MQ=39)
                                                                      ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCACCGCCGCGGCGTTGAAATCggt                                                                                                                                 >  2:127357/1‑95 (MQ=39)
                                                                      ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGta                                                                              >  7:271742/1‑148 (MQ=255)
                                                                      ccGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGta                                                                              >  1:416947/1‑148 (MQ=255)
                                                                       cGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGtat                                                                             >  1:75174/1‑148 (MQ=255)
                                                                           tCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGTCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGa                                                                         <  7:104606/148‑1 (MQ=255)
                                                                           tCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGa                                                                         <  4:244168/148‑1 (MQ=255)
                                                                              gATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAatat                                                                     >  8:141635/1‑149 (MQ=255)
                                                                               aTGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATAtt                                                                    <  8:56044/149‑1 (MQ=255)
                                                                                  aCGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGc                                                                  <  1:408106/148‑1 (MQ=255)
                                                                                     gCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCtat                                                                                      >  6:247841/1‑125 (MQ=255)
                                                                                     gCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCtat                                                                                      <  5:247841/125‑1 (MQ=255)
                                                                                      cATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa                                                             >  7:431347/1‑149 (MQ=255)
                                                                                            ttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAg                                                       >  4:287759/1‑149 (MQ=255)
                                                                                            ttacACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAg                                                       >  2:358187/1‑149 (MQ=255)
                                                                                               cacCGTGAAATTGCTGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAat                                                                       <  5:49131/130‑1 (MQ=255)
                                                                                                acCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                                                                               <  6:14006/121‑1 (MQ=255)
                                                                                                acCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGt                                                                               >  5:14006/1‑121 (MQ=255)
                                                                                                acCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgat                                                   <  8:271742/149‑1 (MQ=255)
                                                                                                  cGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa                                                             <  7:400994/137‑1 (MQ=255)
                                                                                                  cGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTaaaa                                                             >  8:400994/1‑137 (MQ=255)
                                                                                                  cGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatc                                                  >  5:309971/1‑148 (MQ=255)
                                                                                                    tGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCAAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgatcgat                                               >  7:10339/1‑149 (MQ=255)
                                                                                                      aaaTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATcc                                             <  7:349750/149‑1 (MQ=255)
                                                                                                      aaaTTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATcc                                             <  2:348571/149‑1 (MQ=255)
                                                                                                             gggAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCTATCCGGCCAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACg                                      >  5:114700/1‑149 (MQ=255)
                                                                                                                  ggCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTaaa                                 <  3:287759/149‑1 (MQ=255)
                                                                                                                    cATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGa                               <  6:344850/149‑1 (MQ=255)
                                                                                                                     aTTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAt                              >  1:167611/1‑149 (MQ=255)
                                                                                                                     aTTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAt                              >  1:32143/1‑149 (MQ=255)
                                                                                                                      ttAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgat                                                   >  6:278084/1‑127 (MQ=255)
                                                                                                                      ttAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGcgat                                                   <  5:278084/127‑1 (MQ=255)
                                                                                                                      ttAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAt                              >  6:73158/1‑148 (MQ=255)
                                                                                                                      ttAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTAGAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGAtt                             >  8:429585/1‑149 (MQ=255)
                                                                                                                         aaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAcc                          <  5:73158/149‑1 (MQ=255)
                                                                                                                         aaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAc                           <  1:358187/148‑1 (MQ=255)
                                                                                                                         aaTTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAAAATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTAcc                          <  5:91740/149‑1 (MQ=255)
                                                                                                                         aaTTCCAGGATGGCACCGATTTCAAAGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGACTCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAACCGGAAGCGATCGATCCGACACCGGTAAAGATTAcc                          <  3:6196/149‑1 (MQ=255)
                                                                                                                                      cACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGaa                                                        >  5:85706/1‑106 (MQ=255)
                                                                                                                                      cACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGaa                                                        <  6:85706/106‑1 (MQ=255)
                                                                                                                                      cACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGt             >  5:181322/1‑149 (MQ=255)
                                                                                                                                      cACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGt             >  8:365299/1‑149 (MQ=255)
                                                                                                                                      cACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGt             >  4:413130/1‑149 (MQ=255)
                                                                                                                                       aCCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGtt            >  5:139244/1‑149 (MQ=255)
                                                                                                                                       aCCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGt             >  3:426040/1‑148 (MQ=255)
                                                                                                                                       aCCGATATCAACGCCGCGGCGGTGAAAGCGAATATCCACCGGGCCACCGATCCGGCGAATCAAATTAAACGCTATAACCTCTAAAAGAAAATTGGTAAAAAGCAAGCGATCCCAACGACAACGGTAAACATTACCCTCAAACAGCCGt             >  6:78714/1‑148 (MQ=255)
                                                                                                                                         cGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTtct          <  6:181322/149‑1 (MQ=255)
                                                                                                                                           aTTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCa        <  3:413130/149‑1 (MQ=255)
                                                                                                                                           aTTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCa        <  6:237012/149‑1 (MQ=255)
                                                                                                                                                aaCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAAc                   >  7:37430/1‑133 (MQ=255)
                                                                                                                                                aaCGCCGAGGCGGTGTAAGCGAATCTGGACCGGTCCATCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACGCTCAAAc                   <  8:37430/133‑1 (MQ=255)
                                                                                                                                                 aCGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGttt  <  1:321680/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TCGTTTTACCAGGGGCTGTTCGGTCTGGATAAAGAGATGAAACTGAAAAACGTGCTGGCGGAGAGTTATACCGTTTCCGATGACGGCATTACTTACACCGTGAAATTGCGGGAAGGCATTAAATTCCAGGATGGCACCGATTTCAACGCCGCGGCGGTGAAAGCGAATCTGGACCGGGCCAGCGATCCGGCGAATCATCTTAAACGCTATAACCTGTATAAGAATATTGCTAAAACGGAAGCGATCGATCCGACAACGGTAAAGATTACCCTCAAACAGCCGTTCTCAGCGTTT  >  NC_000913/869579‑869872

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