Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	513,462	A→G	6.8%	intergenic (‑109/‑82)	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_002947	513,462	0	A	G	6.8%	214.5 / 4.7	73	intergenic (‑109/‑82)	crp/yhfA	DNA‑binding transcriptional dual regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base G (2/3);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.21e-01

GATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATGTCG  >  NC_002947/513315‑513610                                                                                                                                                    |                                                                                                                                                      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)                                        cACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGATGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGt                                                                                                              <  5:190485/149‑1 (MQ=255)                                         aCGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGt                                                                                                              <  7:116800/148‑1 (MQ=255)                                         aCGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGt                                                                                                              <  7:176567/148‑1 (MQ=255)                                         aCGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTc                                                                                                             >  2:132066/1‑149 (MQ=255)                                           gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCt                                                                                                                                    >  8:184549/1‑125 (MQ=255)                                           gAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCt                                                                                                                                    <  7:184549/125‑1 (MQ=255)                                                acTGTTCGATATGACCCATAGGCTTTTTAAGAGTGTGGCTGACGGAGAGCGCCAGTCATTTGCCGCCAGCTTATCAGACACTACGAGATCGTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTg                                                                                                      <  6:307769/148‑1 (MQ=255)                                                  tGTTCGATATGACGCATAGGCTTTTTTAGAGGGTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                                     >  1:215794/1‑148 (MQ=255)                                                    ttCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                                     <  6:271950/146‑1 (MQ=255)                                                    ttCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGa                                                                                                     >  5:271950/1‑146 (MQ=255)                                                             aCGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTcc                                                                                          <  2:215794/148‑1 (MQ=255)                                                               gCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGa                                                                                        >  6:43900/1‑148 (MQ=255)                                                                 aTAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaaaa                                                                                     >  7:65088/1‑149 (MQ=255)                                                                 aTAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaaaa                                                                                     >  3:127092/1‑149 (MQ=255)                                                                     gCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                                 <  8:65088/149‑1 (MQ=255)                                                                     gCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                                 >  7:68877/1‑149 (MQ=255)                                                                       ttttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCca                                                                                                                            >  4:225247/1‑104 (MQ=255)                                                                       ttttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCca                                                                                                                            <  3:225247/104‑1 (MQ=255)                                                                       ttttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAg                                                                                 >  5:178312/1‑147 (MQ=255)                                                                         ttttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGa                                                                               <  1:132066/147‑1 (MQ=255)                                                                           ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaa                                                                                       >  3:210076/1‑137 (MQ=255)                                                                           ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGaa                                                                                       <  4:210076/137‑1 (MQ=255)                                                                           ttAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTcc                                                                           >  2:6991/1‑149 (MQ=255)                                                                             agagGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCgg                                                                         >  1:336140/1‑149 (MQ=255)                                                                                 gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGATACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                                >  5:93856/1‑122 (MQ=255)                                                                                 gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGATACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                                <  6:93856/122‑1 (MQ=255)                                                                                 gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                                <  3:203729/122‑1 (MQ=255)                                                                                 gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACAccc                                                                                                >  4:203729/1‑122 (MQ=255)                                                                                 gTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAg                                                                     <  6:178312/149‑1 (MQ=255)                                                                                    ggCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAATACGACACATCTCGTCCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGAt                                                                   >  7:251191/1‑148 (MQ=255)                                                                                     gCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGa                                                                 <  5:43900/149‑1 (MQ=255)                                                                                          cGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCa                                                            <  2:336140/149‑1 (MQ=255)                                                                                              gagCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGt                                                                       >  7:87173/1‑134 (MQ=255)                                                                                              gagCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGt                                                                       <  8:87173/134‑1 (MQ=255)                                                                                                 cgcCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATc                                                      >  8:82021/1‑148 (MQ=255)                                                                                                   ccAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAg                                                   >  3:104668/1‑149 (MQ=255)                                                                                                    cAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGt                                                  <  4:127092/149‑1 (MQ=255)                                                                                                         aTTAGGCGCCAGATTATCAGACACTACGAGATGTTTCGGTACTAATCCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGACGGCACGTATCCAGTGGGcc                                             <  5:333288/149‑1 (MQ=255)                                                                                                             ggCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGt                                                                                                              <  7:263111/80‑1 (MQ=255)                                                                                                             ggCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGt                                                                                                              >  8:263111/1‑80 (MQ=255)                                                                                                             ggCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGt                                          <  7:289923/148‑1 (MQ=255)                                                                                                                gcCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAgg                                      >  1:130971/1‑149 (MQ=255)                                                                                                                      ttATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGt                                >  5:259101/1‑149 (MQ=255)                                                                                                                          cAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTTGTTCCCACAGCCGCACAAAGTCTAAGATGACACACTTTTCCGAAAATCAGGAGTCCGGGTAAATGAAGGCACGAAACCATTGGCCCGGTGAGGCAATgatac                             >  1:76465/1‑145 (MQ=255)                                                                                                                            gACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCg                          >  5:8068/1‑149 (MQ=255)                                                                                                                            gACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCg                          >  7:166005/1‑149 (MQ=255)                                                                                                                             acacTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCgg                         <  7:82021/149‑1 (MQ=255)                                                                                                                               acTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCg                       <  6:259101/149‑1 (MQ=255)                                                                                                                                 tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCg                                            <  1:106779/126‑1 (MQ=255)                                                                                                                                 tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCg                                            >  2:106779/1‑126 (MQ=255)                                                                                                                                 tACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCAAGTGGGCCGGTGAGGCAATGTTCCTCGGCGaa                     >  1:266389/1‑149 (MQ=255)                                                                                                                                  aCGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAAt                    >  1:142346/1‑149 (MQ=39)                                                                                                                                   cGAGATGTTTCGGTACTAATCCGACAGATCTCGTGCTT‑TTCCCACAGCCGCACTAGCTCTAAGCTGACTCCCTTTTCCGTAAATCAGGAGTCCGGGTTGATGAAGGCACGTATCGAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAAt                    <  2:76465/148‑1 (MQ=255)                                                                                                                                   cGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGaa                                                                >  7:281104/1‑104 (MQ=255)                                                                                                                                   cGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGaa                                                                <  8:281104/104‑1 (MQ=255)                                                                                                                                   cGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTc                                                                                           >  3:53126/1‑77 (MQ=18)                                                                                                                                   cGAGATGTGTCGGTTTAGTACCGAAACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTc                                                                                           <  4:53126/77‑1 (MQ=18)                                                                                                                                    gagaTGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCg                  >  7:190476/1‑149 (MQ=255)                                                                                                                                     agaTGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCgg                 <  2:130971/149‑1 (MQ=255)                                                                                                                                         gTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCAAAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAg             >  1:248887/1‑149 (MQ=255)                                                                                                                                              ggTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTc                              <  5:309080/127‑1 (MQ=255)                                                                                                                                              ggTACTAAACCGACACATCTCGTGCTT‑ATCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTTGAAGAAGGCACGTATCCACTGCGCCGGTGAGGCAATGTTc                              >  6:309080/1‑127 (MQ=255)                                                                                                                                                tACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCAc                                                           <  2:141511/96‑1 (MQ=255)                                                                                                                                                tACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCAc                                                           >  1:141511/1‑96 (MQ=255)                                                                                                                                                   tAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATgtc   >  6:207181/1‑149 (MQ=255)                                                                                                                                                   tAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACCTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCgg                 >  8:319411/1‑135 (MQ=255)                                                                                                                                                   tAAACCGACACATCTCGTGCTT‑TTCCCACACCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCgg                 <  7:319411/135‑1 (MQ=255)                                                                                                                                                    aaaCCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATgtcg  >  6:35452/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GATGTTCTTTATCTTGGCGGGTAGGGCGGAGGCAACCATCACGAAATCCTGTTCGATATGACGCATAGGCTTTTTTAGAGGTTGGCTGACGGAGAGCGCCAGTCATTTGGCGCCAGCTTATCAGACACTACGAGATGTTTCGGTACTAAACCGACACATCTCGTGCTT‑TTCCCACAGCCGCACAAGGTCTAAGCTGACACCCTTTTCCGAAAATCAGGAGTCCGGGTAGATGAAGGCACGTATCCAGTGGGCCGGTGAGGCAATGTTCCTCGGCGAATCGGGGAGCGGCCATGTCG  >  NC_002947/513315‑513610

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: ‑