Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,652,125 C→T 6.1% G75S (GGC→AGC)  gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,652,1250CT 6.1% 1023.9 / 34.4 324G75S (GGC→AGC) gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (153/151);  new base T (10/10);  total (163/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 333 total aligned reads displayed.

TGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAGGC  >  NC_002947/1651976‑1652274
                                                                                                                                                     |                                                                                                                                                     
tgcagtgTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGcc                                                                                                                                                       <  2:6334504/150‑1 (MQ=255)
  cagtgTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCca                                                                                                                                                     <  1:4395473/150‑1 (MQ=255)
            tttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAg                                                                                                                                           >  2:5882844/1‑150 (MQ=255)
            tttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAg                                                                                                                                           <  1:2300841/150‑1 (MQ=255)
             ttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAgg                                                                                                                                          <  2:652131/150‑1 (MQ=255)
             ttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAgg                                                                                                                                          <  2:651918/150‑1 (MQ=255)
             ttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAgg                                                                                                                                          <  1:6863904/150‑1 (MQ=255)
                ctggctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCTGGGCAGCAccgtgtgg                                                                                                                                       <  1:1238185/150‑7 (MQ=255)
                ctggctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCa                                                                                                                                       <  2:5499715/150‑1 (MQ=255)
                 tggctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACATGGTGCTGCTGGGCAGCAccgtgtggg                                                                                                                                      <  1:2567409/150‑8 (MQ=255)
                 tggctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCAt                                                                                                                                      >  1:1014449/1‑150 (MQ=255)
                 tggctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCAt                                                                                                                                      >  2:3685869/1‑150 (MQ=255)
                     tggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAg                                                                                                                                  >  1:6534622/1‑150 (MQ=255)
                      gggCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAgg                                                                                                                                 <  1:6311349/150‑1 (MQ=255)
                         ccGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGccc                                                                                                                              <  1:1439973/150‑1 (MQ=255)
                            aTCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAgg                                                                                                                           <  2:5856637/150‑1 (MQ=255)
                            aTCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAgg                                                                                                                           <  2:2506477/150‑1 (MQ=255)
                             tCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAggg                                                                                                                          <  1:5016784/150‑1 (MQ=255)
                              cAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGc                                                                                                                         <  1:5349148/150‑1 (MQ=255)
                                         gTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGtt                                                                                                              <  1:3047188/150‑1 (MQ=255)
                                          tGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTg                                                                                                             <  1:6043496/150‑1 (MQ=255)
                                               aCGGCAGCAGCAACTCCACGCCGATGGCGACCCAGAGAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGCACCGGGCCCACACGGGGCTGCCCAGCAGCACCAGGCCATACCGGCCCAGGGCGATCTGGCGTTGGGCca                                                                                                        <  2:3944651/150‑1 (MQ=255)
                                               aCGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCca                                                                                                        <  1:5399665/150‑1 (MQ=255)
                                               aCGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCca                                                                                                        <  1:4895413/150‑1 (MQ=255)
                                               aCGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCca                                                                                                        >  1:2415069/1‑150 (MQ=255)
                                                cGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcag                                                                                                       <  1:4821827/150‑1 (MQ=255)
                                                cGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcag                                                                                                       <  2:2668646/150‑1 (MQ=255)
                                                  gcagcagcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagca                                                                                                     <  2:480148/150‑1 (MQ=255)
                                                        gcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagcag                                                                                                 >  1:791408/1‑148 (MQ=255)
                                                        gcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagcag                                                                                                 <  2:791408/148‑1 (MQ=255)
                                                        gcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGa                                                                                               <  2:5684104/150‑1 (MQ=255)
                                                         caACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGAc                                                                                              >  1:2116454/1‑150 (MQ=255)
                                                         caACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGAc                                                                                              <  2:1145769/150‑1 (MQ=255)
                                                           acacCACGCCGATAGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGc                                                                                            >  2:291170/1‑150 (MQ=255)
                                                             accacGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGccc                                                                                          >  1:522911/1‑150 (MQ=255)
                                                                acGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGTCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGaa                                                                                       >  2:5288271/1‑150 (MQ=255)
                                                                acGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGaa                                                                                       >  1:2226600/1‑150 (MQ=255)
                                                                 cGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGaaa                                                                                      >  1:6201932/1‑150 (MQ=255)
                                                                  gCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAg                                                                                     >  2:4417240/1‑150 (MQ=255)
                                                                   ccGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGTCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                                                                                    <  1:3955292/150‑1 (MQ=255)
                                                                   ccGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGTCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                                                                                    <  1:5288271/150‑1 (MQ=255)
                                                                   ccGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAgg                                                                                    <  2:2582163/150‑1 (MQ=255)
                                                                    cGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAggg                                                                                   >  2:6807103/1‑150 (MQ=255)
                                                                        ggCGAACCAGATAGTCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCaaa                                                                               >  2:1501472/1‑150 (MQ=255)
                                                                        ggCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCaaa                                                                               <  1:4588170/150‑1 (MQ=255)
                                                                          cGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGc                                                                             <  2:4721902/150‑1 (MQ=255)
                                                                           gAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGcc                                                                            >  1:4720312/1‑150 (MQ=255)
                                                                                 gATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATAcc                                                                      >  1:4638681/1‑150 (MQ=255)
                                                                                 gATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATAcc                                                                      >  2:6245455/1‑150 (MQ=255)
                                                                                      ccagccGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGa                                                                 >  2:3876226/1‑150 (MQ=255)
                                                                                      ccagccGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGa                                                                 >  2:455591/1‑150 (MQ=255)
                                                                                         gccGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCa                                                              <  1:6729777/150‑1 (MQ=255)
                                                                                           cGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAgg                                                            <  2:2344058/150‑1 (MQ=255)
                                                                                             tGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgta                                                          <  1:4071060/150‑1 (MQ=255)
                                                                                                cACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgtagta                                                       <  1:3920896/150‑1 (MQ=255)
                                                                                                 aCCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTCGTAc                                                      >  2:5592653/1‑150 (MQ=255)
                                                                                                 aCCAGTTGCCACCAGTCCAGGCCTACGTCGAACCGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTAc                                                      >  2:6806157/1‑150 (MQ=255)
                                                                                                   cAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAg                                                    <  1:6241974/150‑1 (MQ=255)
                                                                                                    aGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGc                                                   <  2:1298988/150‑1 (MQ=255)
                                                                                                     gTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCg                                                  >  1:5925911/1‑150 (MQ=255)
                                                                                                      ttGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGAt                                                                <  1:1992820/135‑1 (MQ=255)
                                                                                                        gCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATc                                               >  1:5182711/1‑150 (MQ=255)
                                                                                                         ccaccaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATcc                                              >  1:6577538/1‑150 (MQ=255)
                                                                                                           accaGTCCAGGCCTACTTCGAACAGGGCTCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGc                                            >  1:5729020/1‑150 (MQ=255)
                                                                                                           accaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTCGTACAGCGATCCGc                                            >  1:3592720/1‑150 (MQ=255)
                                                                                                           accaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGc                                            >  2:4867307/1‑150 (MQ=255)
                                                                                                           accaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGc                                            >  1:2609791/1‑150 (MQ=255)
                                                                                                           accaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGc                                            >  1:2857012/1‑150 (MQ=255)
                                                                                                             caGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGccc                                          >  2:3311043/1‑150 (MQ=255)
                                                                                                                  cAGGCCAACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCAAACAGGCCCAGGGCGATCAGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTg                                     >  1:952768/1‑150 (MQ=255)
                                                                                                                    ggCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGc                                   >  2:5581391/1‑150 (MQ=255)
                                                                                                                     gCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCt                                  >  1:5460357/1‑150 (MQ=255)
                                                                                                                     gCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCt                                  >  2:3371454/1‑150 (MQ=255)
                                                                                                                       cTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCa                                <  2:4987779/150‑1 (MQ=255)
                                                                                                                          cTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCt                             <  1:4202897/150‑1 (MQ=255)
                                                                                                                           ttCGAACAGGGCCCCCACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                            <  2:2309414/150‑1 (MQ=255)
                                                                                                                           ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCCAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                            >  2:4294676/1‑150 (MQ=255)
                                                                                                                           ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                            >  2:1057016/1‑150 (MQ=255)
                                                                                                                           ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                            >  1:3177506/1‑150 (MQ=255)
                                                                                                                           ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                            >  2:2124572/1‑150 (MQ=255)
                                                                                                                             cGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGa                          >  1:1788976/1‑150 (MQ=255)
                                                                                                                             cGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGa                          <  1:562256/150‑1 (MQ=255)
                                                                                                                                aCAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATAc                       >  1:4328119/1‑150 (MQ=255)
                                                                                                                                 cAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATAcc                      <  1:1882806/150‑1 (MQ=255)
                                                                                                                                  aGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATAccg                     >  1:2957798/1‑150 (MQ=255)
                                                                                                                                    ggCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGcc                                                                            >  2:2851015/1‑93 (MQ=255)
                                                                                                                                       ccACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGc                >  2:1383049/1‑150 (MQ=255)
                                                                                                                                         gcctggTGCTGCTGGGCAGCACCGGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCca              >  2:5114290/5‑150 (MQ=255)
                                                                                                                                          cctggTGCTGCTGGGCAGCACCGGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcag             >  2:2227829/4‑150 (MQ=255)
                                                                                                                                           ctggTGCTGCTGGGCAGCACCGGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcagc            >  1:6409835/3‑150 (MQ=255)
                                                                                                                                           acGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcagc            >  1:3335316/1‑150 (MQ=255)
                                                                                                                                            tggTGCTGCTGGGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcagca           >  2:1968425/2‑150 (MQ=255)
                                                                                                                                            cGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcagca           <  2:1788976/150‑1 (MQ=255)
                                                                                                                                             ggTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCATAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcagcag          <  1:2767584/150‑1 (MQ=255)
                                                                                                                                              gTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAgg         <  2:2999619/150‑1 (MQ=255)
                                                                                                                                               tgctgcCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAggg        <  1:4548410/150‑1 (MQ=255)
                                                                                                                                                gctgcCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGc       <  2:6577538/150‑1 (MQ=255)
                                                                                                                                                  tgcCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCa     >  2:6419993/1‑150 (MQ=255)
                                                                                                                                                  tgcCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCa     >  2:4513938/1‑150 (MQ=255)
                                                                                                                                                    cccAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAgg   >  1:4234817/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAGGC  >  NC_002947/1651976‑1652274

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