Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,760	0	T	C	23.9%	137.6 / 22.6	67	K250K (AAA→AAG)	citG	triphosphoribosyl‑dephospho‑CoA synthase
Reads supporting (aligned to +/- strand):  ref base T (2/49);  new base C (16/0);  total (18/49)
Fisher's exact test for biased strand distribution p-value = 1.36e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646622‑646896                                                                                                                                            |                                                                                                                                         cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCtt                                                                                                                                         >  1:302404/1‑140 (MQ=255) cTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCtt                                                                                                                                         >  1:612567/1‑140 (MQ=255)  ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCttt                                                                                                                                        <  1:13997/140‑1 (MQ=255)  ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCttt                                                                                                                                        <  2:800294/140‑1 (MQ=255)    aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCttttt                                                                                                                                      <  2:743373/140‑1 (MQ=255)    aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCttttt                                                                                                                                      <  2:56359/140‑1 (MQ=255)          ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCaata                                                                                                                                <  1:34682/140‑1 (MQ=255)           tAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GTCCCCCTTTTTCaataa                                                                                                                               >  2:390898/1‑140 (MQ=255)                  gtgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAT‑GTCCCCCTTTTTCaataaaattttt                                                                                                                        >  2:79112/1‑133 (MQ=255)                   tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAT‑GCCCCCCTTTTTCAATAATTTTTtcg                                                                                                                       >  1:823213/1‑140 (MQ=255)                   tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GCCCCCCTTTTTTAATAAAGTTTgtg                                                                                                                       >  2:127708/1‑138 (MQ=255)                    gCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GCCCCCCTTTTTTAATAATGTTTtcgc                                                                                                                      >  2:510944/1‑140 (MQ=255)                     cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGcc                                                                                                                     <  1:510944/140‑1 (MQ=255)                     cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGcc                                                                                                                     <  1:182110/140‑1 (MQ=255)                     cTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GCCCCCCTTTTTTAATAATGTTTTCGcc                                                                                                                     >  1:905038/1‑140 (MQ=255)                       aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                                   <  1:532536/140‑1 (MQ=255)                       aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                                   <  1:196095/140‑1 (MQ=255)                       aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTc                                                                                                                   <  2:612567/140‑1 (MQ=255)                         aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTcgc                                                                                                                 <  2:479894/140‑1 (MQ=255)                            ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCt                                                                                                              <  2:324923/140‑1 (MQ=255)                                 tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGc                                                                                                         <  2:302404/140‑1 (MQ=255)                                 tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGc                                                                                                         <  1:127708/140‑1 (MQ=255)                                 tAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGc                                                                                                         <  1:104769/140‑1 (MQ=255)                                     gATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgc                                                                                                     <  1:991074/140‑1 (MQ=255)                                          gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAT‑GCCCCCCTTTTTCAAAAAATTTTTCGCCCCCCCCCGTGGCCAGCGCCgg                                                                                                >  2:1359049/1‑140 (MQ=255)                                            aGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcc                                                                                              <  1:563120/140‑1 (MQ=255)                                              gTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGcccc                                                                                            <  1:94973/140‑1 (MQ=255)                                                cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAT‑GCCCCCTTTTTTTAATAATATTTTCGGCTCCCGCTGTTTGCAGCGCAGGaccccc                                                                                          >  2:1132605/1‑140 (MQ=255)                                                cAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAT‑GCCCCCCTTTTTTTATAATGTTTTCGCCCCCCCCCGTAGGCCCCGCCCGcccccc                                                                                          >  2:885774/1‑140 (MQ=255)                                                   cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAT‑GCCCCCCTTTTTTAATAATATTTTCGCCCCCCCCCGTAGCCACCGCCGGccccccccc                                                                                       >  2:1084431/1‑137 (MQ=255)                                                   cACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGCGGGTCGGATGGCCCCCTTTTTTC‑ATAAATTTTTTGCCTCCCGCTGTAGCCCACGCCGGccccccccc                                                                                       >  2:1327702/1‑137 (MQ=255)                                                     cTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGGTGGAT‑GCCCCCCTTTTTCAATAATTTTTTCGCCCCGCGCTGGTGCCAGCGCGGGccccccccccc                                                                                     >  2:1116872/1‑135 (MQ=255)                                                           ccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACcgcg                                                                               <  2:99501/140‑1 (MQ=255)                                                             tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTGGGAT‑GCCCCCCCTTTTTCATAATTTTTTCGCCCCGCGCCGTAGCCCGCGCCGGCCCCCCCCGCccccccccc                                                                             >  2:404478/1‑135 (MQ=255)                                                                aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                                          <  2:1079350/140‑1 (MQ=255)                                                                aCTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATg                                                                          <  2:1071485/140‑1 (MQ=255)                                                                 aggagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCTTCGCCACCGCGCGATGc                                                                         <  1:429308/138‑1 (MQ=255)                                                                            ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCGCCCCTCGCCTCCGCGCGATGCAACGTTGGTAt                                                              <  1:51485/140‑1 (MQ=255)                                                                            ttcgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTAt                                                              <  1:653053/140‑1 (MQ=255)                                                                              cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAT‑GCCCCCCTTTTTTAATAATGTTTTCGCCCCCCGCTGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCCCTGCCACGGTGGTGTCg                                                            >  2:135959/1‑140 (MQ=255)                                                                              cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTTGGAT‑GCCCCCCTTTTTCAAAAATTTTTTTGCCTCCCGCCGGAGGCAGCGCCGGCCCCCCCCCCCCCCCCCCGATGCAACGTTGGTATcc                                                            >  2:546828/1‑139 (MQ=255)                                                                              cgATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAT‑GCCCCCCTTTTTTAATAATTTTTTTGCCTCCCGCTGTAGCCCGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAACGGTGGTAttt                                                            >  2:429308/1‑138 (MQ=255)                                                                                  acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGt                                                        <  1:79112/140‑1 (MQ=255)                                                                                      tCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGAt                                                    <  1:1359049/140‑1 (MQ=255)                                                                                        gcTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCg                                                  <  1:74309/139‑1 (MQ=255)                                                                                        ccTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GCCCCCCTTTTTCAATAATTTTTTCGCCCCCCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGTTGGTATCGCCGCTGATcc                                                  >  2:416531/1‑139 (MQ=255)                                                                                           gTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GCCCCCCTTTTTCAATAAATTTTTTGCCTCCCGCTGTGGCCAGCGCCGGCCCCCCCCCCCCCCCCGCGCTGCCACGGTGGGATTGCCGTTGATTGccc                                               >  2:570891/1‑139 (MQ=255)                                                                                                caggCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                          <  1:864736/137‑1 (MQ=255)                                                                                                aCTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGt                                          <  2:602390/140‑1 (MQ=255)                                                                                                  tGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAg                                        <  1:926078/140‑1 (MQ=255)                                                                                                  tGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAg                                        <  1:1116872/140‑1 (MQ=255)                                                                                                  tGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAg                                        <  1:1084431/140‑1 (MQ=255)                                                                                                  tGCCGGAGATAATCGAGATCGGCGGGGGGTCGGAT‑GCCCCCCTTTTTCAATAATATTTTCGCCCCGCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGCTGCAACCGTGGTATCGCCGGTGGTCGCCCTCAGTAg                                        >  2:1388012/1‑140 (MQ=255)                                                                                                  tGCCGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GCCCCCCTTTTTTAATAAAATTTTCGCCCCCCGCCGGAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGCGGCAACGTTGGGATCGCCCCTGATCGCCCTCCGTAg                                        >  2:684649/1‑140 (MQ=255)                                                                                                    ccGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGa                                      <  1:130594/140‑1 (MQ=255)                                                                                                    ccGGAGATAATCGAGATCGGCGGGGGGTTGGAT‑GCCCCCCTTTTTTAATAATGTTTTCGCCTCGCGCCGTTGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGGTGCAACGTTGGGTTCGCCGTTGGTCCCCCTCaggaggg                                      >  2:74309/1‑139 (MQ=255)                                                                                                    ccGGAGATAATCGAGATCGGCGGGGGGTCGGAT‑GCCCCCCTTTTTCAATAATTTTTTTGCCTCCCGCCGTAGCCCGCGCAGGCCCCCCCCCCCCCCCCCCGGTGCAACGTTGGTATCGCCGCTGATCGCCCTCCGTAggg                                      >  2:281136/1‑139 (MQ=255)                                                                                                    ccGGAGATAATCGAGATCGGCGGGGGGTCGAAT‑GCCCCCCTTTTGTAATAATTTTTTCGCCCCGCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCGCGATGCAACGGTGGGATCGCCGGTGATCGCCCTCAGTAGGa                                      >  2:1040737/1‑140 (MQ=255)                                                                                                      ggAGATAATCGAGATCGGCGGGGGGTTGGAT‑GCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCCGTAGGCAGCGCAGGCCCCCCCCCCCCCCCCCCGATGCCAAGTTGGGATCGCCCTTGATCGCCCTCAGTAGGAgg                                    >  2:472353/1‑139 (MQ=255)                                                                                                       gagaTAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                                   <  1:1327702/140‑1 (MQ=255)                                                                                                       gagaTAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                                   <  1:1132605/140‑1 (MQ=255)                                                                                                       gagaTAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                                   <  2:1299021/140‑1 (MQ=255)                                                                                                       gagaTAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCa                                   <  1:1138117/140‑1 (MQ=255)                                                                                                            aaTCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTa                              <  2:1018664/140‑1 (MQ=255)                                                                                                             aTCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCGCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                             <  1:1202976/140‑1 (MQ=255)                                                                                                             aTCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                             <  2:1497915/140‑1 (MQ=255)                                                                                                             aTCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                             <  1:570891/140‑1 (MQ=255)                                                                                                              tCGAGATCGGCGGGGGGTTGGAT‑GCCCCCCTTTTTCAATAATATTTTTGCCTCCCGCTGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGCTGCAACGTTGGTATCGCCGTTGATCGCCCTCAGTAGGGGCAAgggctg                            >  2:864736/1‑136 (MQ=255)                                                                                                              tCGAGATCGGCGGGGGGTTGAAT‑GCCCCCCTTTTTCAATAATGTTTTCGCCTCCCGCTGTAGCCAGCGCCGGCCCCCCCCCCCCCCCCCCGGTGCAACGTTGGTATCGCCGCTGATCGCCCTCAGTAGGAGGAAGGGATc                            >  2:130594/1‑140 (MQ=255)                                                                                                                gagaTCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                          <  2:994029/140‑1 (MQ=255)                                                                                                                       gcggGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc                   <  1:281136/140‑1 (MQ=255)                                                                                                                       gcggGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGc                   <  1:96166/140‑1 (MQ=255)                                                                                                                        cggGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGcc                  <  1:546828/140‑1 (MQ=255)                                                                                                                            ggTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt              <  1:472353/140‑1 (MQ=255)                                                                                                                             gTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTc             <  1:885774/140‑1 (MQ=255)                                                                                                                              ttCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa            <  1:1017983/140‑1 (MQ=255)                                                                                                                              ttCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCa            <  1:699131/140‑1 (MQ=255)                                                                                                                                 gAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGa         >  1:1249163/1‑140 (MQ=255)                                                                                                                                    t‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      <  1:1087081/140‑1 (MQ=255)                                                                                                                                    t‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      <  2:823213/140‑1 (MQ=255)                                                                                                                                    t‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt      <  1:417275/140‑1 (MQ=255)                                                                                                                                         cccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  <  1:1040737/140‑1 (MQ=255)                                                                                                                                         cccTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAAcc  >  2:1457560/1‑140 (MQ=255)                                                                                                                                            |                                                                                                                                         CTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGCGGGGGTTCGAAT‑GCCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTTGGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAACC  >  NC_000913/646622‑646896

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