Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,137,696 A→G 8.8% intergenic (+81/‑209) PP_2753 → / → PP_2754 ABC transporter ATP‑binding protein/OprD family outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,137,6960AG 8.8% 256.9 / 11.6 91intergenic (+81/‑209)PP_2753/PP_2754ABC transporter ATP‑binding protein/OprD family outer membrane porin
Reads supporting (aligned to +/- strand):  ref base A (41/42);  new base G (5/3);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

TGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTTG  >  NC_002947/3137555‑3137843
                                                                                                                                             |                                                                                                                                                     
tgctgcCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGt                                                                                                                                               >  2:1197651/1‑150 (MQ=255)
tgctgcCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGt                                                                                                                                               >  2:1304755/1‑150 (MQ=255)
  ctgcCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGttt                                                                                                                                             >  1:1820920/1‑150 (MQ=255)
       gAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTc                                                                                                                                        <  1:247064/150‑1 (MQ=255)
         gCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCACCGCGGCCCAGGCCGCGTTTTTTTCGc                                                                                                                                      <  1:2234963/150‑1 (MQ=255)
                gcgcgGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTGTAGGCGTCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACg                                                                                                                               >  1:1956736/1‑150 (MQ=255)
                   cgGGATGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCTCAGGCCGCGTTTTTTTCGCCCTGACGaaa                                                                                                                            <  1:361106/150‑1 (MQ=255)
                   cgGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGaaa                                                                                                                            <  2:1387604/150‑1 (MQ=255)
                     ggAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACTCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGaaa                                                                                                                            >  1:1957121/1‑148 (MQ=255)
                     ggAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACTCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGaaa                                                                                                                            <  2:1957121/148‑1 (MQ=255)
                          atatTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCtgt                                                                                                                     <  1:1886179/150‑1 (MQ=255)
                          atatTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCtgt                                                                                                                     <  1:1628114/150‑1 (MQ=255)
                           tatTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCtgtg                                                                                                                    <  2:1820920/150‑1 (MQ=255)
                               cAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCTGGGCCGCGtt                                                                                                                                              <  1:219292/120‑1 (MQ=255)
                               cAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCTGGGCCGCGtt                                                                                                                                              >  2:219292/1‑120 (MQ=255)
                                aGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGcc                                                                                                               >  2:1093068/1‑150 (MQ=255)
                                            tGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGt                                                                                                  >  2:1813818/1‑150 (MQ=255)
                                            tGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGt                                                                                                  >  1:173989/1‑150 (MQ=255)
                                                 ggcAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGc                                                                                             <  1:515229/150‑1 (MQ=255)
                                                 ggcAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGc                                                                                             >  1:465216/1‑150 (MQ=255)
                                                 ggcAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGc                                                                                             <  2:808410/150‑1 (MQ=255)
                                                   cAAGGGCTTACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAg                                                                                           <  2:1875525/150‑1 (MQ=255)
                                                   cAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGaaa                                                                                                                            <  1:502938/118‑1 (MQ=255)
                                                   cAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGaaa                                                                                                                            >  2:502938/1‑118 (MQ=255)
                                                    aaGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGa                                                                                          >  2:698079/1‑150 (MQ=255)
                                                     aGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGAt                                                                                         <  1:1304755/150‑1 (MQ=255)
                                                       ggCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTg                                                                                       >  2:611728/1‑150 (MQ=255)
                                                         ctgactgaGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCt                                                                                     <  2:2178718/150‑1 (MQ=255)
                                                          tgactgaGCGGTGGCGACGGCTGTGGAGCGGCGCGGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCtt                                                                                    >  1:638423/1‑150 (MQ=255)
                                                          tgactgaGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCtt                                                                                    <  1:1093068/150‑1 (MQ=255)
                                                            actgaGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCg                                                                                  <  1:698079/150‑1 (MQ=255)
                                                              tgaGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCa                                                                                            <  1:1657101/138‑1 (MQ=255)
                                                              tgaGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCa                                                                                            >  2:1657101/1‑138 (MQ=255)
                                                              tgaGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCa                                                                                >  1:1965385/1‑150 (MQ=255)
                                                                       gCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGca                                                                       >  2:2041849/1‑150 (MQ=255)
                                                                         gACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAgcg                                                                                                                         <  2:1365141/99‑1 (MQ=255)
                                                                         gACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAgcg                                                                                                                         >  1:1365141/1‑99 (MQ=255)
                                                                          aCGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGcatca                                                                    >  1:978309/1‑150 (MQ=255)
                                                                              cTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGAGCGCCGAGATGGGCAGATTGCTTCGCAGACGGAGCATCATAgt                                                                >  2:2372236/1‑150 (MQ=255)
                                                                               tgtgGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAg                                                                         <  2:1917971/140‑1 (MQ=255)
                                                                               tgtgGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAg                                                                         >  1:1917971/1‑140 (MQ=255)
                                                                               tgtgGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAgtg                                                               >  1:1857541/1‑150 (MQ=255)
                                                                                    aGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAg                                                          <  1:2041849/150‑1 (MQ=255)
                                                                                    aGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAg                                                          >  1:309084/1‑150 (MQ=255)
                                                                                       ggcgCAGTCATACCGCCTTCTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGaa                                                      <  2:1588643/150‑1 (MQ=255)
                                                                                       ggcgAAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGaa                                                      <  1:191805/150‑1 (MQ=255)
                                                                                              tCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTc                                               <  1:611728/150‑1 (MQ=255)
                                                                                                aTACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCaa                                             >  2:605232/1‑150 (MQ=255)
                                                                                                   ccGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGa                                                           >  2:2040975/1‑134 (MQ=255)
                                                                                                   ccGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGa                                                           <  1:2040975/134‑1 (MQ=255)
                                                                                                    cGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATa                                                                  >  1:2249255/1‑126 (MQ=255)
                                                                                                    cGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATa                                                                  >  1:1604703/1‑126 (MQ=255)
                                                                                                    cGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATa                                                                  <  2:1604703/126‑1 (MQ=255)
                                                                                                    cGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATa                                                                  <  2:2249255/126‑1 (MQ=255)
                                                                                                     gCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGATGAATTTGTTCAATTTCt                                        <  2:761526/150‑1 (MQ=255)
                                                                                                      ccTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCtt                                       <  1:520417/150‑1 (MQ=255)
                                                                                                        ttttACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAgttgtt                                                <  1:1857527/139‑1 (MQ=255)
                                                                                                        ttttACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAgttgtt                                                >  2:1857527/1‑139 (MQ=255)
                                                                                                        ttttACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGa                                     >  1:1654213/1‑150 (MQ=255)
                                                                                                         tttACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGaa                                    <  1:1813818/150‑1 (MQ=255)
                                                                                                            aCTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGTATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATaa                                 <  1:1068842/150‑1 (MQ=255)
                                                                                                              tGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACa                               >  2:1568212/1‑150 (MQ=255)
                                                                                                                aGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAAt                             >  1:2090321/1‑150 (MQ=255)
                                                                                                                     gTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGttt                        <  1:605232/150‑1 (MQ=255)
                                                                                                                           ttCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGGGCCATGGTGGCAATTGCTTGAATAACAATGGttttttttt                  <  1:2372236/150‑1 (MQ=255)
                                                                                                                           ttCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGttttttttt                  >  2:2128106/1‑150 (MQ=255)
                                                                                                                             ccATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGc                                                                                             >  1:1162038/1‑74 (MQ=255)
                                                                                                                             ccATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGc                                                                                             <  2:1162038/74‑1 (MQ=255)
                                                                                                                             ccATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAg                                                     <  1:1592146/113‑1 (MQ=255)
                                                                                                                             ccATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAg                                                     >  2:1592146/1‑113 (MQ=255)
                                                                                                                             ccATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTttgt                >  2:395209/1‑150 (MQ=255)
                                                                                                                              cATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGcatca                                                                    >  1:251876/1‑98 (MQ=255)
                                                                                                                              cATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGcatca                                                                    <  2:251876/98‑1 (MQ=255)
                                                                                                                               aTTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTttgttt              >  2:2094651/1‑149 (MQ=255)
                                                                                                                                 tCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCtt                                       >  2:1498669/1‑123 (MQ=255)
                                                                                                                                 tCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCtt                                       <  1:1498669/123‑1 (MQ=255)
                                                                                                                                   aaCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCAt          >  1:1806713/1‑150 (MQ=255)
                                                                                                                                    aCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGa                                                           <  2:1238628/101‑1 (MQ=255)
                                                                                                                                    aCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGa                                                           >  1:1238628/1‑101 (MQ=255)
                                                                                                                                    aCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGa                                                       >  2:1696839/1‑104 (MQ=255)
                                                                                                                                    aCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGa                                                       <  1:1696839/104‑1 (MQ=255)
                                                                                                                                    aCGCGGCCTGGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTtttttgtcctt         >  2:2294900/1‑145 (MQ=255)
                                                                                                                                     cgcgGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATgg        <  2:2090321/150‑1 (MQ=255)
                                                                                                                                      gcgGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGgt       >  2:597584/1‑150 (MQ=255)
                                                                                                                                         gCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGgttgt    <  2:1965385/150‑1 (MQ=255)
                                                                                                                                          ccaaGGCCGTGTTCTTTGCGCCCGGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCGGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGgttgtt   <  1:410170/147‑1 (MQ=255)
                                                                                                                                           ccAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATaa                                 >  1:644902/1‑119 (MQ=255)
                                                                                                                                           ccAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATaa                                 <  2:644902/119‑1 (MQ=255)
                                                                                                                                           ccAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTttgtt               <  1:1709464/137‑1 (MQ=255)
                                                                                                                                           ccAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTttgtt               >  2:1709464/1‑137 (MQ=255)
                                                                                                                                           ccAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGttgttg  <  1:1041251/150‑1 (MQ=255)
                                                                                                                                           ccAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGttgttg  >  1:2094598/1‑150 (MQ=255)
                                                                                                                                             |                                                                                                                                                     
TGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTTG  >  NC_002947/3137555‑3137843

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