Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,445,981	0	T	G	33.3%	125.6 / 13.5	60	G34G (GGT→GGG)	ydbL	DUF1318 family protein
Reads supporting (aligned to +/- strand):  ref base T (11/29);  new base G (19/1);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 7.96e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGAAGTCGCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATA  >  NC_000913/1445747‑1446208                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      aTTGAAGTCGCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGaaa                                                                                                                                                                                                                                   >  1:100021/1‑239 (MQ=255) aTTGAAGTCGCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGaa                                                                                                                                                                                                                                    <  1:794636/238‑1 (MQ=255)          gCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCCGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGCTGAAGCCAGAACTCAGGGGCGGGTAGGGGGAAACTTTTatg                                                                                                                                                                                                                          >  2:54770/1‑237 (MQ=255)          gCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGATAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACg                                                                                                                                                                                                                          >  1:1325215/1‑239 (MQ=255)          gCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGAGAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACg                                                                                                                                                                                                                          >  1:1316851/1‑239 (MQ=255)                    aGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGt                                                                                                                                                                                                                >  1:1116246/1‑239 (MQ=255)                                      atatGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGAt                                                                                                                                                                                                >  1:375819/1‑237 (MQ=255)                                           aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTTCGGTTATCTGGTTGCGCTGAAAACg                                                                                                                                                                                                   >  1:890204/1‑229 (MQ=255)                                           aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTTCGCTGAAAACGGACGGT‑Gaaa                                                                                                                                                                                        >  1:876588/1‑239 (MQ=255)                                           aaCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACg                                                                                                                                                                                                   <  2:890204/229‑1 (MQ=255)                                            aCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑Gaaa                                                                                                                                                                                        <  2:1316851/238‑1 (MQ=255)                                            aCGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑Gaaa                                                                                                                                                                                        <  2:1325215/238‑1 (MQ=255)                                                            tGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGcc                                                                                                                                                                       <  2:1161046/239‑1 (MQ=255)                                                                       tcaAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGc                                                                                                                                                            <  2:1116246/239‑1 (MQ=255)                                                                               gACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGTGGAACCTTTTTCGGGTTTCTGGTTGCGCCGAAAAAGGATGCTGGAAACA‑AGAAATTTGTAGCCGGCATTAATGCCGCACgga                                                                                                                                                    >  2:870282/1‑237 (MQ=255)                                                                                aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑g                                                                                                                                                                                           <  1:1155902/199‑1 (MQ=255)                                                                                aCAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTACGGTTATCTGGGTGCGCCGAAAACGGATGCG‑GAAACAGAGAAATTAGTAGCCCACATTAATGTCGAACGCaa                                                                                                                                                   >  2:11913/1‑239 (MQ=255)                                                                                  aaaGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACa                                                                                                                                                                                      <  1:895455/202‑1 (MQ=255)                                                                                          cGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTAc                                                                                                                                         <  2:375819/239‑1 (MQ=255)                                                                                          cGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGGGAAAAATTTTTCGGGTATCTGGTTTCGCTGAAAACGGGAGGT‑GTAAAAGAGaa                                                                                                                                                                                 >  2:1155902/1‑199 (MQ=255)                                                                                                gctgctTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAAAATTTTACGGTTATTTGGTTGCGCTGAAAATGGGTGCT‑GCAAAGAAGAAAATAGTAGTCTACATTAATTCCGAACGGAAAGTGAggttacaaaaa                                                                                                                                   >  2:15240/1‑229 (MQ=255)                                                                                                     ttGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTgg                                                                                                                               <  1:11913/238‑1 (MQ=255)                                                                                                      tGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGgca                                                                                                                             <  1:458474/239‑1 (MQ=255)                                                                                                              gTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGGAACATTTTACGGTTATCTGGGTGCGTTGTAAACGGATGGC‑GAAACAGAGAAATTAGTTGCCGGCATTAATTGCGGACGGCAAACGAGTTAACAACAAATGGCAAAGCaaaa                                                                                                                     >  2:935734/1‑239 (MQ=255)                                                                                                                     tCTTTTCTGATGTGATGATGAAGAAAACATTACTTCTTTGGTCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAAtgt                                                                                                               <  2:546582/238‑1 (MQ=255)                                                                                                                        tttCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCg                                                                                                           <  2:100021/239‑1 (MQ=255)                                                                                                                                          agaaAACATTACTTCTTTGCGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGTGAAAAAGTTTACGGTTATCTGGTTGCGCTGAAAACGGAAGCT‑GAAACAAACACATTAGTAGCCGACATTAAAGCCGAAAGGAAAGCGAGTTAACAACAACCGGCAAAGCAAAATAAAGTGGTGGGAGACGGAATTGCGaaa                                                                                         >  2:1112512/1‑239 (MQ=255)                                                                                                                                              aaCATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCGGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCg                                                                                     <  1:347218/239‑1 (MQ=255)                                                                                                                                                           tgtgCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAAtg                                                                                                                <  2:950020/199‑1 (MQ=255)                                                                                                                                                           tgtgCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGGTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAAtg                                                                                                                >  1:950020/1‑199 (MQ=255)                                                                                                                                                             tgCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATGTGGTTGCGCTGAAAGCGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATcg                                                                                              <  1:1034887/215‑1 (MQ=255)                                                                                                                                                             tgCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTTCGGTTATTTGGTTGCGCTGAAAACGGATGCT‑GAAAAAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATACCGAAAATCGCCGGGCCAAAACTCg                                                                      >  2:866145/1‑239 (MQ=255)                                                                                                                                                                                   agcagcAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGa                                               <  1:647303/239‑1 (MQ=255)                                                                                                                                                                                     cagcaATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGgca                                                                                                                             <  2:383461/160‑1 (MQ=255)                                                                                                                                                                                     cagcaATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTTGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGgca                                                                                                                             >  1:383461/1‑160 (MQ=255)                                                                                                                                                                                        caATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGaa                                                                                                                                                                                                                                    <  2:68537/55‑1 (MQ=255)                                                                                                                                                                                        caATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGaa                                                                                                                                                                                                                                    >  1:68537/1‑55 (MQ=255)                                                                                                                                                                                          aTGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTACGGTTATCTGGGTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAACGAGTTACCAAAAAATGGCAAAGCAAAAAAATGGGTCGGGAGATGATATCGCGAAACTCGCCGGGCAAAAACTCGTAGCTC‑GGGCCAAAACGGGGGGATATgga                                        >  1:400334/1‑237 (MQ=255)                                                                                                                                                                                           tGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGATCGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGa                                               <  1:830263/231‑1 (MQ=255)                                                                                                                                                                                           tGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTAc                                       <  2:83900/239‑1 (MQ=255)                                                                                                                                                                                              aaTGGCATTGACTCGGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACaa                                     <  1:870282/238‑1 (MQ=255)                                                                                                                                                                                               aTGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGGAGGGGAAAAATTTTTCGGTTATGTTGGTGGGCTGAAAAAGGGTGGT‑GGAAAAGAGAAATTTGTAGCCGACAATAATGCCGAACGCAAAGCGAGTTACCAAAAACTGGCAAAGCAAAAAAATTAGGCGGGAGATGATATATCGAAAATGGCCGGGCAAAAAATCGTGGCAC‑GGGGCAAAACCGGGGGAGATGTTCaaag                                   >  2:20238/1‑237 (MQ=255)                                                                                                                                                                                                     ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATaaa                             <  1:866145/239‑1 (MQ=255)                                                                                                                                                                                                     ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAAAATTTTACGGGTTTTTGGTTGCGCTGAAAACGGGTGCT‑GAAACAGAGAAAATAGTAGGCGGCATTAATGGCGAACGTAAAGCGAGTTACCAAAAACTGGCAAAGAAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAACCCGTAGATCGGGGCAAAACCGGGGGAA‑ATTTACAAAGGATaaa                             >  2:1382671/1‑239 (MQ=255)                                                                                                                                                                                                     ttGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGGGAAAAAATTTACGGTTATCTGGGTGGGCTGAAAACGGATGCT‑GAAAAAAAGAAAATAGAAGTCGACATTAATGACCAAAGGAAAGCGAGTTAACAAAAAATGGCAAAGCAAAATAATGTGTCGGGAAAGGATATCGCGAAAATCGCCCGGCAAAAACTCGTGGCCC‑GGGGCAAATCGGGGGAAAATGTACAAgaag                                 >  2:830263/1‑232 (MQ=255)                                                                                                                                                                                                       gACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACg                           <  1:935734/239‑1 (MQ=255)                                                                                                                                                                                                        aCTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGcc                                                                                                                                                                       <  1:1067525/99‑1 (MQ=255)                                                                                                                                                                                                         ctctGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATTCT‑GAAACAGAGAAATTAGTAGcc                                                                                                                                                                       <  1:816553/98‑1 (MQ=255)                                                                                                                                                                                                                 gAAGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGt                 >  1:534333/1‑239 (MQ=255)                                                                                                                                                                                                                  aaGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTACGGTTATCTGGTTGCGCTGGAAACGGGAGCT‑GAAAAAGAGAAAATAGTAGCCGACATTaaag                                                                                                                                                             >  2:1067525/1‑97 (MQ=255)                                                                                                                                                                                                                   aGCCAGAACTCAGGGGCGGGTAGGGGAAACATTTTTCGGTTTTTTGGTTGCGCTGAAAAAGGATGCC‑GAAACAGAGGAATTAGTAGACGAACTTAAta                                                                                                                                                             >  2:816553/1‑97 (MQ=255)                                                                                                                                                                                                                    gCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACaa                                     <  1:45014/216‑1 (MQ=255)                                                                                                                                                                                                                        gAACTCAGGGGTGGGTAGGGGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGCAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAat                                             >  1:1360943/1‑204 (MQ=255)                                                                                                                                                                                                                        gAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGcc                                                                                                                                                                       <  2:544652/83‑1 (MQ=255)                                                                                                                                                                                                                        gAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAat                                             <  2:1360943/204‑1 (MQ=255)                                                                                                                                                                                                                        gAACTCAGGGGCGGGTAGGGGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGcc                                                                                                                                                                       >  1:544652/1‑83 (MQ=255)                                                                                                                                                                                                                         aaCTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATa                               <  1:805664/217‑1 (MQ=255)                                                                                                                                                                                                                              aGGGGCGGGGAGGGGAAACATTTTTCGGGTATCCGGTTGCCCTGAAAACGGATGCT‑GGAACAGAGAAATTAGTTGCCCACATTAATGCCCAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATCTTTCGGTAGATGATATCCCGAAACTCGCCGGGCAAAAACTCGGAGCTC‑CGGCCCAACCGGGGGGATATGTACCAGGGATCAACg                           >  2:45014/1‑216 (MQ=255)                                                                                                                                                                                                                                gggCGGGTAGGTGAAACATTTTACGGTTATCTGCTTGCGCTGAAAAAGGGTCTT‑GAAACAGAGAAATTAGTAGCCGCAAATTATGCCGAACGTAAAGCGAGTTTCCAACAAATGGCAAAGCAAAAAAATGTGTCCGTAGATGATAATCCGAGACTCGCCGGGCCAAAGATCGTAGCTC‑GGGGCAAAACGGGGGGATATGTGCCAGTGATAAAAGGTTAATAGGGGCGGAAAATTTGATa  >  2:1046901/1‑239 (MQ=255)                                                                                                                                                                                                                                   cGGGGAGGGGAAACATTTTACGGTTATCTGGTTGCGCTGAAAAAGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACACGCCAGGGAAAAACACGTAGCTC‑GGGGCAAAACCGGGGAATATGTACAAGGGATAAACGGGAAAt                     >  2:805664/1‑217 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                      ATTGAAGTCGCTGCACCTAAGGAACCGATCACTATCAATATGAACGTTAAAATTGAGCATGAGATCATCATCAAGGCAGACAAAGATGTCGAAGAGCTGCTTGAAACTCGTAGCGATCTTTTCTGAGGTGATGATGAAGAAAACATTACTTCTTTGTGCGTTTCTTGTTGGGCTGGTAAGCAGCAATGTAATGGCATTGACTCTGGATGAAGCCAGAACTCAGGGGCGGGTAGGTGAAACATTTTACGGTTATCTGGTTGCGCTGAAAACGGATGCT‑GAAACAGAGAAATTAGTAGCCGACATTAATGCCGAACGTAAAGCGAGTTACCAACAACTGGCAAAGCAAAATAATGTGTCGGTAGATGATATCGCGAAACTCGCCGGGCAAAAGCTCGTAGCTC‑GGGCCAAACCGGGGGAATATGTACAAGGGATAAACGGTAAATGGGTGCGAAAATTTTAATA  >  NC_000913/1445747‑1446208

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