Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	67,900	C→G	V146V (GTG→GTC)	araA ←	L‑arabinose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	67,900	0	C	G	97.8%	144.6 / ‑3.2	45	V146V (GTG→GTC)	araA	L‑arabinose isomerase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (24/20);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGAT  >  NC_000913/67673‑68126                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCTTAcg                                                                                                                                                                                                             >  2:295283/1‑251 (MQ=255) cTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATAcg                                                                                                                                                                                                             >  1:293938/1‑251 (MQ=255)          ttCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAAt                                                                                                                                                                                                    >  2:335058/1‑251 (MQ=255)              ccaccTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGa                                                                                                                                                                                                  >  1:98392/1‑249 (MQ=255)                                        cAGGTATTGACGGAAAAACCGAACTTAATCTGTGCGACAATTTTACCCCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGt                                                                                                                                                                      <  1:228358/251‑1 (MQ=255)                                           gTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTg                                                                                                                                                                   >  2:92483/1‑251 (MQ=255)                                               tGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttc                                                                                                                                                               >  1:196914/1‑251 (MQ=255)                                                gagggaaaaaCCGACTTTGAGCTGGGCGGCAACTTTACCACCACCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttca                                                                                                                                                              <  1:270628/244‑1 (MQ=255)                                                    gagaAACCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGCGCCCGACCCTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGAAACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcaggtt                                                                                                                                                          >  1:367616/1‑251 (MQ=255)                                                     agaAACCGAACTTGATCTGTACGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTggttcaggttc                                                                                                                                                         >  2:366733/1‑251 (MQ=255)                                                       aaaCCGAACTTGATCTGTGCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCAt                                                                                                                                                       >  1:210499/1‑251 (MQ=255)                                                                 tGATCTGGGCGGCAACTTTATCCCCACGGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCat                                                                                                                                              <  1:3728/250‑1 (MQ=255)                                                                  gATCTGTGCGGCAACTTTACCCCCAGCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCata                                                                                                                                             <  1:130338/250‑1 (MQ=255)                                                                       gtgCGGCATCTTAACCGCCACCGGTGACCGCCTCTTCACGCATGTTATCCCCAAAGCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCg                                                                                                                                          <  1:14564/248‑1 (MQ=255)                                                                        tgCGGCAACTTTATCACCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGATATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATAc                                                                                                                                      >  2:185572/1‑251 (MQ=255)                                                                                             ggTGACCGCCACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGt                                                                                                                 >  1:126962/1‑251 (MQ=255)                                                                                                gcccGCCCCTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCAGGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGGCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGa                                                                                                              <  2:161347/249‑1 (MQ=255)                                                                                                 aCCGCCACTTCAAGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGCAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATACAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGa                                                                                                              >  2:175311/1‑250 (MQ=255)                                                                                                      cACTTCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTggg                                                                                                        <  1:335897/251‑1 (MQ=255)                                                                                                         ttCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGgtgt                                                                                                     >  1:156802/1‑251 (MQ=255)                                                                                                          tCACGCATGTTATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGtgtg                                                                                                    >  1:269250/1‑251 (MQ=255)                                                                                                                    tATCGCCAAAACGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTgcagc                                                                                          >  1:372308/1‑251 (MQ=255)                                                                                                                     agcGCCAAAACTGCAGACTTTCATATGACGGGTATCCTGGTTAGAAACCGCCGGACGCATCCAGGGGCCGATACGCTCAGGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGTCCGCCATGTGCGGTCTGGTTCAGGTTCATAAAGTCCAGATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTTTAATTgcagc                                                                                          <  2:306699/248‑1 (MQ=255)                                                                                                                     aTCGCCAAAACGGCAGACTTTCAGATGACGG‑TATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGGCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTgcagca                                                                                         <  1:242256/250‑1 (MQ=255)                                                                                                                       cGCCAAAACGGCAGTCTGTCAGATGACGGGGATCGTGTTTAGAAACCGCCGGTCGCATCAAGGAGCCGATACGCTCATGGGCTTTTTCATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGACCTGGGTGTGGAATTGCAGCAAc                                                                                       <  2:110064/251‑1 (MQ=255)                                                                                                                        gCCAAAACGGCAGACTTTCAGATGCCGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTGGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACg                                                                                      <  1:66590/251‑1 (MQ=255)                                                                                                                              actgCAGCCTTTCAGCAGACGGGTATCCTGTTTAGAACCCGCCTGACCCATCACGGCGCAGATAAGCTCATGGGCTTGTTCATCCTGCCAGTGACCGGTACCGCCGGCATGTTGCTGCCGCATACGCGCGCCACTGAAGCCGACCTCGCCACCGCCCTGTGCAGTCTGGTTCAGGCTCATAAAGTCCATATCGATACTGTCCCACGTCAGCGCCGCGTTGGCCTGGGTGTGGACTTGCAGCAATGGTttgt                                                                                >  2:363726/4‑251 (MQ=255)                                                                                                                                   aGACTTTCAGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGCTACGCTCATGGGCTTGTTTATCCTGCCAGTGCCCGGTAACGACGGCATGTTGCTGACGCATCCGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCCTATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTTCTGCAACGGTTTGTTGAGc                                                                           >  2:223024/1‑251 (MQ=255)                                                                                                                                           aGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTTAGCATGGTCAg                                                                   <  1:83978/251‑1 (MQ=255)                                                                                                                                           aGATGACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAg                                                                   <  1:188849/251‑1 (MQ=255)                                                                                                                                              ggACGGGTATCCTGTTTAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGcc                                                                <  2:123889/250‑1 (MQ=255)                                                                                                                                                             ttAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACAtttt                                                 >  2:288013/1‑251 (MQ=255)                                                                                                                                                              tAGAAACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTg                                                >  1:241195/1‑251 (MQ=255)                                                                                                                                                               aaaaaCCGCCTGCCGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTgg                                               <  2:224380/247‑1 (MQ=255)                                                                                                                                                                        cTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCgggg                                         >  1:307955/1‑248 (MQ=255)                                                                                                                                                                                 ccAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGc                              <  1:93665/250‑1 (MQ=255)                                                                                                                                                                                      aGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCcac                        >  1:258264/1‑251 (MQ=255)                                                                                                                                                                                        ccGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCacac                      >  2:69971/1‑251 (MQ=255)                                                                                                                                                                                             aCGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACaccacc                  <  1:316319/250‑1 (MQ=255)                                                                                                                                                                                              cGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCaccag                <  1:81523/251‑1 (MQ=255)                                                                                                                                                                                                  aaTGGGCTTGTTGACCCGGCCCGTGACCGGGAACGACGGCAGTTTGCGGACGAATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGGGCAGTCTGGTTCAGGTTCATAAATTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACAccaccacaaca            <  2:195063/250‑6 (MQ=255)                                                                                                                                                                                                         ttGTTTATCCTGCCAGAGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGTCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAAc     >  2:324751/1‑251 (MQ=255)                                                                                                                                                                                                           gTTTATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGCTCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAgcgccaca    >  2:206595/1‑246 (MQ=255)                                                                                                                                                                                                           gTTTACCCTGCAAGTGACCGTTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCAGAGCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTGGATCCACATTTTGGCAGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACg    <  2:300153/250‑1 (MQ=255)                                                                                                                                                                                                             ttATCCTGCCAGTGACCGGTAACGACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGACCTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAGCTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGGGCAGCCACACCACCAGACCAGCGCAACGAt  <  1:248110/250‑1 (MQ=255)                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   CTGATGGAGTTCACCACCTGCACCAGATCGCCAACCGCCCAGGTATTGACGGAGAAACCGAACTTGATCTGTGCGGCAACTTTATCGCCATCGGTGACCGCCACTTCACGCATGTTATCGCCAAATCGGCAGACTTTCAGATGACGGGTATCCTGTTTAGAGACCGCCTGACGCATCCAGGAGCCGATACGCTCATGGGCTTGTTTATCCTGCCAGTGACCGGTAACCACGGCATGTTGCTGACGCATACGCGCGCCAATGAAGCCGAACTCGCGACCGCCATGTGCAGTCTGGTTCAGGTTCATAAAGTCCATATCGATACTGTCCCACGGCAGCGCCGCGTTGAACTGGGTGTGGAATTGCAGCAACGGTTTGTTGAGCATGGTCAGGCCGTTGATCCACATTTTGGCCGGGGAGAAGGTGTGCAGCCACACCACCAGACCAGCGCAACGAT  >  NC_000913/67673‑68126

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