Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,912,153	T→C	5.3%	S142G (AGC→GGC)	htpX ←	putative endopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,912,153	0	T	C	5.3%	291.4 / 4.7	94	S142G (AGC→GGC)	htpX	putative endopeptidase
Reads supporting (aligned to +/- strand):  ref base T (44/44);  new base C (2/3);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

TGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCA  >  NC_000913/1912006‑1912299                                                                                                                                                    |                                                                                                                                                   tGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGACCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCt                                                                                                                                                    >  7:268559/1‑148 (MQ=255)    cGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGAt                                                                                                                                                 >  4:353590/1‑148 (MQ=255)      ttGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCTCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGAttt                                                                                                                                               >  3:256982/1‑148 (MQ=255)        gCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTTGCTGATTTCg                                                                                                                                             >  2:345673/1‑148 (MQ=255)         ctctctTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCg                                                                                                                                             <  3:204382/147‑1 (MQ=255)         ctctctTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGt                                                                                                                                            <  8:281503/148‑1 (MQ=255)            tctTCACCTTCATCACGATTTCCGCCCATAAAACCCTCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTg                                                                                                                                           <  4:256982/146‑1 (MQ=255)             ctTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAg                                                                                                                                         >  7:258726/1‑147 (MQ=255)             ctTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAg                                                                                                                                         >  4:265254/1‑147 (MQ=255)             ctTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCCGATTTCGTGAg                                                                                                                                         >  7:236020/1‑147 (MQ=255)             ctTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCCGATTTCGTGAg                                                                                                                                         >  7:241462/1‑147 (MQ=255)                cACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGc                                                                                                                                        <  8:288338/145‑1 (MQ=255)                   cTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTAc                                                                                                                                  >  5:44492/1‑148 (MQ=255)                   cTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTAc                                                                                                                                  >  7:334041/1‑148 (MQ=255)                    ttCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTAcc                                                                                                                                 >  1:327870/1‑148 (MQ=255)                    ttCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTAcc                                                                                                                                 >  4:21396/1‑148 (MQ=255)                         caCGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTc                                                                                                                            >  8:312691/1‑148 (MQ=255)                         caCGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTc                                                                                                                            >  4:68342/1‑148 (MQ=255)                           cGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCgg                                                                                                                          >  1:142615/1‑148 (MQ=255)                            gATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGc                                                                                                                         <  2:15043/148‑1 (MQ=255)                                  cGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATc                                                                                                                   <  6:71896/148‑1 (MQ=255)                                  cGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATc                                                                                                                   <  3:68342/148‑1 (MQ=255)                                  cGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCCGATTTCGTGAGCAATTACCGCCTCGGCTTCATc                                                                                                                   <  8:236020/148‑1 (MQ=255)                                  cGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCCGATTTCGTGAGCAATTACCGCCTCGGCTTCATc                                                                                                                   <  8:241462/148‑1 (MQ=255)                                   gCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATcc                                                                                                                  <  2:91373/148‑1 (MQ=255)                                        tAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCt                                                                                                             >  8:47523/1‑148 (MQ=255)                                         aaaaCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCg                                                                                                                           <  7:357407/133‑1 (MQ=255)                                         aaaaCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCg                                                                                                                           >  8:357407/1‑133 (MQ=255)                                         aaaaCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCt                                                                                                             >  8:101358/1‑147 (MQ=255)                                         aaaaCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCt                                                                                                             >  8:343634/1‑147 (MQ=255)                                         aaaaCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGCTGACCAAATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCt                                                                                                             >  4:80521/1‑147 (MQ=255)                                            aCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTg                                                                                                                                           >  2:8205/1‑114 (MQ=255)                                            aCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTg                                                                                                                                           <  1:8205/114‑1 (MQ=255)                                            aCCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATg                                                                                                         >  8:25103/1‑148 (MQ=255)                                               cgcgGCAAGCTGCGCGAGAATACGGGAAATAAAGGTAACGAAGGTGTTCACCCCGCCCTGAATGCGCGTCATGGTGCCCATATCAGCATTGGCGATGTGGGGGCTTTCGTGAGCAATTACCGCCTCGGGTTCTTCCGGGCTCATGTTc                                                                                                      >  8:172020/1‑148 (MQ=255)                                               cgcgGCAAGCTGCGCCAGAATACGGGAAATAAAGATATCGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAAt                                                                                                                                     <  2:58331/117‑1 (MQ=255)                                               cgcgGCAAGCTGCGCCAGAATACGGGAAATAAAGATATCGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAAt                                                                                                                                     >  1:58331/1‑117 (MQ=255)                                               cgcgGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTc                                                                                                      >  8:99288/1‑148 (MQ=255)                                                gcgGCAAGCTGCGCCAGAATACGGGAAATAAATATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCCGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCt                                                                                                     <  4:207679/148‑1 (MQ=255)                                                gcgGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCt                                                                                                     <  7:78013/148‑1 (MQ=255)                                                     aaGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTgcagc                                                                                                <  8:12247/148‑1 (MQ=255)                                                        cTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCaaa                                                                                             <  4:41202/148‑1 (MQ=255)                                                          gcgcCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAAcc                                                                                           >  7:92787/1‑148 (MQ=255)                                                             ccAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGt                                                                                        <  8:334041/148‑1 (MQ=255)                                                             ccAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGt                                                                                        >  6:276805/1‑148 (MQ=255)                                                              cAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTg                                                                                       <  3:353590/148‑1 (MQ=255)                                                                 aaTACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCt                                                                                     >  5:195157/1‑147 (MQ=255)                                                                 aaTACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTg                                                                                    <  8:258726/148‑1 (MQ=255)                                                                 aaTACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTg                                                                                    <  8:125742/148‑1 (MQ=255)                                                                  aTACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGa                                                                                   >  3:39986/1‑148 (MQ=255)                                                                      gggAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAg                                                                                >  3:145147/1‑147 (MQ=255)                                                                       gcaaaTAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACAATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGGTCTGCAGCAAACCGGTGCTGAAAGc                                                                               >  1:14794/3‑147 (MQ=255)                                                                         aaaTAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAAc                                                                            >  2:330339/1‑148 (MQ=255)                                                                              aaGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCagag                                                                       <  7:183490/148‑1 (MQ=255)                                                                              aaGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCaga                                                                        <  3:265254/147‑1 (MQ=255)                                                                                 aTAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAgg                                                                    <  8:216062/148‑1 (MQ=255)                                                                                          gtgtTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGgc                                                           <  2:142615/148‑1 (MQ=255)                                                                                             ttCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCACCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGgcgcg                                                        >  8:217044/1‑148 (MQ=255)                                                                                              tCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGgcgcgc                                                       >  4:209802/1‑148 (MQ=255)                                                                                               caccacGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCaa                                                                             <  1:125528/125‑1 (MQ=255)                                                                                               caccacGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCaa                                                                             >  2:125528/1‑125 (MQ=255)                                                                                                 ccacGCCCTGAATCAGCGTCATGGTGACCATATCAGCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCAcc                                                    <  8:268559/148‑1 (MQ=255)                                                                                                  cacGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCg                                                   >  3:92126/1‑148 (MQ=255)                                                                                                   acGCCCTGAATCAGCGTCATGGTGACCATATCCCCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACCGCAAACAGAGAGGCATCACGGCGCGCCCCgg                                                  >  4:39213/1‑148 (MQ=255)                                                                                                   acGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTCCCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCTGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCgg                                                  <  3:39213/148‑1 (MQ=255)                                                                                                    cGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTCCCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGt                                                 <  5:110185/148‑1 (MQ=255)                                                                                                      cccTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTg                                               <  7:343634/148‑1 (MQ=255)                                                                                                        cTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGc                                              >  8:287544/1‑147 (MQ=255)                                                                                                         tGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCa                                             >  2:315761/1‑147 (MQ=255)                                                                                                         tGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCa                                             >  6:23212/1‑147 (MQ=255)                                                                                                           aaTCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCaaaa                                          <  4:255976/148‑1 (MQ=255)                                                                                                           aaTCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCaaa                                           <  6:297737/147‑1 (MQ=255)                                                                                                            aTCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAg                                         >  4:180130/1‑148 (MQ=255)                                                                                                               aGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGt                                      <  8:92787/148‑1 (MQ=255)                                                                                                               aGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGt                                      <  7:25103/148‑1 (MQ=255)                                                                                                               aGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAAATACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGATAGGCATCACGGCGCGCACCGGTTGCAAAAGCg                                       <  7:172020/147‑1 (MQ=255)                                                                                                                gCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGt                                      >  3:282662/1‑147 (MQ=255)                                                                                                                gCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGt                                      <  4:282662/147‑1 (MQ=255)                                                                                                                 cGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTg                                    <  7:287544/148‑1 (MQ=255)                                                                                                                 cGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTg                                    <  1:156438/148‑1 (MQ=255)                                                                                                                  gTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCATTTACCGCCTCGGCTTCATCCGGGCTCATGGTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGa                                   >  5:321340/1‑148 (MQ=255)                                                                                                                  gTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGa                                   <  6:44492/148‑1 (MQ=255)                                                                                                                     aTGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGt                                <  3:76306/148‑1 (MQ=255)                                                                                                                            ccATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGgcgc                         <  2:144923/148‑1 (MQ=255)                                                                                                                                atCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATgg                     >  8:79309/1‑148 (MQ=255)                                                                                                                                 tCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATgg                     <  3:180130/147‑1 (MQ=255)                                                                                                                                 tCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGt                    <  1:315761/148‑1 (MQ=255)                                                                                                                                   aCCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTa                   <  3:21396/147‑1 (MQ=255)                                                                                                                                     cATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGAt                <  4:39986/148‑1 (MQ=255)                                                                                                                                       ttGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCg                             <  6:262859/133‑1 (MQ=255)                                                                                                                                       ttGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCg                             >  5:262859/1‑133 (MQ=255)                                                                                                                                        tGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGc             >  3:316076/1‑148 (MQ=255)                                                                                                                                        tGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGc             >  4:221905/1‑148 (MQ=255)                                                                                                                                                  gCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGc   <  5:23212/148‑1 (MQ=255)                                                                                                                                                   cTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCCCTTGCGGCa  >  5:208746/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TGCCGTTGCTCTCTTCACCTTCATCACGATTTCCGCCCATAAAACCCGCGGCAAGCTGCGCCAGAATACGGGAAATAAAGATAACGAAGGTGTTCACCACGCCCTGAATCAGCGTCATGGTGACCATATCACCATTGGCGATGTGGCTGATTTCGTGAGCAATTACCGCCTCGGCTTCATCCGGGCTCATGTTCTGCAGCAAACCGGTGCTGACAGCAACCAGAGAGGCATCACGGCGCGCACCGGTTGCAAAAGCGTTGATGTCCGGCGCATGGTAGATAGCCACTTGCGGCA  >  NC_000913/1912006‑1912299

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