Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,590,908	(C)5→7	coding (372/2151 nt)	btsT ←	pyruvate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,590,903	1	.	C	90.3%	201.2 / 13.5	62	coding (377/2151 nt)	btsT	pyruvate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base . (0/6);  new base C (37/19);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 2.88e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
*	NC_000913	4,590,903	2	.	C	87.1%	193.4 / 21.1	62	coding (377/2151 nt)	btsT	pyruvate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base . (0/8);  new base C (37/17);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 3.20e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

TAAGAAACAGCCAAACAGCGCGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑A‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTGGTCGGAAC  >  NC_000913/4590767‑4591033                                                                                                                                            ||                                                                                                                                   tAAGAAACAGCCAAACAGCGCGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑A‑‑ccc                                                                                                                                 >  1:619882/1‑140 (MQ=255)   aGAAACAGCCAAACAGCGCGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑A‑‑ccccc                                                                                                                               >  2:946823/1‑140 (MQ=255)   aGAAACAGCCAAACAGCGCGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑A‑‑ccccc                                                                                                                               >  1:324580/1‑140 (MQ=255)     aaaCAGCCAAACAGCGCGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑A‑‑ccccccc                                                                                                                             >  1:178438/1‑138 (MQ=255)              aaCAGCGCGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACGCCC‑‑CCCCCGCCagcagc                                                                                                                      <  2:178438/140‑1 (MQ=255)                  gcgcGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACGCCC‑‑CCCCCGCCAGCAGCcaca                                                                                                                  <  2:357693/140‑1 (MQ=255)                  gcgcGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCcaca                                                                                                                  <  2:212506/140‑1 (MQ=255)                  gcgcGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCcaca                                                                                                                  <  2:753343/140‑1 (MQ=255)                   cgcgATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAg                                                                                                                 <  2:347391/140‑1 (MQ=255)                   cgcgATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAg                                                                                                                 >  1:347391/1‑140 (MQ=255)                     cgATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCg                                                                                                               >  1:375956/1‑140 (MQ=255)                         aGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑AC‑CCCCCGCCAGCAGCCACAGCGTGCca                                                                                                          <  1:946823/140‑1 (MQ=255)                             cccGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCcaggc                                                                                                       >  1:766421/1‑140 (MQ=255)                             cccGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCcaggc                                                                                                       >  1:465350/1‑140 (MQ=255)                              ccGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCcaggca                                                                                                      <  1:310769/140‑1 (MQ=255)                                       gTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCcaggca                                                                                                      <  1:673378/131‑1 (MQ=255)                                       gTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCcaggca                                                                                                      >  2:673378/1‑131 (MQ=255)                                       gTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCa                                                                                             >  1:143262/1‑140 (MQ=255)                                       gTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCa                                                                                             >  1:17756/1‑140 (MQ=255)                                       gTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCa                                                                                             >  2:502851/1‑140 (MQ=255)                                            aTCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGTAGGTAGCCCATCTgc                                                                                        >  1:767860/1‑140 (MQ=255)                                            aTCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTgc                                                                                        >  1:662955/1‑140 (MQ=255)                                                        atcatcTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAAcc                                                                            >  1:468339/1‑140 (MQ=255)                                                         tcatcTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCg                                                                           >  2:975698/1‑140 (MQ=255)                                                             ctcACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGAcc                                                                       >  2:897504/1‑140 (MQ=255)                                                                aCCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACccac                                                                    >  2:210354/1‑140 (MQ=255)                                                                  cAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACccacca                                                                  >  1:1009115/1‑140 (MQ=255)                                                                   aaGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAg                                                                 <  1:664328/140‑1 (MQ=255)                                                                   aaGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAg                                                                 <  1:975698/140‑1 (MQ=255)                                                                       gaTGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGa                                                             >  2:634076/1‑140 (MQ=255)                                                                          gcgcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACGCCC‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGAc                                                            <  2:486971/138‑1 (MQ=255)                                                                          gcgcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGAc                                                            >  1:486971/1‑138 (MQ=255)                                                                          gcgcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCa                                                          >  2:163404/1‑140 (MQ=255)                                                                          gcgcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCa                                                          >  1:149177/1‑140 (MQ=255)                                                                          gcgcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCa                                                          >  1:129759/1‑140 (MQ=255)                                                                          gcgcCACTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCa                                                          >  1:837986/1‑140 (MQ=255)                                                                            gcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGc                                                        <  2:662955/140‑1 (MQ=255)                                                                            gcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGc                                                        >  2:30283/1‑140 (MQ=255)                                                                            gcCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGc                                                        >  1:356191/1‑140 (MQ=255)                                                                            cccATTACGGCGAGAGGAGAAAAACAGCCCCATAAAGTCCTGAACCGCACCGGCCAGCACGCCC‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGTAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGc                                                        <  2:767860/139‑1 (MQ=255)                                                                             ccATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACGCCC‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCa                                                       <  2:375956/140‑1 (MQ=255)                                                                                ttACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCg                                                    >  1:335034/1‑140 (MQ=255)                                                                                     gCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCCTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAAACGGACCCACCAGCGGACCCGCACCGGCGaa                                               >  2:980450/1‑139 (MQ=255)                                                                                       gagaGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAg                                             >  2:41222/1‑140 (MQ=255)                                                                                        agagGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACGCCC‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAg                                             <  2:420484/139‑1 (MQ=255)                                                                                        agagGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAg                                             >  1:420484/1‑139 (MQ=255)                                                                                            gagaTAAACAGCCCCATAAAGCCCTGAACCGCCCCGGCCAGCACGCCC‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAgcggcg                                        <  2:149177/140‑1 (MQ=255)                                                                                            gagaTAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAgcggcg                                        <  2:356191/140‑1 (MQ=255)                                                                                                  aaCAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAgtgg                                  >  2:63156/1‑140 (MQ=255)                                                                                                     aGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGa                               <  1:897504/140‑1 (MQ=255)                                                                                                      gccccATAAAGTCCTGACCCGCACCGGCCAGCCCGCCC‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGAc                              <  2:129759/137‑1 (MQ=255)                                                                                                      gCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGAc                              >  1:831274/1‑140 (MQ=255)                                                                                                       cACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGAcc                             >  2:453755/1‑140 (MQ=255)                                                                                                       cACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGAcc                             >  1:322721/1‑140 (MQ=255)                                                                                                       cACCATAAAGTCCTGAACCGCACCGACCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGAcc                             >  1:305373/1‑140 (MQ=255)                                                                                                        aCCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCa                            >  1:292059/1‑140 (MQ=255)                                                                                                         ccATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCaa                           <  2:831274/140‑1 (MQ=255)                                                                                                           aTAAAGTCCTGAACCGCACCGGCCAGCACG‑‑AC‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAaaca                        <  1:666854/140‑1 (MQ=255)                                                                                                            tAAAGTCCTGAACCGCACCGACCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAaaca                        <  2:305373/140‑1 (MQ=255)                                                                                                             aaaGTCCTGACCCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAaacaa                       <  2:766421/140‑1 (MQ=255)                                                                                                             aaaGTCCTGAACCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAaacaa                       <  2:143262/140‑1 (MQ=255)                                                                                                                      aaCCGCACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACg              >  2:41208/1‑140 (MQ=255)                                                                                                                           cACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTgg         >  2:32062/1‑140 (MQ=255)                                                                                                                           cACCGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTgg         >  2:219024/1‑140 (MQ=255)                                                                                                                              cGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTGGTCg      <  1:32062/140‑1 (MQ=255)                                                                                                                              cGGCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTGGTCg      >  1:127331/1‑140 (MQ=255)                                                                                                                                gCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTGGTCGGa    <  2:468339/140‑1 (MQ=255)                                                                                                                                gCCAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTGGTCGGa    <  2:17756/140‑1 (MQ=255)                                                                                                                                  cAGCACG‑‑ACCCCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTGGTCGGAAc  <  2:335034/140‑1 (MQ=255)                                                                                                                                            ||                                                                                                                                   TAAGAAACAGCCAAACAGCGCGATAGTCCCCGGTACTGGTCCCATCTCTTCTTTGATCATCTCACCAAGAGATGCGCCATTACGGCGAGAGGAGATAAACAGCACCATAAAGTCCTGAACCGCACCGGCCAGCACG‑‑A‑‑CCCCCGCCAGCAGCCACAGCGTGCCAGGCAGGTAGCCCATCTGCGCGGCGAGAACCGGACCCACCAGCGGACCAGCACCGGCGATAGCGGCGAAGTGGTGACCAAACAACACGTAACGGTTGGTCGGAAC  >  NC_000913/4590767‑4591033

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