Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,147,792	T→A	S58T (TCG→ACG)	plsX →	putative phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,147,792	0	T	A	97.3%	109.4 / NA	37	S58T (TCG→ACG)	plsX	putative phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (14/22);  minor base C (1/0);  total (15/22)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

CATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACC  >  NC_000913/1147546‑1148024                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         cATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTAcgcg                                                                                                                                                                                                                                      >  1:222161/1‑251 (MQ=255)   tCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTc                                                                                                                                                                                                                                    >  2:480051/1‑251 (MQ=255)      gCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGTAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTg                                                                                                                                                                                                                                  <  1:33620/250‑1 (MQ=255)              aaGCTTATTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGGTGCCTTGCCACGCCTAGCCCCGGCGTTAGATGTCATGGGAGGGGATGTTGGCCCTTCCGTGAAAGTGCCTGCACCATTGCAGGCACTGAATTCTAATTCGTAACTCACTATTCTTTTCGTCGGCATTTCCGTCGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAtattat                                                                                                                                                                                                                          <  2:24955/250‑1 (MQ=255)                                   cAGGCAACAGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCTTGAGATGTCATGGGAGGGGATTGTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATACGCAACTCACTCTGCTTTTAGTCGGCAATTCCGACGCCAGCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGcc                                                                                                                                                                                                    <  2:133939/251‑1 (MQ=255)                                     ggcaccgggggAAAGACCAAACCGGGCGGGGACGATACCTTGACACGCCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACCCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGGCAGTTATCGc                                                                                                                                                                                                     <  2:607390/241‑1 (MQ=255)                                     ggCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTATTCTTTTAGTCGGCATTTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGc                                                                                                                                                                                                     >  1:607390/1‑248 (MQ=255)                                         ccgggggAAAGACCAAACCGGTCGGTGACGATACCTTGACTCGTCTAACCCTGGCGGTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGAAAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGCCGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGAt                                                                                                                                                                                              <  2:535478/248‑1 (MQ=255)                                             gggAAAGACAAAACCGGGCGGCGACGATACCTTGACCCGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGccc                                                                                                                                                                                          <  1:472098/251‑1 (MQ=255)                                                 aaGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGcc                                                                                                                                                                                      <  1:427220/251‑1 (MQ=255)                                                                ggcgACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCcg                                                                                                                                                                       >  2:37191/1‑251 (MQ=255)                                                                      gttACCTTCACACGGCTACCCCTGGCGTTAGATGTAATTGGAGGGTATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAg                                                                                                                                                                 <  2:141879/249‑1 (MQ=255)                                                                               acacGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTggg                                                                                                                                                          >  1:55249/1‑249 (MQ=255)                                                                               acacGTCTAACCCTAGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGAAAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTggg                                                                                                                                                          >  1:176475/1‑249 (MQ=255)                                                                                      tAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATg                                                                                                                                                 >  1:442185/1‑251 (MQ=255)                                                                                                     aTGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTgg                                                                                                                                    >  1:157487/1‑249 (MQ=255)                                                                                                                            ccTTCCGTGAACGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGATGGTGAAAGAATGTCGAGCGCa                                                                                                           <  2:419539/251‑1 (MQ=255)                                                                                                                               tCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGAAAGc                                                                                                        >  2:559101/1‑251 (MQ=255)                                                                                                                               tCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCAGCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAg                                                                                                         <  2:236536/250‑1 (MQ=255)                                                                                                                                       aGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCa                                                                                                <  1:377431/251‑1 (MQ=255)                                                                                                                                         tGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGt                                                                                              >  1:35863/1‑251 (MQ=255)                                                                                                                                            cTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCAATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCCAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGcc                                                                                           >  1:318111/1‑251 (MQ=255)                                                                                                                                                  caTTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAAt                                                                                     <  2:132749/251‑1 (MQ=255)                                                                                                                                                       cAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCgg                                                                                <  1:37191/251‑1 (MQ=255)                                                                                                                                                               gAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCt                                                                          <  1:210405/249‑1 (MQ=255)                                                                                                                                                                            cAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTc                                                                                                                                                               >  1:246925/1‑151 (MQ=255)                                                                                                                                                                            cAACGCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTGTGAACAACGTACGCGGCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCGTCGCAAGCTATCCGCGCCAGTc                                                                                                                                                               <  2:246925/151‑1 (MQ=255)                                                                                                                                                                                       tcttTTAGTCGGCAATTCAAACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGCACAACGTACGCGTCTGCAGATTCTTCCTGCGCAGTCAGTTATCGCCAGTGATGCCAGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAACATTATTACTc                                                 >  2:385088/1‑250 (MQ=255)                                                                                                                                                                                       cgtttaaGTCGGCAATTCCGACGCCATCCCGCCATTACTGCCTACAGCTGACTGTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCCGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGGTCAATGCGCGTGGCCCTGGGGCTGGTGAAGGAAGGTCGATCGCAATCCTGTGTCAGTGCCGGTAATACCGGGGTGCTGATGGGGCTGGCAAAATTATTACt                                                  <  2:413839/243‑1 (MQ=255)                                                                                                                                                                                          tttAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAg                                              <  1:148409/250‑1 (MQ=255)                                                                                                                                                                                                 ggCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGCACAACGTACGCGTCTGCAGCTTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGAAACATTATTACTCAAGCCCCTgg                                       >  1:585059/1‑250 (MQ=255)                                                                                                                                                                                                    aatcccgccGCCATCACGCCATTACTTGCTACAGCTTACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCGATCCGCGCCATTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAggt                                   <  2:513143/243‑2 (MQ=255)                                                                                                                                                                                                         cGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTg                              <  1:255952/251‑1 (MQ=255)                                                                                                                                                                                                                  cACGCCATTACTTGCTAGAGCTGACTTTGAACAACGTCCGCGTCTGCAGATTATTACTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCACTGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCACAATTATTACTCAAGCCCATGGAGGGGATTGAGCGTCCgg                     >  2:315328/1‑251 (MQ=255)                                                                                                                                                                                                                           aCTTGCTAAAGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATGGAGCGTCCGGCGCTGGTg             <  1:480051/250‑1 (MQ=255)                                                                                                                                                                                                                                   aaGCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTa    <  1:242782/251‑1 (MQ=255)                                                                                                                                                                                                                                     gCTGACTTTGAACAACGTACGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTAcc  <  1:476912/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                         CATCTGCGTGATGAAGCTTAGTGAGGATTTTCCCCAGGCAACTGGGGAAAGACCAAACCGGGCGGCGACGATACCTTGACACGTCTAACCCTGGCGTTAGATGTCATGGGAGGGGATTTTGGCCCTTCCGTGACAGTGCCTGCAGCATTGCAGGCACTGAATTCTAATTCGCAACTCACTCTTCTTTTAGTCGGCAATTCCGACGCCATCACGCCATTACTTGCTAAAGCTGACTTTGAACAACGTTCGCGTCTGCAGATTATTCCTGCGCAGTCAGTTATCGCCAGTGATGCCCGGCCTTCGCAAGCTATCCGCGCCAGTCGTGGGAGTTCAATGCGCGTGGCCCTGGAGCTGGTGAAAGAAGGTCGAGCGCAAGCCTGTGTCAGTGCCGGTAATACCGGGGCGCTGATGGGGCTGGCAAAATTATTACTCAAGCCCCTGGAGGGGATTGAGCGTCCGGCGCTGGTGACGGTATTACC  >  NC_000913/1147546‑1148024

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