Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,663,466	G→T	6.6%	V9L (GTG→TTG)	suhB →	inositol monophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,663,466	0	G	T	6.6%	231.9 / 4.8	75	V9L (GTG→TTG)	suhB	inositol monophosphatase
Reads supporting (aligned to +/- strand):  ref base G (37/33);  new base T (3/2);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.51e-01

AGCATTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATC  >  NC_000913/2663319‑2663614                                                                                                                                                    |                                                                                                                                                     agctTTGTCTTTCTTCATCTCAAGATTATTCACGCATCTTATCATAAAACGATGACCGATGCCGATCTCGCAGCTCTACTCTGCGCCGTTTTCCCGTTCTTTACCATCCAGTGCGAGAGACCGATGCATCCGATGCTGAACATCGCCGt                                                                                                                                                     <  1:284944/149‑1 (MQ=255) agcaTTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGt                                                                                                                                                     <  8:197392/149‑1 (MQ=255) agcaTTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGAACACGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACAGATGCATCCGATGCTGAACATCGCCGt                                                                                                                                                     >  3:354699/1‑149 (MQ=255) agcaTTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCCATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGATAGACCGATGCATCCGATGCTGAACATCGCCGt                                                                                                                                                     >  6:192927/1‑149 (MQ=255)        tctttcttCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCAATGCATCCGATGCTGAACATCGCCGTgcgcgca                                                                                                                                              <  1:294622/149‑1 (MQ=255)                        gATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGt                                                                                                                               >  4:45299/1‑148 (MQ=255)                               cACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaa                                                                                                                        <  1:28795/148‑1 (MQ=255)                                aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaatt                                                                                                                      >  7:182231/1‑149 (MQ=255)                                aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                       <  3:416/148‑1 (MQ=255)                                aCGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGtaatttaat                                                                                                                       <  3:154242/148‑1 (MQ=255)                                    aTCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGAGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGcc                                                                                                                   >  2:23641/1‑148 (MQ=255)                                          tCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAc                                                                                                             >  7:131619/1‑148 (MQ=255)                                            aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                   <  8:94762/124‑1 (MQ=255)                                            aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                   <  4:141381/124‑1 (MQ=255)                                            aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                   >  3:141381/1‑124 (MQ=255)                                            aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCTTGCGCGCTGCGCGCACGGc                                                                                                                                   >  7:94762/1‑124 (MQ=255)                                            aTAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                          <  1:138147/149‑1 (MQ=255)                                             tAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTAt                                                                                                          >  2:368272/1‑148 (MQ=255)                                                  cGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcc                                                                                                    <  3:34946/149‑1 (MQ=255)                                                   gAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAccc                                                                                                   >  7:90226/1‑149 (MQ=255)                                                    aaGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAcccc                                                                                                  <  1:377845/149‑1 (MQ=255)                                                            tGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGGAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTg                                                                                          >  5:54775/1‑149 (MQ=255)                                                                gATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGa                                                                                      <  8:131619/149‑1 (MQ=255)                                                                gATCTCGCAGCTATACTCTGCGCCGTTTTCCCGTTCATTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGCTAATTTAATTGCCAAAAACTATGAAACCCAGGACCCTCTAg                                                                                       >  3:50918/1‑148 (MQ=255)                                                                      gctgctATAATCTGCGCCGTGTTCCCGTTCTTTAACATCCAGTGCGAGAGACCGATGCATCAGATGCTGAACATCGCCGGGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCTAAAAACTATGAAACCCCGGACGCTGTAGAAGCGa                                                                                 <  5:192927/148‑1 (MQ=255)                                                                           tataCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAg                                                                            >  7:235088/1‑148 (MQ=255)                                                                             taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaaa                                                                         <  1:347292/149‑1 (MQ=255)                                                                             taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGaaa                                                                         <  8:69534/149‑1 (MQ=255)                                                                             taCTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGa                                                                           >  3:103736/1‑147 (MQ=255)                                                                                tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAgg                                                                       >  4:266760/1‑148 (MQ=255)                                                                                tctGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCAAGAAAgg                                                                       >  6:344506/1‑148 (MQ=255)                                                                                  tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAc                                                                                                     >  8:337740/1‑116 (MQ=255)                                                                                  tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAAc                                                                                                     <  7:337740/116‑1 (MQ=255)                                                                                  tGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGc                                                                      <  8:182231/147‑1 (MQ=255)                                                                                  tGCGCCGTTTTCACGTTCTTTAACATCCAGTGAGAGATACCGATGCATCCGATGCTGAACATCGCCGTGCGAGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGACAGAAAGCCAg                                                                    >  1:122004/1‑149 (MQ=255)                                                                                    cgcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCTATGCATCCTATGCTGAACATCGCCGTGCGCGCAGCGAGCAAGGCGGGTAATTTAATTGCCAAATACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCa                                                                     <  1:209202/146‑1 (MQ=255)                                                                                     gcCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTa                                                                                                           >  4:123753/1‑107 (MQ=255)                                                                                             cccGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAACACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCg                                                         >  7:301603/1‑149 (MQ=255)                                                                                                gTTCTTTAACATCCAGTGAGAGAGACCGATTCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                      <  7:177032/149‑1 (MQ=255)                                                                                                        aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGt                                               <  1:116680/148‑1 (MQ=255)                                                                                                        aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                              >  1:379206/1‑149 (MQ=255)                                                                                                        aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                              <  7:339060/149‑1 (MQ=255)                                                                                                        aCATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                              <  1:23641/149‑1 (MQ=255)                                                                                                          aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGa                                                                                                        >  5:92670/1‑89 (MQ=255)                                                                                                          aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGa                                                                                                        <  6:92670/89‑1 (MQ=255)                                                                                                          aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCATAAAGTCAGTAACGATTTCGTGACCAACGTAg                                             >  5:497/1‑148 (MQ=255)                                                                                                          aTCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCCGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACg                                                >  7:57633/1‑145 (MQ=255)                                                                                                           tCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGa                                            >  1:296245/1‑148 (MQ=255)                                                                                                                tGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAATATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGc                                      <  3:45299/149‑1 (MQ=255)                                                                                                                  agagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCt                                     <  3:266760/148‑1 (MQ=255)                                                                                                                  agagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTg                                    >  7:81551/1‑149 (MQ=255)                                                                                                                  agagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTAGTGACCAACGTAGATAAAGCTg                                    >  8:134701/1‑149 (MQ=255)                                                                                                                   gagagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGc                                   >  7:199085/1‑149 (MQ=255)                                                                                                                    agagaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGCCGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGAGTTCGTGACCACCGTAGATAAGGCTGcc                                  >  8:170554/1‑149 (MQ=255)                                                                                                                      agaCCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCg                                 >  6:174307/1‑148 (MQ=255)                                                                                                                          cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGACACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCg                            >  2:109174/1‑149 (MQ=255)                                                                                                                          cGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTTGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCg                            >  4:70776/1‑149 (MQ=255)                                                                                                                            aTGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGt                          >  8:38576/1‑149 (MQ=255)                                                                                                                                aTCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTATTTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                      <  4:166889/117‑1 (MQ=255)                                                                                                                                aTCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                      >  3:166889/1‑117 (MQ=255)                                                                                                                                 tCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTa                     >  7:283006/1‑149 (MQ=255)                                                                                                                                   cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTAt                    <  7:57289/148‑1 (MQ=255)                                                                                                                                   cGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATc                   <  8:90226/149‑1 (MQ=255)                                                                                                                                         tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                              <  5:338151/116‑1 (MQ=255)                                                                                                                                         tGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTa                                              >  6:338151/1‑116 (MQ=255)                                                                                                                                          gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                      >  8:56203/1‑107 (MQ=255)                                                                                                                                          gAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGa                                                      <  7:56203/107‑1 (MQ=255)                                                                                                                                           aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTATTTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGATAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa            <  4:354699/148‑1 (MQ=255)                                                                                                                                           aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTa                     <  1:217977/139‑1 (MQ=255)                                                                                                                                           aaCATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGa            <  2:296245/148‑1 (MQ=255)                                                                                                                                           aaCATCGCCGTGCCCGCAGCGCGCAAGGAGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCAGTATAACCCAGCCAGAAAGGCACTAACGAATTCGTGAAAAACGTAGATAAAGCTGACGAAGCGTTGATTATAGatcgcaa           >  2:217977/1‑143 (MQ=255)                                                                                                                                             cATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTc         <  8:235088/149‑1 (MQ=255)                                                                                                                                              aTCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAAATATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCg        >  4:285905/1‑149 (MQ=255)                                                                                                                                               ccGCCTTGCGCGCTGCGCGCACGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAATGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCg        >  4:289596/2‑148 (MQ=255)                                                                                                                                                 gCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTTGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTa      <  3:70776/148‑1 (MQ=255)                                                                                                                                                    gTGCGCGCAGCGCGCAAGGCGGGTAAGTTAATTGCCATAAACTATGAAACCCCGGACGCGGGAGAACCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGAGTCGTAAATc  >  5:110559/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     AGCATTGTCTTTCTTCATCTAAAGATTATTCACGCATCTTATCATAAAACGAAGACAGATGCCGATCTCGCTGCTATACTCTGCGCCGTTTTCCCGTTCTTTAACATCCAGTGAGAGAGACCGATGCATCCGATGCTGAACATCGCCGTGCGCGCAGCGCGCAAGGCGGGTAATTTAATTGCCAAAAACTATGAAACCCCGGACGCTGTAGAAGCGAGCCAGAAAGGCAGTAACGATTTCGTGACCAACGTAGATAAAGCTGCCGAAGCGGTGATTATCGACACGATTCGTAAATC  >  NC_000913/2663319‑2663614

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