Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,442,582	0	T	G	44.8%	100.4 / 13.2	68	G273G (GGT→GGG)	mnmC	fused 5‑methylaminomethyl‑2‑thiouridine‑forming enzyme methyltransferase and FAD‑dependent demodification enzyme
Reads supporting (aligned to +/- strand):  ref base T (6/31);  new base G (30/0);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.97e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTAAAC  >  NC_000913/2442433‑2442732                                                                                                                                                      |                                                                                                                                                          ggATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTa                                                                                                                                                          <  1:420396/151‑1 (MQ=255)  gATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTAt                                                                                                                                                         <  1:147679/151‑1 (MQ=255)       aCGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCa                                                                                                                                                    <  2:724322/151‑1 (MQ=255)       aCGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGAtt                                                                                                                                                        >  2:299158/1‑147 (MQ=255)        cGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAg                                                                                                                                                   <  2:492113/151‑1 (MQ=255)           tGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCcgcg                                                                                                                                                 >  2:176795/1‑150 (MQ=255)            gCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGGGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGcccacggg                                                                                                                                               >  1:755427/1‑145 (MQ=255)              aaaaaCGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCgt                                                                                                                                              <  1:436660/150‑1 (MQ=255)                 aaCGTAAGGGATTTGGGCGCAAACGGGAAATGCTTTGAGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGGTTGGCAGGGGgttgtt                                                                                                                                          >  1:423808/1‑151 (MQ=255)                       aGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCCGCGCGTTGTTGTCg                                                                                                                                      >  1:903582/1‑149 (MQ=255)                         ggCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTgg                                                                                                                                  <  1:1012654/151‑1 (MQ=255)                         ggCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTgg                                                                                                                                  <  2:45945/151‑1 (MQ=255)                           cTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGGTTGGCAGCGCGGTGTTGTTGCTcgcg                                                                                                                                >  1:334729/1‑151 (MQ=255)                           cTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGGTTTCCAGCGCGGGGTTGTTGCTgg                                                                                                                                  >  1:863147/1‑149 (MQ=255)                                 gcgcAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGTCAGCGCGGTGTTGTCGCTGGCGCTAtt                                                                                                                           >  2:950834/1‑150 (MQ=255)                                         gggAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGCGCTATTACGGCGCgg                                                                                                                  >  2:756873/1‑151 (MQ=255)                                          ggAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGCGCTATTACGGCGCgg                                                                                                                  >  2:596675/1‑150 (MQ=255)                                          ggAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGGTGGCGGCGGCGCTATTACgg                                                                                                                       >  1:697455/1‑145 (MQ=255)                                          ggAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGGTTTCCAGCGCGGTGGTGGCGCTTGGGGCATTTCGGcgcg                                                                                                                   >  1:39643/1‑149 (MQ=255)                                          ggAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGACCGCGCGTTGGTGTTTCTGGGACTAt                                                                                                                            >  1:962169/1‑140 (MQ=255)                                              aTGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTTCCAGCGCGGTGTTGTCGCTGGCGCTATTACGGCGCGGCCgg                                                                                                              >  2:422527/1‑150 (MQ=255)                                              aTGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGc                                                                                                             <  2:1408640/151‑1 (MQ=255)                                                gCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCa                                                                                                            <  1:411607/150‑1 (MQ=255)                                                gCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCa                                                                                                            <  1:422521/150‑1 (MQ=255)                                                gCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCa                                                                                                            <  1:244613/150‑1 (MQ=255)                                                    tGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTTCCCGCCCGGTTTTGTCGCTGGCGCTCTTTCCGCGCGGCCGGGAggg                                                                                                         >  1:454691/1‑148 (MQ=255)                                                       ggggTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCCCCGTGGTGTTGCTGGCGCTATTACGGCGCGGCTGGCAGGGaa                                                                                                       >  1:1133012/1‑148 (MQ=255)                                                        gggTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCCGCGCGGTGTTTTCGCTGGCGCTATTACGGCGCgggcggggggtaa                                                                                                       >  1:702422/1‑136 (MQ=255)                                                         ggTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCtt                                                                                                  <  2:1266951/151‑1 (MQ=255)                                                            gATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTAt                                                                                               <  1:449025/151‑1 (MQ=255)                                                                 aaCAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTCCAGCGCGGTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGGAGGGAAAGCtttttt                                                                                              >  2:105880/1‑147 (MQ=255)                                                                     gACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGAt                                                                                      <  1:1238605/151‑1 (MQ=255)                                                                        aTTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGTTTGCCAGCGCGGTGGTGTCGCTGGCGCTATTACGGCGCCGGTGGGCGGGAAAGCTTTTTTTTGCGGATg                                                                                     >  1:1338816/1‑149 (MQ=255)                                                                           aCCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGcc                                                                               <  2:437152/151‑1 (MQ=255)                                                                           aCCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAACGCGGTGTTGTCGCCGGCGCTATTACCGCGCGGGCGGGAGGGAACGCCTTATTTCGGGGGTGA‑ggg                                                                                >  2:737490/1‑149 (MQ=255)                                                                             cGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTTGCCGCGCGTTGTTGTCGCTGGCGCGATTTCGGCGCGGCTGGGCGGGAAAGCTTTTTTTCGCGGGTGG‑GGgccc                                                                             >  1:393728/1‑151 (MQ=255)                                                                              gCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGGGCTATTACCGCGCGGGTGGGCGGGGAAGCTTTTTTTCGCGGGTGA‑GGgcccc                                                                            >  1:389779/1‑151 (MQ=255)                                                                               cTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGGTGGTGTTGCTTGCGCTCTTTAGGGGCGGCTGGGAGGGAAACCCTTTTTTCGCGGGGGA‑GGgcccc                                                                            >  2:277003/1‑150 (MQ=255)                                                                                tCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGTTTGCCCGCGCGTTGTTGTCGCCCGGGCTCTTACGGCGCGGGTGGGGGGGTAAGGTTTTTTGCGCGGGGGA‑GGccccccc                                                                          >  1:803308/1‑149 (MQ=255)                                                                                 ccccTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCg                                                                           <  1:104150/149‑1 (MQ=255)                                                                                  cccTGCTCCGCGCCGTGGTTTAACCGTACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTGCCAGCGCGGTGTTGTCGCTGGCGCTATTTCGGCGCGGGCGGGAGGGAACGCTTTTTTTCGCGGGTGA‑GGGCCCCCCAc                                                                        >  2:218585/1‑151 (MQ=255)                                                                                        tacgcgCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGCCGCATTCCCAACCCGCTGCTGTCGCCCGCGCCATTTCGGCGCCGCCCGCAGGCAAACCCTTATTTCGAGGGTCA‑CGCCCCCCCCCTggggg                                                                  >  1:436449/3‑149 (MQ=255)                                                                                          cgcgCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑t                                                              <  2:145209/151‑1 (MQ=255)                                                                                            cgcCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGGGCTATTACGGCGCGGCTGGGAGGGAACGCTTTTTTTCGCGGGTGA‑GGGCCCCCCACTggggg                                                                  >  2:510504/1‑145 (MQ=255)                                                                                               cGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCgg                                                         <  1:1101565/151‑1 (MQ=255)                                                                                                    tttAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCCGCGCGTTGGTGTCGCTGGCGCTATTACGGGGGGGGGGGGGGGGAAAGCTTTTTTTTGCGGGTGG‑GGgcccccccccggg                                                                    >  1:321651/1‑129 (MQ=255)                                                                                                     ttAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATc                                                    <  1:938909/150‑1 (MQ=255)                                                                                                            cACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGCGCGGTGTTGTCGCTGGCGCTATTACGGGGCGGGTGGGAGGGAAAGCTTTTTTTCGCGGGGGA‑GGCCCCCCCCCTGGGG‑GCTTTCCGGC‑ATCGCCCgggg                                            >  2:432226/1‑151 (MQ=255)                                                                                                             aCGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGGCAACCCGGTGGTGGTGCTGGGGCTATTAAGGGGCGGGTGGCAGGGAAACCCTTTTTTTGCGGGTGG‑GGGCCCCCCCCCGGGG‑GG‑TTCCgggcaactcccggggg                                           >  1:1272847/1‑137 (MQ=255)                                                                                                              cGGGTAGCAGCAAACGGGAAGCGGCGATTATCGGCGGGGGGATTTCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGGAAAGCTTTTTTTCGCGGGTGA‑GGGCCCCCCACTGGGG‑GG‑TTCCCGCCATCGCCCgggggg                                          >  2:48213/1‑150 (MQ=255)                                                                                                              cGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGgc                                          <  2:755553/151‑1 (MQ=255)                                                                                                              cGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGgc                                          <  1:714268/151‑1 (MQ=255)                                                                                                              cGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGgc                                          <  1:515062/151‑1 (MQ=255)                                                                                                              cGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGGATTGCCAGAGCGGTGTTGTGGCGGGGGCTATTACAGCGCGGGTGGGCGGGAAAGCTTTTTTTCGCGGGTGA‑GGGCCCCCCCCCGGGG‑GG‑CTCCGGCACTCCCCCgggggg                                          >  1:1007879/1‑150 (MQ=255)                                                                                                                 gcagcagcaAACGGGAAGCGGCGATTATCGGCGGTGGGAATGCCAGCGCGTTGTTGTGGCTGGGGCTATTACGGGGCGGCTGGGAGGGAAATCCTTTTTTTGCGGGTGA‑GGCCCCCGCCCCGGGG‑GC‑TTCCGGGGATCGCCCggggg                                           >  2:739970/1‑147 (MQ=255)                                                                                                                       gcaAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATcc                                 <  2:1239166/151‑1 (MQ=255)                                                                                                                         aaaCGGGAAGCGGCGATTATCGGCGGTGGTATTTCCAACGCGCTGGTGTCCCTGGGGCTTTTACGGCGCGGCCCGCAGGGAAAGGGTTATTTCGCGGGCGC‑GGCCCCCGCACTGGGTGGC‑TTCCGGAAAACCCCCggggggggggatc                                   >  2:886068/1‑140 (MQ=255)                                                                                                                            cGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGttat                            <  2:687138/151‑1 (MQ=255)                                                                                                                             gggAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGttatt                           <  1:858151/151‑1 (MQ=255)                                                                                                                               gAAGCGGCGATTATCGGCGGTGGTATTTCCAGCGCGGTGTTGTCGCTGGCGCTATTACGGCGCGGGCGGCCGGGAAAGCCTTTTTTTGCGGGTGAGGGCCCCCCCCCGGGGG‑GT‑TTCCGG‑AAACTCCCGGGGGGGCTGTATTCGttttt                           >  2:1013825/1‑149 (MQ=255)                                                                                                                               gAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTaa                         <  1:1150586/151‑1 (MQ=255)                                                                                                                                    ggcgATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCaaa                    <  1:345312/151‑1 (MQ=255)                                                                                                                                       gATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAacac                 <  2:1335080/151‑1 (MQ=255)                                                                                                                                        aTTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACg                <  1:576223/151‑1 (MQ=255)                                                                                                                                        aTTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACg                <  2:926589/151‑1 (MQ=255)                                                                                                                                              ggcggTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAgg          <  1:773843/151‑1 (MQ=255)                                                                                                                                               gcggTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAgg          >  2:726800/1‑150 (MQ=255)                                                                                                                                               gcggTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGgc         >  2:1234841/1‑151 (MQ=255)                                                                                                                                                  gtggtATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCt      >  1:1028012/1‑151 (MQ=255)                                                                                                                                                      tATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTAAAc  >  1:1280595/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                          GGATTCACGATGCAAAAACGTAAGGGCTTTGGGCGCAAACGGGAAATGCTTTGCGGGGTGATGGAACAGACATTACCGCTCCCCTGCTCCGCGCCGTGGTTTAACCGCACGGGCAGCAGCAAACGGGAAGCGGCGATTATCGGCGGTGGTATTGCCAGCGCGTTGTTGTCGCTGGCGCTATTACGGCGCGGCTGGCAGGTAACGCTTTATTGCGCGGATGA‑GGCCCCCGCACTGGGT‑GC‑TTCCGGCAATCGCCAGGGGGCGCTGTATCCGTTATTAAGCAAACACGATGAGGCGCTAAAC  >  NC_000913/2442433‑2442732

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