Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,926,366	0	A	C	21.3%	145.7 / 13.3	61	G30G (GGT→GGG)	mioC	flavoprotein MioC
Reads supporting (aligned to +/- strand):  ref base A (36/12);  new base C (0/13);  total (36/25)
Fisher's exact test for biased strand distribution p-value = 7.92e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACGCGTTG  >  NC_000913/3926228‑3926505                                                                                                                                            |                                                                                                                                            cTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAAc                                                                                                                                            >  1:665536/1‑140 (MQ=255)      ccTGCAATGCTTCATATAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGcc                                                                                                                                       >  1:761521/1‑140 (MQ=255)      ccTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGcc                                                                                                                                       <  1:835626/140‑1 (MQ=255)       cTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCct                                                                                                                                      >  2:134512/1‑140 (MQ=255)              gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAg                                                                                                                               >  1:72719/1‑140 (MQ=255)              gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAg                                                                                                                               >  2:385580/1‑140 (MQ=255)              gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCCTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAg                                                                                                                               >  2:149432/1‑140 (MQ=255)                  cATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCttt                                                                                                                           >  2:89669/1‑140 (MQ=255)                   aTAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGGGGAGCTGATAACCAGCCAGATCCCGGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCTTTTCGGTGGAAAAACCCGCCTCTTCCAGCtttt                                                                                                                          <  2:911224/140‑1 (MQ=255)                          ggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGCCCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAg                                                                                                                   <  2:494017/140‑1 (MQ=255)                          ggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAg                                                                                                                   >  1:1155010/1‑140 (MQ=255)                           gagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAgg                                                                                                                  >  2:1098377/1‑140 (MQ=255)                             aaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtg                                                                                                                <  2:1084391/139‑1 (MQ=255)                              aaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtgt                                                                                                               >  1:1016226/1‑140 (MQ=255)                               aatgTTGTCCGGAATATCTCCGGCCCCGTGGGTGGAGCTGATACCCAGCCAGACCCGTGAGCCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCTTTTGGGGGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt                                                                                                              <  2:1155010/137‑1 (MQ=255)                               aaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATACCCAGCCAGACCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt                                                                                                              <  2:352519/140‑1 (MQ=255)                               aaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt                                                                                                              <  1:134512/140‑1 (MQ=255)                               aaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt                                                                                                              <  1:677222/140‑1 (MQ=255)                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                        >  2:652328/1‑140 (MQ=255)                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                        >  2:766393/1‑140 (MQ=255)                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                        >  2:788680/1‑140 (MQ=255)                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                        >  2:798188/1‑140 (MQ=255)                                     gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGAGTTCTGCTa                                                                                                        >  1:439090/1‑140 (MQ=255)                                      tCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGACCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGCCCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTAc                                                                                                       <  2:205292/140‑1 (MQ=255)                                               tctcGGGCACGGGGGGTGGAGCTGAAACCCAGCCAGATCCCGGAGGCAGGTAAATCTTCAAACAGC‑GGCCCGTGCAGCTTTCGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGc                                                                                              <  2:133132/140‑1 (MQ=255)                                                ctcCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCa                                                                                             >  1:674999/1‑140 (MQ=255)                                                ctcCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCa                                                                                             >  2:536586/1‑140 (MQ=255)                                                       aCCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGa                                                                                      >  2:353003/1‑140 (MQ=255)                                                               gggAGCTGAAACCCAGCCACACCCTGGAGGAAGTTAAATCTTTAAACAGC‑GGCCGGTGGGGTTTTGGGGAGATAAATCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTCCTGCTACATATTCGGCACCGCCGAGGGtgctg                                                                              <  2:630497/139‑1 (MQ=255)                                                                 gagcgaaTACCCACCCAGTTCCGTGAGGCAGTTAAATCTTCAACCGGC‑GCCCGGTGCGGCTTTTGGGGGTTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCCCCGCCGAGGGTGCTGcc                                                                            <  2:439090/134‑1 (MQ=255)                                                                 gAGCTAAAACCCAGCCAGATCCTGGAGGAAGGTAATTCTTCT‑ACAGCGGGCCCGGGCGGGGTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGcc                                                                            <  2:656282/140‑1 (MQ=255)                                                                   gCTGAAACCCACCCAGACCCTTGAGGCAGGAAAATTTTCAAACGGC‑GGACCGTGCGGCGTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGc                                                                          <  2:260227/140‑1 (MQ=255)                                                                       aaaaCCAGCCAGATCCCGGAGGAAGTTAAATCTTCTAACAGC‑GGCCCGTGCGGGTTTTGGGGGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTTGGCACCGCCGAGGGTGCTGCCGCTGAt                                                                      <  2:127580/138‑1 (MQ=255)                                                                       aTACCCAGCCAGACCCCTGAGGCAGGTAAATCTTCAAACAGC‑GGACCGTGCAGCGTTTCGGTGGAAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                                      <  2:1088365/140‑1 (MQ=255)                                                                       aTAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                                      >  2:50465/1‑140 (MQ=255)                                                                        aaaCCAGCCAGCCCCCGGAGGCAGTAAAATTTTCAAACAGC‑GGCCGGTGCGGCTTTTGGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATa                                                                     <  2:761521/139‑1 (MQ=255)                                                                          aCCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAag                                                                   <  1:353003/140‑1 (MQ=255)                                                                            cagccagcTCCCGGAGGCAGGAAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAagag                                                                 <  2:225412/140‑1 (MQ=255)                                                                            cagccagATTCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAagag                                                                 >  1:965051/1‑140 (MQ=255)                                                                             agccagATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGt                                                                >  1:420913/1‑140 (MQ=255)                                                                              gccagATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTg                                                               >  1:225412/1‑140 (MQ=255)                                                                                   ctCCCGGAGGCAGGAAATTCTTCAAACGGC‑GGACGGTGCAGTGTTTGGGTGTAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATc                                                          <  2:674999/139‑1 (MQ=255)                                                                                   aTCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATc                                                          >  2:387520/1‑140 (MQ=255)                                                                                              aGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAAcc                                               >  2:918292/1‑140 (MQ=255)                                                                                               ggTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCg                                              >  1:117277/1‑140 (MQ=255)                                                                                               ggTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCg                                              >  1:238455/1‑140 (MQ=255)                                                                                                gTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                             >  1:273977/1‑140 (MQ=255)                                                                                                aaaattcttcaAACAGC‑GGCCCGTGCGGCTTTTCGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                             <  2:420913/135‑1 (MQ=255)                                                                                                  aattcttcaAACGGG‑GGCCGGTGCGGCTTTTGGGTGAAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                           <  2:665536/137‑1 (MQ=255)                                                                                                  aaaTCTTCTAACAGG‑GGCCCGGGCAGCGTTTGGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                           <  1:513775/140‑1 (MQ=255)                                                                                                  aaaTCTTCAAACGGC‑GGACCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                           <  1:865272/140‑1 (MQ=255)                                                                                                   atttttctAACCGC‑GGCCCGGGGAGCTTTTGGGTGTAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCAc                                          <  2:411089/136‑1 (MQ=255)                                                                                                    aTCTCCAAACAGG‑GGACCGGGCAGCGTTGCGGGGGTAAAGCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAAAAGAGTGATATCTGCCATAAACCGCCAcc                                         <  1:62271/140‑1 (MQ=255)                                                                                                           aaCAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCAGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatta                                  >  2:58467/1‑140 (MQ=255)                                                                                                             acgggggACGGTGCGGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTCCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTaa                                 <  2:72719/135‑1 (MQ=255)                                                                                                              ggC‑GGACCGTGCGGCGTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAaga                               <  1:975647/139‑1 (MQ=255)                                                                                                               gC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAagag                              >  2:151566/1‑140 (MQ=255)                                                                                                               gC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAagag                              >  1:211035/1‑140 (MQ=255)                                                                                                               aggggcccGTGCGGTTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAaga                               <  2:496660/134‑1 (MQ=255)                                                                                                                c‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGt                             >  2:879988/1‑140 (MQ=255)                                                                                                                c‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGt                             >  2:432820/1‑140 (MQ=255)                                                                                                                  ggACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTg                            >  1:451849/1‑140 (MQ=255)                                                                                                                   gACGGTGCGGCTTTTCGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTgg                           <  1:149432/140‑1 (MQ=255)                                                                                                                   gACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTgg                           >  1:276437/1‑140 (MQ=255)                                                                                                                         gCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTAtt                     >  1:698048/1‑140 (MQ=255)                                                                                                                           aGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGt                   <  1:918292/140‑1 (MQ=255)                                                                                                                            gCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTa                  <  1:536586/140‑1 (MQ=255)                                                                                                                              gTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACg                >  1:351038/1‑140 (MQ=255)                                                                                                                                  gggtgtaaaacCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCtgt            <  2:1016226/133‑1 (MQ=255)                                                                                                                                      ggtAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAc        <  1:50465/140‑1 (MQ=255)                                                                                                                                       gtAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAcg       >  2:635038/1‑140 (MQ=255)                                                                                                                                         aaaaCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAcgcg     >  1:30748/1‑140 (MQ=255)                                                                                                                                         aaaaCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAcgcg     <  1:1098377/140‑1 (MQ=255)                                                                                                                                          aaaCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAATGCGt    <  1:766393/140‑1 (MQ=255)                                                                                                                                            aCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACGCGTTg  <  1:798188/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                            CTGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCTGAGGCAGGTAAATCTTCTAACAGC‑GGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACGCGTTG  >  NC_000913/3926228‑3926505

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