Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,351,660	0	G	A	33.0%	111.7 / 77.1	91	Q456* (CAG→TAG)	arcB	aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base G (29/32);  new base A (12/18);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.86e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTGGATCGTCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACGGCGTTA  >  NC_000913/3351514‑3351802                                                                                                                                                   |                                                                                                                                               gTGGATCGTCAACGTAAAGGTTGAACCTTTGCCCTGTGCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTgg                                                                                                                                               >  8:128152/1‑148 (MQ=255)    gATCGTCTACGTAAAGGATTAACCTTTGCCCTGTTCGCTGGTAACGGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGATTACCGCCATGACTGTCTTTCACCTGGTaa                                                                                                                                            <  5:47623/148‑1 (MQ=255)    gATCGTCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTaa                                                                                                                                            <  2:175073/148‑1 (MQ=255)    gATCGTCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGAAACCGTAATATCGCGGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAAt                                                                                                                                           <  5:77304/149‑1 (MQ=255)     aTCGTCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATa                                                                                                                                          <  7:228490/149‑1 (MQ=255)     aTCGTCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATa                                                                                                                                          >  2:198331/1‑149 (MQ=255)         tCAACGTAAAGGTTGAACCTTTGCGCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATg                                                                                                                                      >  3:315462/1‑149 (MQ=255)         tCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATg                                                                                                                                      >  5:338296/1‑149 (MQ=255)         tCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATg                                                                                                                                      >  7:243276/1‑149 (MQ=255)           aaCGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATgg                                                                                                                                     <  6:297980/148‑1 (MQ=255)           aaCGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGc                                                                                                                                    <  3:262785/149‑1 (MQ=255)              gTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGa                                                                                                                                  >  7:18477/1‑148 (MQ=255)                 aaGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGaaaa                                                                                                                               <  8:199521/148‑1 (MQ=255)                 aaGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGaaaa                                                                                                                               <  2:366154/148‑1 (MQ=255)                     ttGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTa                                                                                                                          <  4:297853/149‑1 (MQ=255)                       gAACCTTTGCCCTGTTCGATGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCATACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTAt                                                                                                                         <  7:204221/148‑1 (MQ=255)                          ccTTTGGCGTGTTCGCTGGTAACCGTAATTTCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCAGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACTTGGCGAAAATTTTATCCAg                                                                                                                     <  2:170433/149‑1 (MQ=255)                           cTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCGATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGc                                                                                                                    >  8:377362/1‑149 (MQ=255)                           cTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACGTGGCGAAAATTTTATCCAGc                                                                                                                    <  6:237495/149‑1 (MQ=255)                              tGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCa                                                                                                                 >  3:354753/1‑149 (MQ=255)                                cccTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATc                                                                                                               <  2:340912/149‑1 (MQ=255)                                 ccTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATcc                                                                                                              >  5:279734/1‑149 (MQ=255)                                     ttCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGc                                                                                                           <  4:124671/148‑1 (MQ=255)                                     ttCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGc                                                                                                           <  1:260433/148‑1 (MQ=255)                                     ttCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGc                                                                                                           <  3:347618/148‑1 (MQ=255)                                     ttCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCg                                                                                                          <  5:139486/149‑1 (MQ=255)                                          tGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATc                                                                                                                        <  7:272887/130‑1 (MQ=255)                                            gTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGcc                                                                                                   <  7:21581/149‑1 (MQ=255)                                            gTAACCGTAATATCGCCGCACATATTATTCGCCAGACGACGAGAAACGTCCAGCCCAATACCGGTGCCGGTTGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGc                                                                                                    <  6:365100/148‑1 (MQ=255)                                              aaCCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATc                                                                                                                        <  8:155867/126‑1 (MQ=255)                                              aaCCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATc                                                                                                                        >  7:155867/1‑126 (MQ=255)                                                 cGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATAc                                                                                              <  4:354753/149‑1 (MQ=255)                                                 cGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATAc                                                                                              >  2:222673/1‑149 (MQ=255)                                                   tAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATAc                                                                                              >  1:385661/1‑147 (MQ=255)                                                   tAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGAt                                                                                                >  3:179068/1‑145 (MQ=255)                                                   tAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGAt                                                                                                <  4:179068/145‑1 (MQ=255)                                                    aaTATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAAtt                                                                                                                             >  8:24672/1‑115 (MQ=255)                                                    aaTATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAAtt                                                                                                                             <  7:24672/115‑1 (MQ=255)                                                     atatCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCaga                                                                                          <  1:231236/149‑1 (MQ=255)                                                      tatCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCaga                                                                                          <  2:137274/148‑1 (MQ=255)                                                      tatCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCa                                                                                            >  1:137274/1‑146 (MQ=255)                                                           ccgccCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTc                                                                                    <  1:198331/149‑1 (MQ=255)                                                           ccgccCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTc                                                                                    <  7:128152/149‑1 (MQ=255)                                                           ccgccCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAATATTTTATCCAGCTCATCCTGCGGAATGACGATACCAGAGTCTTc                                                                                    >  6:362680/1‑149 (MQ=255)                                                               ccATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACt                                                                                <  4:199589/149‑1 (MQ=255)                                                                  tatTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTc                                                                              >  6:283616/1‑148 (MQ=255)                                                                   atTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGaaaa                                                                                                                               <  4:214031/98‑1 (MQ=255)                                                                    tttttCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAACACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCaa                                                                            >  1:250716/1‑148 (MQ=255)                                                                    tttttCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGCCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCaa                                                                            >  4:104734/1‑148 (MQ=255)                                                                          gCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATgca                                                                      <  4:315462/148‑1 (MQ=255)                                                                          gCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATgca                                                                      <  6:338296/148‑1 (MQ=255)                                                                             aGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTgg                                                                                                                                               <  1:37590/72‑1 (MQ=255)                                                                             aGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTgg                                                                                                                                               >  2:37590/1‑72 (MQ=255)                                                                               acgacgAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCatat                                                                >  3:164533/1‑149 (MQ=255)                                                                                 gacgaGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCg                                                              >  4:62987/1‑149 (MQ=255)                                                                                    gagaAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCt                                                           <  2:385661/149‑1 (MQ=255)                                                                                    gagaAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGCAATGCCGATACCACAGTCATCCACTACAAAATGCAGCATAGCGCCt                                                           >  2:139019/1‑149 (MQ=255)                                                                                    gagaAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCt                                                           >  4:371297/1‑149 (MQ=255)                                                                                      gaAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCtt                                                          >  8:89518/1‑148 (MQ=255)                                                                                       aaaCGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCa                                                        >  1:254164/1‑149 (MQ=255)                                                                                       aaaCGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCa                                                        <  3:62987/149‑1 (MQ=255)                                                                                         aCGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATc                                                      >  2:204695/1‑149 (MQ=255)                                                                                            gCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTa                                                   <  7:89518/149‑1 (MQ=255)                                                                                                gACCAATACCGGTGCCGGTGGCAGGTTTACCTCCATGACTCTCTTTCACCTAGTAATACATGGCGAAAATATTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAgcg                                                >  4:158396/1‑148 (MQ=255)                                                                                                 aCCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAgcgca                                              >  1:213884/1‑149 (MQ=255)                                                                                                   cAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGcc                                                            <  8:87513/133‑1 (MQ=255)                                                                                                   cAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGcc                                                            >  7:87513/1‑133 (MQ=255)                                                                                                      tACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAgcgcacgcgc                                         <  1:222673/149‑1 (MQ=255)                                                                                                        ccggtgccggtgGCAGGTTTACCGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAgcgcacgcgcac                                       <  7:177420/149‑1 (MQ=255)                                                                                                         cggtgccggtgGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATgcag                                                                     >  3:16719/1‑118 (MQ=255)                                                                                                         cggtgccggtgGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATgcag                                                                     <  4:16719/118‑1 (MQ=255)                                                                                                           gtgccggttGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGt                                    <  7:377362/149‑1 (MQ=255)                                                                                                               cggtgGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAAcc                                >  8:107929/1‑149 (MQ=255)                                                                                                                    gCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGcc                           >  7:276294/1‑149 (MQ=255)                                                                                                                            aCCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTgg                    <  8:81589/148‑1 (MQ=255)                                                                                                                             ccGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGt                  >  1:372286/1‑149 (MQ=255)                                                                                                                             ccGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGt                  <  3:203201/149‑1 (MQ=255)                                                                                                                              cGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCt                      >  2:212650/1‑144 (MQ=255)                                                                                                                              cGCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCt                      <  1:212650/144‑1 (MQ=255)                                                                                                                               gCCATGACTGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTg                 >  7:294392/1‑148 (MQ=255)                                                                                                                                      cTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGcgcacgcgcacg                                      >  6:217917/1‑120 (MQ=255)                                                                                                                                      cTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGcgcacgcgcacg                                      <  5:217917/120‑1 (MQ=255)                                                                                                                                      cTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGAc         <  8:243276/149‑1 (MQ=255)                                                                                                                                      cTGTCTTTCACCTGGTAACACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGAc         <  2:250716/149‑1 (MQ=255)                                                                                                                                       tGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAAt              >  6:47271/1‑143 (MQ=255)                                                                                                                                       tGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAAt              <  5:47271/143‑1 (MQ=255)                                                                                                                                       tGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACg        >  4:314518/1‑149 (MQ=255)                                                                                                                                       tGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTACACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACg        >  1:144528/1‑149 (MQ=255)                                                                                                                                       tGTCTTTCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACg        >  8:395693/1‑149 (MQ=255)                                                                                                                                           tttCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACGGCGt    >  2:114179/1‑149 (MQ=255)                                                                                                                                           tttCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACGGCGt    >  6:221538/1‑149 (MQ=255)                                                                                                                                             tCACCTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACGGCGt    <  7:395693/147‑1 (MQ=255)                                                                                                                                              ctccTAGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACGGCGTTa  <  8:18477/146‑1 (MQ=255)                                                                                                                                              cACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACGGCGTTa  <  1:139019/148‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                               GTGGATCGTCAACGTAAAGGTTGAACCTTTGCCCTGTTCGCTGGTAACCGTAATATCGCCGCCCATATTTTTCGCCAGACGACGAGAAACGGCCAGACCAATACCGGTGCCGGTGGCAGGTTTACCGCCATGACTGTCTTTCACCTGGTAATACATGGCGAAAATTTTATCCAGCTCATCCTGCGGAATGCCGATACCAGAGTCTTCCACTTCAAAATGCAGCATATCGCCTTCATCGTAGCGCACGCGCACGGTAACCTGGCCTTGCTGGGTGAATTTGACGGCGTTA  >  NC_000913/3351514‑3351802

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