Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 513,459 A→T 6.7% intergenic (‑106/‑85) 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,4590AT 6.7% 220.1 / 5.2 74intergenic (‑106/‑85)crp/yhfADNA‑binding transcriptional dual regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/34);  new base T (2/3);  total (37/37)
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.22e-01

CGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAG  >  NC_002947/513314‑513599
                                                                                                                                                 |                                                                                                                                             
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                                                                                                                               <  2:310505/140‑1 (MQ=255)
                      tAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACAtct                                                                                                                               >  1:310505/1‑140 (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                                                                                                                         <  7:184549/125‑1 (MQ=255)
                                           gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCt                                                                                                                         >  8:184549/1‑125 (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                                                                          >  7:65088/1‑149 (MQ=255)
                                                                 aTAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaaaa                                                                          >  3:127092/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                                                                            <  4:210076/137‑1 (MQ=255)
                                                                           ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaa                                                                            >  3:210076/1‑137 (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                                                                                                   >  8:263111/1‑80 (MQ=255)
                                                                                                             ggCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGt                                                                                                   <  7:263111/80‑1 (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               >  7:166005/1‑149 (MQ=255)
                                                                                                                            gACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCg               >  5:8068/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)
                                                                                                                                              ggTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTc                   <  5:309080/127‑1 (MQ=255)
                                                                                                                                              ggTACTAAACCGACACATCTCGTGCTTATCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTTGAAGAAGGCACGTATCCACTGCGCCGGTGAGGCAATGTTc                   >  6:309080/1‑127 (MQ=255)
                                                                                                                                                tACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCAc                                                <  2:141511/96‑1 (MQ=255)
                                                                                                                                                tACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCAc                                                >  1:141511/1‑96 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CGATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTTTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAG  >  NC_002947/513314‑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: