Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,718,298	0	A	C	25.6%	94.5 / 15.2	83	G77G (GGT→GGG)	yfiF	putative methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (44/17);  new base C (6/15);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 6.53e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTT  >  NC_000913/2718152‑2718444                                                                                                                                                   |                                                                                                                                                   gCTTCCTGACGGCTCTGGAAAAGTGCATGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCAcc                                                                                                                                                  <  1:650886/149‑1 (MQ=255)  cTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCAccc                                                                                                                                                 <  1:948876/149‑1 (MQ=255)       ggACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGgcgc                                                                                                                                            >  1:233646/1‑149 (MQ=255)       ggACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTATCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATATAACTTTTACCACTGATGCCACCGTGATCGGCCTTTGCCGGCCTCTCATCACCCGgcgc                                                                                                                                            >  1:1000951/1‑149 (MQ=255)           ggCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCgg                                                                                                                                        >  2:1821574/1‑149 (MQ=255)             ctctGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGCCGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCgggc                                                                                                                                      >  1:1777149/1‑148 (MQ=255)             ctctGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTGCCGCAAGCTGAAGACGAAACACTTCCGGATCAATAAAACTTTTAGCACTGATGCCACCGTGATCGGCATTTTCCGGCGTCTCATCACCCGGAGCGCtggt                                                                                                                                      >  1:1779495/1‑148 (MQ=255)               ctGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCAGCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGaaa                                                                                                                                    >  1:109479/1‑149 (MQ=255)                tGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAAc                                                                                                                                   >  1:354669/1‑149 (MQ=255)                tGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGAGGGTTTCTTCCGCACGATTACGACGCAACACTTCCGGATCAATAAATCTTTTACAACTTATGGCACCGTGATCGGAATTTTACTGCGTCTCATCACACGGCTCGCGGGAAAc                                                                                                                                   >  2:1209617/1‑149 (MQ=255)                  gAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCg                                                                                                                                 >  1:559225/1‑149 (MQ=255)                    aaaGTGCCTGACATTCATTTTCGCTGTATAAGATGGTTTAAACCGCACGCTGAAGAATCAACACTTACGGATCAATAAAACTTTTACCACTGATGCCAACGTGATCGGCCCTTTCCGGCGTCTCATCACCCGGCACGCGGGAAACCGTg                                                                                                                               <  2:1219005/149‑1 (MQ=255)                     aaGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCATCACTTCCGGATCAATAAAACTTTTAGGACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTgc                                                                                                                              <  1:161745/149‑1 (MQ=255)                     aaGTCCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGAAACACTTCCGGATCAATAAAACTTTTACCACAGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTgc                                                                                                                              <  1:1000943/149‑1 (MQ=255)                          ccTGACATGCATTTTCGCCGTAGCCGCGGGTTTCTTCCGCACGCTGCCGACGCAACACTTCCGGCTCAATAAAACTTTGACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAGCCGTGCGCCAc                                                                                                                         >  1:1104745/1‑149 (MQ=255)                           cTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGCCGACGCAACACTTCCGGATCAATAAAACTTTTACCACTCATGCCACCGTGATCGGCCTTTTCCCTCGTCTCATCCCCCGGCGCGCGGGCAGCCGTGCGCCACg                                                                                                                        >  2:1473866/1‑149 (MQ=255)                            tGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCAGAAGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCGCgg                                                                                                                       >  1:682813/1‑149 (MQ=255)                                  gCATTTTCTCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCGACGGCGATTc                                                                                                                 >  1:966251/1‑149 (MQ=255)                                  gCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAAGACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATc                                                                                                                 >  1:733081/1‑149 (MQ=255)                                   cATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCACCGGCGAATCt                                                                                                                >  1:442303/1‑149 (MQ=255)                                    aTTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACTACGCAACACTTCGGGATCAATAAAACTTTTACCACTGATGCCCCCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCtt                                                                                                               <  1:59815/149‑1 (MQ=255)                                    aTTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATACTACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGACACCGTGCGCCACGGCGAATCtg                                                                                                               >  1:394417/1‑148 (MQ=255)                                    aTTTACGCCGTAGACGCGGGTTTCTTCCCCACGCTGACGACGCAACACTACCGGATCAATAAAACTTTTCCCACTGATGCCCCCGTGATCGCCCTTTTCCGGCGTCACATCACCCGGCGCGCGGGAAACCGGGCGCCACGGCGAATCtt                                                                                                               <  1:785362/149‑1 (MQ=255)                                      tttCCCCGTAGACGGGGGTTTCTTCCGCACGCTGCCGACGCAACCCTTCCGGAACAATAAAACTTTACCCACTCATGCCACCGTGAGCGGCCTTTGCCGGCGTCTCACCACCCGGCGCGCGGTAAACCGTGCGCCACGGCGAATCTTcc                                                                                                             <  1:624170/149‑1 (MQ=255)                                          gCCGTAGCCGCGGGTTTCCTCCGCACGCTGACGACGCAACACGTCCGGATCAATAACCCTGTTCCCACTGATGCCACCGTGCTCGGCCGGTTCCGGCGTCTCATCACCCGGCGCGCGGGAAGCCGTGCGCCCCGGCGCATCTTCGCAcg                                                                                                         >  1:1625375/1‑148 (MQ=255)                                          gCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGCATCTTCCCACt                                                                                                         >  2:1387564/1‑149 (MQ=255)                                          gCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACt                                                                                                         >  2:1599796/1‑149 (MQ=255)                                            cGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCTCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCagtgg                                                                                                       >  2:1106438/1‑145 (MQ=255)                                             gtcgccgcgGGTTGCTCCCGCACGCGGACGACGCAACACTTCCGGAGCAATAAAACTTTTACAACTGAGGCCCCCGTGAGCGGCCTTTTCGGGCGGCGCCTCCCCGGGCGCGCGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGc                                                                                                      <  1:105542/144‑1 (MQ=255)                                             gTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATATAACTTTTACCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCCCTTGc                                                                                                      >  2:223894/1‑149 (MQ=255)                                              tAGAGGCGGGTTTCTTCCGCACGCTTACCACGCAACACTTCCGTATCAATAAAACTTTTACCACTGACGCCCCCGTTATCGGCCTTTTCCGGCTTCTCATCCCCCGGCGCGCGGTAAACCGTGCGCCACGGCGAATCTTCCCACTTGCg                                                                                                     <  2:1588602/149‑1 (MQ=255)                                                  cgcgGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCa                                                                                                 >  1:1494961/1‑149 (MQ=255)                                                   gcgGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCCCGGCGAATCTTCCCACTTGCGGTCAc                                                                                                >  1:1849567/1‑149 (MQ=255)                                                    ggggTTTATTCCGCACGCTGACGACGCAACCCTTCCGGAGCAATAAAACTTTTCCCACTGATGCCCCCGTTCTCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACg                                                                                               <  1:540794/148‑1 (MQ=255)                                                     gggTTTCTTCCGCACGCTGACGCCGCAACCCTTCCGGATCAATAAAACTTATACCACTGATCCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACgg                                                                                              <  2:301146/149‑1 (MQ=255)                                                            ttcCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGAc                                                                                       >  2:301318/1‑149 (MQ=255)                                                                gCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTg                                                                                   >  1:873629/1‑149 (MQ=255)                                                                gCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTg                                                                                   >  1:71616/1‑149 (MQ=255)                                                                gCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTACCACTTACGGTCACGGGGCTGACTCTg                                                                                   >  2:1052270/1‑149 (MQ=255)                                                                 cACGCTGAGTCCGAAACACTTCCGTACCAATAAATCTTTTACCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGGCTCATCACCCGGCGGGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGt                                                                                  <  2:1825183/149‑1 (MQ=255)                                                                 cACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCCCCGTGATCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGGCACGGGGCTGCCGCTGt                                                                                  >  2:930690/1‑149 (MQ=255)                                                                    gCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTg                                                                               >  1:565765/1‑149 (MQ=255)                                                                      tGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTTATCGGCCTTTTCCGGCTTCTCATCACCCGGCGCGAGGGAAACCGAGCGCCACGGCGAATCTTCCCACTTGATGTCACGGGGCTGACTCTGTTTGta                                                                             >  2:463187/1‑148 (MQ=255)                                                                        acgacgCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGATTGTCAt                                                                           >  1:58582/1‑149 (MQ=255)                                                                          gacgCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCCCCGTGATCGGCCTGTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGCCTCGGGTTGTCATcg                                                                         >  2:1021837/1‑149 (MQ=255)                                                                          gacgCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATcg                                                                         >  1:559623/1‑149 (MQ=255)                                                                               aaCACTTCCGGATCAATAAACCTTTTACAACTGATGCCACCGTCATCGGCCTATTCCGGCGTCTCATCACACGGCGCGCCGGAAACCGTTCGGCACGGCGAATCTACCCACTTGCCGTCACGGGGCTGACACTGTTTGTCATCgcggcc                                                                    >  1:1387128/1‑147 (MQ=255)                                                                                 acctcacGGATCAATATATGTATTCCCAGTGATGCCACCGTTCGGGGCGTTTTCGGGCGTCGAAGACCCCGGCGCGGGGGAACGCGTGGGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCgcgggcgg                                                                  <  1:1728373/143‑1 (MQ=255)                                                                                  acTTCCGCATCAAAAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGcgggcggg                                                                 <  2:1125751/149‑1 (MQ=255)                                                                                    ttCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCgggcggggc                                                               >  2:474018/1‑149 (MQ=255)                                                                                    ttCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCgggcggggc                                                               >  1:236638/1‑149 (MQ=255)                                                                                          ctCAATAAACCTATTCCCCCAGCGGCCCCCGTGAGCGGCCTGTTCCGGGGTCTCATCCCCCGGCGCGCGGTAACCCGTGCCCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGc                                                         <  2:1804272/148‑1 (MQ=255)                                                                                          aTCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGc                                                         >  1:688660/1‑149 (MQ=255)                                                                                           tcaatcaaCCTTTTCCCCCAGATGCCACCGTCCTCGGCCTTTTCCGGAGTGTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGcc                                                        <  2:1198826/149‑1 (MQ=255)                                                                                              aTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAGCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGTGCGACGGCCAcc                                                     >  2:111822/1‑149 (MQ=255)                                                                                              aTAAAACTATTACCACTGCTGCCACCGTGATCGGCCTTTACCGGCGTCTCATCACCCGGCGCGCGGGATACCGTGCGCGACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGCCGGGGCGACGGCCAcc                                                     <  1:1052270/149‑1 (MQ=255)                                                                                                aaaaCTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGCATCTTCCCCCTTGCGGTTACGCGGCTGACTCTGTGTGTCATCGCGGGCGGGGCGACGGCCACCGt                                                   >  1:1080041/1‑149 (MQ=255)                                                                                                 aaaCTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTc                                                  >  1:1322294/1‑149 (MQ=255)                                                                                                  aaCTTTTCCCCCTGAGCCCCCCGTGATCGCCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCa                                                 <  1:817110/149‑1 (MQ=255)                                                                                                     ttttACCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTCGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCa                                              <  1:242664/149‑1 (MQ=255)                                                                                                     ttttACCACTGATGACTCCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGCACACCGTGCGCCACGGCGATGCTTCCCACGTGCGGGCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGcc                                              >  2:1018244/1‑148 (MQ=255)                                                                                                     ttatACCCCTGATGCCACCTTGACGGGCCTTTACCGGGGACTCATCACCCGGTGCGCGGGAACCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCa                                              <  2:1000951/146‑1 (MQ=255)                                                                                                      tttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCCCGGGGCTGCCTCTGGTTGTCATCGCGGGCGGGGCGACGGCCCCCGTCAGCAc                                             >  1:1533096/1‑149 (MQ=255)                                                                                                        tACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCCCCGTCAGCACga                                           >  2:1808911/1‑149 (MQ=255)                                                                                                           cACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCCCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGAt                                        >  2:828361/1‑149 (MQ=255)                                                                                                               gATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTc                                    <  1:1106438/149‑1 (MQ=255)                                                                                                               ccgccccccGTGAGCGCCCTATCCCGCCGTATAATCCCCCGGCGCGCGGCAAACGGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTc                                    <  1:1387564/142‑1 (MQ=255)                                                                                                               catcccACCGTTCCCGGCCTTTCCCGGCGTCCCACCCCCGGGCGCGCGGCAAACGGGCCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTc                                    <  1:1599796/145‑1 (MQ=255)                                                                                                                ctccccccGTGCGCGGCCTTGACCGCCGGCGCCTCCCCCGCCGGGGGGTAAACCGTGCCCCACGGCGAATCTTCCCACTTGCGGCCCCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTcc                                   <  2:1625375/143‑1 (MQ=255)                                                                                                                 tGCCACCGTGCGCGCCCTTTTCCGGCGTACCCTACCCCGCCGCGCGGCAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCt                                  <  1:341005/149‑1 (MQ=255)                                                                                                                  gCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCCCTTGCGGGCACGGGGCTGACTCTGTTTGTCATCGCGGGCAGGGCGCCGGCCACCGTCAGCCCGAGATTTTCCTg                                 >  1:761071/1‑149 (MQ=255)                                                                                                                  ccccccGTGACCGCCGTTTTCCGCCTTCTCAACCCCGGGCGCGCGGCACCCCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTg                                 <  2:1578673/145‑1 (MQ=255)                                                                                                                   ccACCGTGATCGGCCTTTTCCGGCGTCTCATCCCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTgg                                >  2:1494534/1‑149 (MQ=255)                                                                                                                    cACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGa                               >  2:251644/1‑149 (MQ=255)                                                                                                                            tCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCCCCGTCAGCACGAGAATTTACTGGAcgcccgcc                       >  2:1589653/1‑149 (MQ=255)                                                                                                                              ggCCTTACCCGGCGCCGCACCCCCCGGCGCGCGGGAACCGGTGCGCCACGGCGAATCTTCCCCCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCGTGGACGCCCGcccc                     <  2:888721/149‑1 (MQ=255)                                                                                                                               gcgttttCCGGCGTCTCGTCCCCGGGCGCGCGGGAACCCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCCCGGGGCGGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGccccc                    <  2:589417/146‑1 (MQ=255)                                                                                                                                    taccGGCTTCCCACCCCCCGGCGCGCGGGAAACCGGGCGCCACGGCGAATCATCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTcc               <  1:1021837/147‑1 (MQ=255)                                                                                                                                       cGGCGTCTCATCACCCGGCGCGCGGGTAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGGTGACTCTGTTTGTCCTCGCGGGCGGGGCGGAGGGCGCCGTCTGCCGGTGATTGTCGTGGACGCCCGCCCCCTGTGCCgg            >  1:1031917/1‑149 (MQ=255)                                                                                                                                       cGGCGCCGCCTCACCCGGCGCGCGGTAACCGGGGCGCCACGGCGAATCTTCCCACTTGCGGTCGCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCgg            <  1:1821574/149‑1 (MQ=255)                                                                                                                                         gCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTa          >  2:526915/1‑149 (MQ=255)                                                                                                                                          cGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTAc         >  2:846037/1‑149 (MQ=255)                                                                                                                                            tctcATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTAcgc       <  2:377878/149‑1 (MQ=255)                                                                                                                                                 tcaCCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCCCCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGtt  >  1:988450/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                   GCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGACGCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCGGCGCGCGGGAAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCTGGACGCCCGCCCCCTTTCCCGGTACGCGGGTT  >  NC_000913/2718152‑2718444

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