Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,833,205	G→A	A260T (GCC→ACC)	xthA →	exonuclease III

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,833,205	0	G	A	98.2%	203.7 / ‑4.3	57	A260T (GCC→ACC)	xthA	exonuclease III
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (36/20);  total (37/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

TGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTTAG  >  NC_000913/1833061‑1833352                                                                                                                                                 |                                                                                                                                                    tGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCg                                                                                                                                               >  2:624106/1‑151 (MQ=255)    tttGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCt                                                                                                                                            >  2:840996/1‑151 (MQ=255)         ttACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCg                                                                                                                                       >  2:501941/1‑151 (MQ=255)           acagCTCAAAAGGTTTTGACGGTAAAAGAGGTCTGCGCCTCGACCTGCTGCTCGCGAGCCAACCGCTGGCAGAATAATGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGAc                                                                                                                                     >  2:219624/4‑151 (MQ=255)             cGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCt                                                                                                                                   >  2:617441/1‑151 (MQ=255)                tCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTcc                                                                                                                                >  2:291223/1‑151 (MQ=255)                 cAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTccg                                                                                                                               <  1:501995/151‑1 (MQ=255)                  aaaaGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTccgc                                                                                                                              >  1:537707/1‑151 (MQ=255)                  aaaaGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTccgc                                                                                                                              >  1:635128/1‑151 (MQ=255)                   aaaGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTccgcc                                                                                                                             >  1:839663/1‑151 (MQ=255)                    aaGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTccgccg                                                                                                                            <  2:635063/151‑1 (MQ=255)                     aGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCcgccgc                                                                                                                           >  1:445227/1‑151 (MQ=255)                     aGGGTTTGACGATAACCGTGGTCTGCGCATCGTCCTGCTTCTCGCAAGCTAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGGCTATGAAATCCGCAGCGTGGCAAAACCGTCCGATAACACCCCCCTCTGGGCGGCCTTCcgccgc                                                                                                                           >  2:706024/1‑151 (MQ=255)                           tGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAAttt                                                                                                                     >  1:513711/1‑151 (MQ=255)                           tGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAAttt                                                                                                                     >  2:230612/1‑151 (MQ=255)                            gACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAAttt                                                                                                                     <  1:916031/150‑1 (MQ=255)                                     cGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCtctc                                                                                                           >  2:972841/1‑151 (MQ=255)                                      gTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTcc                                                                                                          >  2:270020/1‑151 (MQ=255)                                              cgcATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCaa                                                                                                  <  2:91320/151‑1 (MQ=255)                                               gcATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCATCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTag                                                                                                                   <  2:354625/133‑1 (MQ=255)                                               gcATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCATCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTag                                                                                                                   >  1:354668/1‑133 (MQ=255)                                                  tCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTg                                                                                              >  2:887729/1‑151 (MQ=255)                                                        tgctgcTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTc                                                                                          >  1:380978/1‑149 (MQ=255)                                                           tgctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAgg                                                                                       >  1:398223/1‑149 (MQ=255)                                                            gctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGgaga                                                                                    <  1:463235/151‑1 (MQ=255)                                                            gctCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGgag                                                                                     <  2:815521/150‑1 (MQ=255)                                                                gccagccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAAtt                                                                                 >  1:275639/1‑150 (MQ=255)                                                                 ccagccaACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTaa                                                                               <  1:292010/151‑1 (MQ=255)                                                                      caACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGTCTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                          <  2:627877/151‑1 (MQ=255)                                                                       aaCCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                          >  1:303942/1‑150 (MQ=255)                                                                        aCCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTgag                                                                        <  1:980276/151‑1 (MQ=255)                                                                        aCCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTg                                                                          >  1:816773/1‑149 (MQ=255)                                                                                cAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATc                                                                 >  1:387323/1‑150 (MQ=255)                                                                                 aGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATcaa                                                               <  2:716665/151‑1 (MQ=255)                                                                                 aGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATca                                                                <  2:448733/150‑1 (MQ=255)                                                                                  gAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATcaac                                                              >  2:540663/1‑151 (MQ=255)                                                                                    aTGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATcaaca                                                             >  2:669774/1‑150 (MQ=255)                                                                                      gTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAAc                                                          <  1:500642/151‑1 (MQ=255)                                                                                        tGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACt                                                         >  1:223914/1‑150 (MQ=255)                                                                                             aGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAg                                                    <  1:436566/150‑1 (MQ=255)                                                                                             aGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAg                                                    <  2:18128/150‑1 (MQ=255)                                                                                              gAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGt                                                   >  1:582425/1‑150 (MQ=255)                                                                                               aaaCCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTaa                                                 <  1:94522/151‑1 (MQ=255)                                                                                                  ccGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATg                                               >  1:962672/1‑150 (MQ=255)                                                                                                    ggCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGAtt                                            >  2:443527/1‑151 (MQ=255)                                                                                                              aTGAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTtgcctgc                                  >  2:437755/1‑151 (MQ=255)                                                                                                                gAAATCCGCAGCATGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCg                                >  2:437248/1‑151 (MQ=255)                                                                                                                       gcagcaTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTg                         <  2:839585/151‑1 (MQ=255)                                                                                                                        cagcaTGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGt                        <  1:299803/151‑1 (MQ=255)                                                                                                                             tGGAAAAACCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATAc                   >  2:625582/1‑151 (MQ=255)                                                                                                                                 aaaaCCGTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGtc               <  2:831786/151‑1 (MQ=255)                                                                                                                                       gTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCg         <  1:270056/151‑1 (MQ=255)                                                                                                                                       gTCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCg         >  1:573228/1‑151 (MQ=255)                                                                                                                                        tCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGt        >  1:516542/1‑151 (MQ=255)                                                                                                                                        tCCGATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGt        >  2:133658/1‑151 (MQ=255)                                                                                                                                           gATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGtttt     >  2:77881/1‑151 (MQ=255)                                                                                                                                           gATCACACCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGtttt     >  2:795754/1‑151 (MQ=255)                                                                                                                                              cacaCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGttttt    >  2:567029/1‑149 (MQ=255)                                                                                                                                              cacaCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGttttt    >  1:28513/1‑149 (MQ=255)                                                                                                                                              cacaCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTTAg  <  2:704671/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    TGGTTTGATTACCGCTCAAAAGGTTTTGACGATAACCGTGGTCTGCGCATCGACCTGCTGCTCGCCAGCCAACCGCTGGCAGAATGTTGCGTAGAAACCGGCATCGACTATGAAATCCGCAGCATGGAAAAACCGTCCGATCACGCCCCCGTCTGGGCGACCTTCCGCCGCTAATTTAGCAGCTCTCCTGGCTCAAACTGGGTCAGGAGAATTAACCTTGAGAAAAATCAACAAACTGTCAGTAATGATTTGTTGCCTGCCGTCCTTTGTTATACCGTCTCTGCGTTTTTAG  >  NC_000913/1833061‑1833352

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