Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 513,453 T→G 6.6% intergenic (‑100/‑91) crp ← / → yhfA DNA‑binding transcriptional dual regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947513,4530TG 6.6% 223.6 / 7.1 76intergenic (‑100/‑91)crp/yhfADNA‑binding transcriptional dual regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/35);  new base G (2/3);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

CAGCTTGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAG  >  NC_002947/513306‑513599
                                                                                                                                                   |                                                                                                                                                   
cagcTTGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGTCGCCAGCTTATCAGACACTACGAGATGttt                                                                                                                                                    <  5:284139/149‑1 (MQ=255)
  gcTTGTCGATGTTCTTTATCTTGTCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCg                                                                                                                                                  <  4:248878/149‑1 (MQ=255)
  gcTTGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTct                                                                                                                                                  >  6:190485/1‑148 (MQ=255)
  gcTTGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCg                                                                                                                                                  <  5:248315/149‑1 (MQ=255)
   cTTGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCgg                                                                                                                                                 <  7:53848/149‑1 (MQ=255)
    ttGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAAGCATCACGATATCGTGTTAGATTTGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCGAGCTTATCAGACACTACGATATGTATCgg                                                                                                                                                 >  5:307769/1‑148 (MQ=255)
    ttGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATACGCTTTTTTAGAGGTAGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACTAGATGTTTCgg                                                                                                                                                 >  6:333288/1‑148 (MQ=255)
        cGATGTTCTTTATCTTGCAGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGAGCCAGCTTATCAGACACTACGAGATGTTTCGGTACt                                                                                                                                             >  8:116800/1‑148 (MQ=255)
         gATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTa                                                                                                                                            <  7:65791/148‑1 (MQ=255)
                    tCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACAt                                                                                                                                 <  2:144624/148‑1 (MQ=255)
                      ttGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTGTCGGTTTAGTACCGAAACAtctc                                                                                                                              <  6:300344/149‑1 (MQ=255)
                        ggCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTTGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTGTCGGTTTAGTACCGAAACAtctc                                                                                                                              <  3:171781/147‑1 (MQ=255)
                        ggCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACAtc                                                                                                                                <  6:335319/145‑1 (MQ=255)
                              tAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACAtct                                                                                                                               >  1:310505/1‑140 (MQ=255)
                              tAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACAtct                                                                                                                               <  2:310505/140‑1 (MQ=255)
                                  gCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTAc                                                                                                                                              <  3:307598/121‑1 (MQ=255)
                                  gCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTAc                                                                                                                                              >  4:307598/1‑121 (MQ=255)
                                                cACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGATGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                   <  5:190485/149‑1 (MQ=255)
                                                 aCGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                   <  7:116800/148‑1 (MQ=255)
                                                 aCGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                   <  7:176567/148‑1 (MQ=255)
                                                 aCGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTc                                                                                                  >  2:132066/1‑149 (MQ=255)
                                                   gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCt                                                                                                                         >  8:184549/1‑125 (MQ=255)
                                                   gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCt                                                                                                                         <  7:184549/125‑1 (MQ=255)
                                                        acTGTTCGATATGACCCATAGGCTTTTTAAGAGTGTGGCTGACGGAGAGCGCCAGTCATTTGCCGCCAGCTTATCAGACACTACGAGATCGTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTg                                                                                           <  6:307769/148‑1 (MQ=255)
                                                          tGTTCGATATGACGCATAGGCTTTTTTAGAGGGTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                          >  1:215794/1‑148 (MQ=255)
                                                            ttCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                          <  6:271950/146‑1 (MQ=255)
                                                            ttCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                          >  5:271950/1‑146 (MQ=255)
                                                                     aCGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTcc                                                                               <  2:215794/148‑1 (MQ=255)
                                                                       gCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGa                                                                             >  6:43900/1‑148 (MQ=255)
                                                                         aTAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaaaa                                                                          >  3:127092/1‑149 (MQ=255)
                                                                         aTAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaaaa                                                                          >  7:65088/1‑149 (MQ=255)
                                                                             gCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                      <  8:65088/149‑1 (MQ=255)
                                                                             gCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                      >  7:68877/1‑149 (MQ=255)
                                                                               ttttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCca                                                                                                                 >  4:225247/1‑104 (MQ=255)
                                                                               ttttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCca                                                                                                                 <  3:225247/104‑1 (MQ=255)
                                                                               ttttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                      >  5:178312/1‑147 (MQ=255)
                                                                                 ttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGa                                                                    <  1:132066/147‑1 (MQ=255)
                                                                                   ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaa                                                                            >  3:210076/1‑137 (MQ=255)
                                                                                   ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaa                                                                            <  4:210076/137‑1 (MQ=255)
                                                                                   ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTcc                                                                >  2:6991/1‑149 (MQ=255)
                                                                                     agagGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCgg                                                              >  1:336140/1‑149 (MQ=255)
                                                                                         gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGATACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                     <  6:93856/122‑1 (MQ=255)
                                                                                         gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGATACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                     >  5:93856/1‑122 (MQ=255)
                                                                                         gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                     >  4:203729/1‑122 (MQ=255)
                                                                                         gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                     <  3:203729/122‑1 (MQ=255)
                                                                                         gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAg                                                          <  6:178312/149‑1 (MQ=255)
                                                                                            ggCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAATACGACACATCTCGTCCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGAt                                                        >  7:251191/1‑148 (MQ=255)
                                                                                             gCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGa                                                      <  5:43900/149‑1 (MQ=255)
                                                                                                  cGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCa                                                 <  2:336140/149‑1 (MQ=255)
                                                                                                      gagCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGt                                                            >  7:87173/1‑134 (MQ=255)
                                                                                                      gagCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGt                                                            <  8:87173/134‑1 (MQ=255)
                                                                                                         cgcCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATc                                           >  8:82021/1‑148 (MQ=255)
                                                                                                           ccAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAg                                        >  3:104668/1‑149 (MQ=255)
                                                                                                            cAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGt                                       <  4:127092/149‑1 (MQ=255)
                                                                                                                 aTTAGGCGCCAGATTATCAGACACTACGAGATGTTTCGGTACTAATCCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGACGGCACGTATCCAGTGGGcc                                  <  5:333288/149‑1 (MQ=255)
                                                                                                                     ggCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                   <  7:263111/80‑1 (MQ=255)
                                                                                                                     ggCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                   >  8:263111/1‑80 (MQ=255)
                                                                                                                     ggCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGt                               <  7:289923/148‑1 (MQ=255)
                                                                                                                        gcCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAgg                           >  1:130971/1‑149 (MQ=255)
                                                                                                                              ttATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGt                     >  5:259101/1‑149 (MQ=255)
                                                                                                                                  cAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTGTTCCCACAGCCGCACAAAGTCTAAGATGACACACTTTTCCGAAAATCAGGAGTCCGGGTAAATGAAGGCACGAAACCATTGGCCCGGTGAGGCAATgatac                  >  1:76465/1‑145 (MQ=255)
                                                                                                                                    gACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCg               >  5:8068/1‑149 (MQ=255)
                                                                                                                                    gACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCg               >  7:166005/1‑149 (MQ=255)
                                                                                                                                     acacTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCgg              <  7:82021/149‑1 (MQ=255)
                                                                                                                                       acTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCg            <  6:259101/149‑1 (MQ=255)
                                                                                                                                         tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCg                                 <  1:106779/126‑1 (MQ=255)
                                                                                                                                         tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCg                                 >  2:106779/1‑126 (MQ=255)
                                                                                                                                         tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCAAGTGGGCCGGTGAGGCAATGTTCCTCGGCGaa          >  1:266389/1‑149 (MQ=255)
                                                                                                                                          aCGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAAt         >  1:142346/1‑149 (MQ=39)
                                                                                                                                           cGAGATGTTTCGGTACTAATCCGACAGATCTCGTGCTTTTCCCACAGCCGCACTAGCTCTAAGCTGACTCCCTTTTCCGTAAATCAGGAGTCCGGGTTGATGAAGGCACGTATCGAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAAt         <  2:76465/148‑1 (MQ=255)
                                                                                                                                           cGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGaa                                                     >  7:281104/1‑104 (MQ=255)
                                                                                                                                           cGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGaa                                                     <  8:281104/104‑1 (MQ=255)
                                                                                                                                           cGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTc                                                                                >  3:53126/1‑77 (MQ=18)
                                                                                                                                           cGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTc                                                                                <  4:53126/77‑1 (MQ=18)
                                                                                                                                            gagaTGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCg       >  7:190476/1‑149 (MQ=255)
                                                                                                                                             agaTGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCgg      <  2:130971/149‑1 (MQ=255)
                                                                                                                                                 gTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCAAAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAg  >  1:248887/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGCTTGTCGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAG  >  NC_002947/513306‑513599

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