Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,647,178 A→G 5.7% L74P (CTG→CCG) ‡ gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,647,1780AG 5.7% 438.4 / 10.1 139L74P (CTG→CCG) ‡gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (69/62);  new base G (4/4);  total (73/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

Only 100 of 147 total aligned reads displayed.

CAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAGGCCC  >  NC_002947/1647029‑1647327
                                                                                                                                                     |                                                                                                                                                     
cagtgTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCca                                                                                                                                                       >  1:856091/1‑150 (MQ=255)
 agtgTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCcag                                                                                                                                                      >  2:4472493/1‑150 (MQ=255)
  gtgtTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCcagc                                                                                                                                                     <  2:1077352/150‑1 (MQ=255)
   tgtTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCcagca                                                                                                                                                    <  1:8675/150‑1 (MQ=255)
   tgtTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCcagca                                                                                                                                                    <  1:3981332/150‑1 (MQ=255)
        gCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCAcc                                                                                                                                               >  2:4234438/1‑150 (MQ=255)
          tttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCATCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAg                                                                                                                                             <  1:3948624/150‑1 (MQ=255)
          tttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCTGGGCAGCAccgt                                                                                                                                             <  2:704854/150‑3 (MQ=255)
          tttGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAg                                                                                                                                             >  1:2715066/1‑150 (MQ=255)
           ttGCTGGCTGGGACTATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCTGGGCAGCAccgtg                                                                                                                                            <  2:4235988/150‑4 (MQ=255)
            tGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACATGGTGCTGCTGGGCAGCAccgtgt                                                                                                                                           <  1:4201685/150‑5 (MQ=255)
            tGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGc                                                                                                                                           <  2:4656188/150‑1 (MQ=255)
              ctggctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCa                                                                                                                                         <  2:1552639/150‑1 (MQ=255)
                 gctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATAc                                                                                                                                      >  2:1699060/1‑150 (MQ=255)
                  ctggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACa                                                                                                                                     <  1:4825814/150‑1 (MQ=255)
                   tggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAg                                                                                                                                    >  1:1988311/1‑150 (MQ=255)
                   tggGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAg                                                                                                                                    <  1:5075648/150‑1 (MQ=255)
                      gCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGcc                                                                                                                                 >  1:3523373/1‑150 (MQ=255)
                       ccGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGccc                                                                                                                                >  2:3233770/1‑150 (MQ=255)
                       ccGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGccc                                                                                                                                >  2:3956431/1‑150 (MQ=255)
                           tCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAggg                                                                                                                            <  2:5086597/150‑1 (MQ=255)
                           tCAGCGGGCGAAGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCACAGCAGCACCAGGCCATACAGGCCCAggg                                                                                                                            >  1:3220295/1‑150 (MQ=255)
                             aGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCg                                                                                                                          >  1:257664/1‑150 (MQ=255)
                              gCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGa                                                                                                                         >  2:332087/1‑150 (MQ=255)
                              gCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGa                                                                                                                         <  2:888723/150‑1 (MQ=255)
                               cGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGAt                                                                                                                        <  2:3523373/150‑1 (MQ=255)
                                 ggCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCt                                                                                                                      <  2:4553087/150‑1 (MQ=255)
                                      cGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGt                                                                                                                 <  1:1944906/150‑1 (MQ=255)
                                        tGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTg                                                                                                               <  1:1079535/150‑1 (MQ=255)
                                         gCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                                                                              >  2:2100517/1‑150 (MQ=255)
                                         gCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                                                                              <  2:2412257/150‑1 (MQ=255)
                                         gCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTgg                                                                                                              <  2:3590008/150‑1 (MQ=255)
                                           ccACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGc                                                                                                            >  1:1266118/1‑150 (MQ=255)
                                             aCGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCca                                                                                                          >  2:723999/1‑150 (MQ=255)
                                                    cagcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCcagcagcag                                                                                                   >  2:3424038/1‑150 (MQ=255)
                                                     agcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAgg                                                                                                  >  1:862889/1‑150 (MQ=255)
                                                      gcaACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGa                                                                                                 >  1:4224059/1‑150 (MQ=255)
                                                       caACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGAc                                                                                                >  1:3607617/1‑150 (MQ=255)
                                                         acacCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGc                                                                                              >  1:2921411/1‑150 (MQ=255)
                                                         acacCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGc                                                                                              <  2:2369417/150‑1 (MQ=255)
                                                             cacGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGa                                                                                          >  1:1100045/1‑150 (MQ=255)
                                                                      ggCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCaaa                                                                                 <  1:2621225/150‑1 (MQ=255)
                                                                      ggCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCaaa                                                                                 <  2:195871/150‑1 (MQ=255)
                                                                       gCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAg                                                                                >  2:2664456/1‑150 (MQ=255)
                                                                       gCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAg                                                                                >  2:2664455/1‑150 (MQ=255)
                                                                        cGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGc                                                                               <  2:1597641/150‑1 (MQ=255)
                                                                         gAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGcc                                                                              <  1:3803677/150‑1 (MQ=255)
                                                                          aaCCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCa                                                                             <  2:258435/150‑1 (MQ=255)
                                                                              aGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATAc                                                                         <  1:2697330/150‑1 (MQ=255)
                                                                              aGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATAc                                                                         <  2:1802075/150‑1 (MQ=255)
                                                                               gATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATAcc                                                                        >  1:4272887/1‑150 (MQ=255)
                                                                               gATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATAcc                                                                        >  1:4272886/1‑150 (MQ=255)
                                                                                  agccagccGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGc                                                                     <  2:3831375/150‑1 (MQ=255)
                                                                                   gccagccGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGCCAGCAGCAGGGCGCCCGAAAGGGCAAAGCCAAAACCGGcca                                                                   >  2:1555543/1‑148 (MQ=255)
                                                                                    ccagccGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGa                                                                   >  1:3155553/1‑150 (MQ=255)
                                                                                    ccagccGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGa                                                                   >  1:1263494/1‑150 (MQ=255)
                                                                                          gTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgt                                                             <  1:3845487/150‑1 (MQ=255)
                                                                                          gTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgt                                                             <  1:3845426/150‑1 (MQ=255)
                                                                                          gTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgt                                                             >  2:4060766/1‑150 (MQ=255)
                                                                                          gTGGCAACTGGTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgt                                                             >  2:33703/1‑150 (MQ=255)
                                                                                            ggCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgtag                                                           >  1:174620/1‑150 (MQ=255)
                                                                                              cACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGgtagta                                                         >  2:305436/1‑150 (MQ=255)
                                                                                                       ccaccaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATcc                                                >  2:634057/1‑150 (MQ=255)
                                                                                                       ccaccaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATcc                                                >  1:2832429/1‑150 (MQ=255)
                                                                                                        caccaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCCAAGCCAATACCGGCGATCAGGTAGAACAGCGATCCg                                               >  1:414950/1‑150 (MQ=255)
                                                                                                        caccaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCg                                               >  1:3865238/1‑150 (MQ=255)
                                                                                                         accaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGc                                              >  1:1189341/1‑150 (MQ=255)
                                                                                                         accaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGc                                              >  1:1695415/1‑150 (MQ=255)
                                                                                                         accaGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGc                                              >  2:2936425/1‑150 (MQ=255)
                                                                                                            aGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCa                                           <  1:2310652/150‑1 (MQ=255)
                                                                                                              tCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGc                                         >  1:1404825/1‑150 (MQ=255)
                                                                                                                        cTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCt                               >  1:2207042/1‑150 (MQ=255)
                                                                                                                         ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGCTCAGGTAGTACAGCGATCCGCCCCGCTGGCTCAGCtt                              >  2:2618958/1‑150 (MQ=255)
                                                                                                                         ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              >  1:913256/1‑150 (MQ=255)
                                                                                                                         ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              >  2:2650373/1‑150 (MQ=255)
                                                                                                                         ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              >  1:3103229/1‑150 (MQ=255)
                                                                                                                         ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              >  1:870713/1‑150 (MQ=255)
                                                                                                                         ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              <  1:172802/150‑1 (MQ=255)
                                                                                                                         ttCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              >  1:3103185/1‑150 (MQ=255)
                                                                                                                         ttCGAACAGGGACCACACGGTGCTTCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCATGACGCCAGAAATGGCAAAGCCAATACAGGAGATCAGGTAGTACAGCGATCCGCACAGCTGGCTCAGAtt                              >  2:1256933/1‑150 (MQ=255)
                                                                                                                           cGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGa                            <  1:2946863/150‑1 (MQ=255)
                                                                                                                           cGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGa                            <  2:3997835/150‑1 (MQ=255)
                                                                                                                            gAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGAt                           <  2:2207042/150‑1 (MQ=255)
                                                                                                                            gAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGAt                           >  2:470137/1‑150 (MQ=255)
                                                                                                                                    cccACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCgg                   >  2:3481182/1‑150 (MQ=255)
                                                                                                                                    cccACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCgg                   >  1:4292072/1‑150 (MQ=255)
                                                                                                                                     ccACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGc                  <  2:2764088/150‑1 (MQ=255)
                                                                                                                                     ccACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGc                  <  2:2764122/150‑1 (MQ=255)
                                                                                                                                      cacaCGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGcc                 >  1:2340735/1‑150 (MQ=255)
                                                                                                                                       acacGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              <  2:2275494/136‑1 (MQ=255)
                                                                                                                                       acacGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCtt                              >  1:2275494/1‑136 (MQ=255)
                                                                                                                                       acacGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCca                >  1:1275866/1‑150 (MQ=255)
                                                                                                                                         ctggTGCTGCTGGGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcagc              >  2:4068214/3‑150 (MQ=255)
                                                                                                                                          tggTGCTGCTGGGCAGCACCGGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGACAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCcagca             >  1:2885834/2‑150 (MQ=255)
                                                                                                                                          tggTGCTGCTGGGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCAGGCTCAGCTTGATACCGCCGGCcagca             >  1:1672618/2‑150 (MQ=255)
                                                                                                                                               ctgcCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGcc        <  1:2042391/150‑1 (MQ=255)
                                                                                                                                                 gcCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAg      <  1:4472493/150‑1 (MQ=255)
                                                                                                                                                   ccAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAGGc    >  2:5040816/1‑150 (MQ=255)
                                                                                                                                                    cagcagCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAGGcc   >  2:3263631/1‑150 (MQ=255)
                                                                                                                                                     agcagcACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAGGccc  <  2:174620/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCACCAGTTGCCACCAGTCCAGGCCTACTTCGAACAGGGCCCACACGGTGCTGCCCAGCAGCACCAGGCCATACAGGCCCAGGGCGATCTGGCGTTGGGCCAGCAGCAGGACGCCCGAAAGGGCAAAGCCAATACCGGCGATCAGGTAGTACAGCGATCCGCCCAGCTGGCTCAGCTTGATACCGCCGGCCAGCAGGGCCAGGCCC  >  NC_002947/1647029‑1647327

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: