Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,068,506	G→T	6.7%	A101A (GCC→GCA)	argO ←	arginine transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,068,506	0	G	T	6.7%	225.2 / 2.9	67	A101A (GCC→GCA)	argO	arginine transporter
Reads supporting (aligned to +/- strand):  ref base G (29/33);  new base T (3/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAAC  >  NC_000913/3068370‑3068629                                                                                                                                         |                                                                                                                            cgCGTTTTGGTTCCCCATCACGTTGCCCGCCAAGGCTGCCCAGTCCCACCAAAGTATCCAGGTCACGATCGGGTTTCCGCCAGGTCTCTGCCAACATGGTGGCGATACGTTTCCATCTGCCTTGCTTCCTGACTTCGgcg                                                                                                                          <  2:1594965/139‑1 (MQ=255) aGCGTTGTGGTTCCACATCAAGTTGCCCGCGAAGGCTGCCCAGTACAACAAAAGTATCCAGGTACCCATGCGGATTCAGCCAGGTCACTGCCAACAGGGTGGCGATAAATTTCCCTCTGCCTTGCTTCATGACTTCGgcg                                                                                                                          <  2:1594976/140‑1 (MQ=255)   cGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCt                                                                                                                        >  1:30335/1‑140 (MQ=255)      tttGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGc                                                                                                                     <  2:612003/140‑1 (MQ=255)      tttGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGc                                                                                                                     <  2:611926/140‑1 (MQ=255)      tttGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGc                                                                                                                     <  2:30335/140‑1 (MQ=255)      ttaggTTCCACTTCAAGGTGCCCGCCAAGGCTGCCCTGTACTACAAATGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGTTAATTTTCCTTCTGCCTTGCTTCATGACTTCGGCGCTGGc                                                                                                                     <  2:1801004/137‑1 (MQ=255)      ttaggTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTCCAACAAAAGTATCCTGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGc                                                                                                                     <  2:842335/137‑1 (MQ=255)       ttGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCt                                                                                                                    <  2:1034623/140‑1 (MQ=255)        tGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGTCTTCTGCGCTGGCTa                                                                                                                   >  2:1890149/1‑140 (MQ=255)        tGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCTGCGCTGGCTa                                                                                                                   >  2:836606/1‑140 (MQ=255)          gTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCGAAc                                                                                                                 >  2:1175828/1‑140 (MQ=255)           ttCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATTACTTCGGCGCTGGCTAACt                                                                                                                <  2:1593156/140‑1 (MQ=255)              cacaTCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCTGCGCTGGCTAACTCaa                                                                                                             <  1:836606/140‑1 (MQ=255)              cacaTCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCTGCGCTGGCTAACTCaa                                                                                                             <  1:1890149/140‑1 (MQ=255)              cacaTCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCaa                                                                                                             <  1:1190333/140‑1 (MQ=255)              cacaTCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCaa                                                                                                             <  2:359667/140‑1 (MQ=255)                     aGTTGCCCGCCAAGGCTGCCAAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGTCGCTGGCTAACTCAATATTaat                                                                                                      >  1:39070/1‑138 (MQ=255)                              ccAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGc                                                                                             <  1:763673/140‑1 (MQ=255)                               cAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCt                                                                                            <  2:1818300/140‑1 (MQ=255)                                       cccAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACTTGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACGACTGCGCTGGCTCACTCAATATTACTGCTCATTGCTGTTTGaaa                                                                                    >  2:30292/1‑140 (MQ=255)                                       cccAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTACATGACTTCGGCGCTGGCTAACTCAATAGTACTGCTCATTGCTGTTGTaaa                                                                                    >  2:30211/1‑140 (MQ=255)                                             acaacaAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGcc                                                                              <  1:1886989/140‑1 (MQ=255)                                                 caaAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCaaaa                                                                          >  1:1886791/1‑140 (MQ=255)                                                    aaGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCGGCCTTGCTTCTTGAGTTCGGCGCTGGCTAACTCAATCTTCCTGCTCATTGCTGTTTTATACGCGCCAGTACCa                                                                       >  2:1494682/1‑140 (MQ=255)                                                      gTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCCACATGGACGCGATAATTTTCCATCTCCCTTGCTTCATGACTTCGGCGCTGGATAACTCAATATTACTGCTCATTGCTGTGATAAAAGCGACAAAACCaaa                                                                     >  2:1581821/1‑138 (MQ=255)                                                      gTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATa                                                                     >  1:76644/1‑140 (MQ=255)                                                      gTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATa                                                                     >  1:902322/1‑140 (MQ=255)                                                      gTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATa                                                                     >  1:1913199/1‑140 (MQ=255)                                                      gTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATa                                                                     >  2:525326/1‑140 (MQ=255)                                                      gTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATa                                                                     >  2:525479/1‑140 (MQ=255)                                                              gTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTGTTAAAAGCGCCAAAACCATACCACAGCa                                                             >  2:353779/1‑140 (MQ=255)                                                               tAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCaa                                                            <  2:1373708/140‑1 (MQ=255)                                                                aaccATGCGGATTCAGCAAGGTCACTGCCAACATGGTGGCGATTATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCaag                                                           <  2:699760/137‑1 (MQ=255)                                                                aaaCATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCaag                                                           <  2:324364/140‑1 (MQ=255)                                                                 aaCATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTCACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCaaga                                                          >  2:1217656/1‑140 (MQ=255)                                                                 aaCATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCaaga                                                          >  2:1217673/1‑140 (MQ=255)                                                                       cGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTaa                                                                                                                  <  1:425052/78‑1 (MQ=255)                                                                       cGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTaa                                                                                                                  >  2:425052/1‑78 (MQ=255)                                                                       cGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTa                                                    >  2:222919/1‑140 (MQ=255)                                                                       cGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTa                                                    >  2:1698320/1‑140 (MQ=255)                                                                        ggATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAc                                                   <  1:1605756/140‑1 (MQ=255)                                                                        ggATTCACCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAc                                                   <  2:1759229/140‑1 (MQ=255)                                                                             cAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACCGCAAGAAGGCTAcgccgc                                              <  1:1665720/140‑1 (MQ=255)                                                                             cAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAcgccgc                                              <  2:315961/140‑1 (MQ=255)                                                                             cAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAcgccgc                                              <  2:1788170/140‑1 (MQ=255)                                                                             cAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAcgccgc                                              <  1:1522514/140‑1 (MQ=255)                                                                             cAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAcgccgc                                              <  1:1364182/140‑1 (MQ=255)                                                                             cAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTAcgccgc                                              <  1:1363244/140‑1 (MQ=255)                                                                              aGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACgccgcc                                             >  1:364546/1‑140 (MQ=255)                                                                                    gTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGt                                       >  1:1525781/1‑140 (MQ=255)                                                                                    gTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGt                                       >  2:892742/1‑140 (MQ=255)                                                                                         tGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCa                                  >  2:1096827/1‑140 (MQ=255)                                                                                              aCATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGcc                             >  2:103500/1‑140 (MQ=255)                                                                                                aTGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAg                           <  1:463358/140‑1 (MQ=255)                                                                                                aTGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAg                           <  1:768334/140‑1 (MQ=255)                                                                                                aTGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAg                           <  1:919843/140‑1 (MQ=255)                                                                                                   gtggCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCaa                        >  1:165325/1‑140 (MQ=255)                                                                                                     ggCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAAcc                      >  1:210819/1‑140 (MQ=255)                                                                                                        gATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACg                   <  2:512469/140‑1 (MQ=255)                                                                                                             ctttCCAACTGCCTTCCTTCCTGACTACGGACCTCGCTAACTCAAAATTACTTCTCATTGCTGTTTAAAAAGCGCCAAAACCAAACCACAGCAAGAAGGCTACGCCGCCACAGGTGACCAGCGCCAGCAACCACGGCGAc              <  2:1462616/139‑1 (MQ=255)                                                                                                                tCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGc           <  2:203912/140‑1 (MQ=255)                                                                                                                tCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGc           >  2:1041426/1‑140 (MQ=255)                                                                                                                tCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGc           >  2:804729/1‑140 (MQ=255)                                                                                                                  cATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCAt         <  1:1058222/140‑1 (MQ=255)                                                                                                                        ccTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATaa   >  1:1867447/1‑140 (MQ=255)                                                                                                                         cTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAAc  >  2:1441075/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                            AGCGTTTTGGTTCCACATCAAGTTGCCCGCCAAGGCTGCCCAGTACAACAAAAGTATCCAGGTAAACATGCGGATTCAGCCAGGTCACTGCCAACATGGTGGCGATAATTTTCCATCTGCCTTGCTTCATGACTTCGGCGCTGGCTAACTCAATATTACTGCTCATTGCTGTTTTAAAAGCGCCAAAACCATACCACAGCAAGAAGGCTACGCCGCCCCAGGTGACCAGCGCCAGCAACCACGGCGACTGCATCAATAAC  >  NC_000913/3068370‑3068629

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