Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,828,954	T→C	S4P (TCC→CCC)	xanP →	xanthine permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,828,954	0	T	C	96.5%	288.3 / 0.1	86	S4P (TCC→CCC)	xanP	xanthine permease
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (42/41);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 2.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

AAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCC  >  NC_000913/3828858‑3829050                                                                                                 |                                                                                                 aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCa                                                                                               >  1:227659/1‑100 (MQ=255) aaaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCa                                                                                               <  1:6203/100‑1 (MQ=255)  aaTATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCAc                                                                                              >  2:85378/1‑100 (MQ=255)     atTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCt                                                                                           >  2:188869/1‑100 (MQ=255)      tttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTc                                                                                          >  2:694292/1‑100 (MQ=255)       ttAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCg                                                                                         <  1:623099/100‑1 (MQ=255)        tAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGa                                                                                        >  2:1587319/1‑100 (MQ=255)        naGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGa                                                                                        >  2:623097/2‑100 (MQ=255)         aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAg                                                                                       <  1:588038/100‑1 (MQ=255)         aGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAg                                                                                       <  1:1905185/100‑1 (MQ=255)            tGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCa                                                                                    <  1:1563651/100‑1 (MQ=255)             gAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAg                                                                                   <  2:1538274/100‑1 (MQ=255)              aGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGa                                                                                  <  1:622151/100‑1 (MQ=255)               ggATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGan                                                                                 <  1:1498915/100‑2 (MQ=255)               ggATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGaa                                                                                 >  1:1398287/1‑100 (MQ=255)                  tgtgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAAt                                                                              <  2:1327356/100‑1 (MQ=255)                    tgGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgc                                                                            >  1:375777/1‑100 (MQ=255)                     ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcg                                                                           >  2:27515/1‑100 (MQ=255)                     ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcg                                                                           >  2:303440/1‑100 (MQ=255)                     ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcg                                                                           >  2:1882212/1‑100 (MQ=255)                     ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcg                                                                           <  2:1681804/100‑1 (MQ=255)                     ggCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcg                                                                           >  1:358446/1‑100 (MQ=255)                      gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcgc                                                                          <  1:1132969/100‑1 (MQ=255)                      gCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATgcgc                                                                          <  1:199617/100‑1 (MQ=255)                       cAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCa                                                                         <  1:1276215/100‑1 (MQ=255)                             tttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACAACGGGACTGTTATCGATGTCTGCTCCCACCCTCGAGTCAGAAAATGCGCAACCgg                                                                    <  2:414184/99‑1 (MQ=255)                              ttGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGtt                                                                  >  2:888531/1‑100 (MQ=255)                                gCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                >  1:1879864/1‑100 (MQ=255)                                gCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                >  1:940969/1‑100 (MQ=255)                                gCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgc                                                                >  2:560391/1‑100 (MQ=255)                                 cTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTgcg                                                               >  1:1943876/1‑100 (MQ=255)                                    ttCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCat                                                            <  1:888534/100‑2 (MQ=255)                                     tctcTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGa                                                           <  1:318194/100‑1 (MQ=255)                                      ctctTTCGGTCAGCTAAAACGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattt                                                          <  1:632919/100‑4 (MQ=255)                                       tctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACt                                                         <  1:1154154/100‑1 (MQ=255)                                        ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttc                                                        <  1:1580588/100‑6 (MQ=255)                                        ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttc                                                        <  2:336548/100‑6 (MQ=255)                                        ctTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTc                                                        >  1:545579/1‑100 (MQ=255)                                         tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttct                                                       <  1:1814602/100‑7 (MQ=255)                                         tttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCa                                                       <  1:357999/100‑1 (MQ=255)                                          ttCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttctg                                                      <  1:817697/100‑8 (MQ=255)                                            cGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCaaaa                                                    <  2:1860409/100‑1 (MQ=255)                                            cGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCattttctgac                                                    <  1:335195/100‑10 (MQ=39)                                             ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGGTTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAGATGCGCAACCGGTTGCGCAGACTCAAAAc                                                   >  2:1230205/1‑100 (MQ=255)                                             ggTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAAc                                                   <  1:1330545/100‑1 (MQ=255)                                              gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTGTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACa                                                  >  2:1133947/1‑100 (MQ=255)                                              gTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACa                                                  >  1:636334/1‑100 (MQ=255)                                                cAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGc                                                <  1:1291419/100‑1 (MQ=255)                                                 aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGCGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                               >  1:402598/1‑100 (MQ=255)                                                 aGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                               >  2:805065/1‑100 (MQ=255)                                                  gCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCg                                               >  1:220787/1‑99 (MQ=255)                                                    tAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAAc                                            <  1:169578/100‑1 (MQ=255)                                                        aTGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAt                                        <  1:1479635/100‑1 (MQ=255)                                                         tGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                       >  1:871705/1‑100 (MQ=255)                                                         tGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAtt                                       <  1:1882230/100‑1 (MQ=255)                                                          gCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGAttt                                      >  2:486883/1‑100 (MQ=255)                                                           cccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTa                                     >  2:63844/1‑100 (MQ=255)                                                            ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAc                                    >  2:173725/1‑100 (MQ=255)                                                            ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAc                                    >  2:543884/1‑100 (MQ=255)                                                            ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAc                                    >  1:1040074/1‑100 (MQ=255)                                                            ccGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTAc                                    <  1:1821909/100‑1 (MQ=255)                                                             cGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGTAAATGCGCAGCCGGTGGCGCAGACACAAAACAGCGAACTGATTTAcc                                   >  2:1377496/1‑100 (MQ=255)                                                              gCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCg                                  <  2:96147/100‑1 (MQ=255)                                                                 ttGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCt                               <  1:1988611/100‑1 (MQ=255)                                                                  tGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                              <  2:823955/100‑1 (MQ=255)                                                                  tGCTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCtt                              >  1:1113878/1‑100 (MQ=255)                                                                    cTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGa                            <  1:1133952/100‑1 (MQ=255)                                                                    cTACCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGAATCAAAACAGCGAACTGATTTACCGTCttta                            >  2:1702876/1‑98 (MQ=255)                                                                      aCCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAg                          <  1:188870/100‑1 (MQ=255)                                                                      aCCACGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAg                          <  2:227658/100‑1 (MQ=255)                                                                         aCGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATc                       <  1:1964931/100‑1 (MQ=255)                                                                          cGGGATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGATCTGATTTACCGTCTTGAAGATCg                      <  2:220786/100‑1 (MQ=255)                                                                            ggATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTc                    >  1:1710688/1‑100 (MQ=255)                                                                             gATTGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTcc                   <  1:684319/100‑1 (MQ=255)                                                                               ttGTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGACCTGATTTACCGTCTTGAAGATCGTccgc                 <  2:409798/100‑1 (MQ=255)                                                                                 gTTTTCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTccgccg               >  1:789736/1‑100 (MQ=255)                                                                                   tttCGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCt             >  1:76420/1‑100 (MQ=255)                                                                                      cGATGTCTGTTCCCACCCTCGAGTCAGAAAATGTGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc          >  2:574806/1‑100 (MQ=255)                                                                                      cGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAATCAGCGAACTGCTTTACCGTCTTGAAGATTGTCCGCCGCTTcc          <  1:1377505/100‑1 (MQ=255)                                                                                      cGATGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTcc          >  2:1042655/1‑100 (MQ=255)                                                                                         tGTCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCa       >  1:629071/1‑100 (MQ=255)                                                                                           tCTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCaaa     <  1:1230211/100‑1 (MQ=255)                                                                                            cTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc    <  2:402597/100‑1 (MQ=255)                                                                                            cTGTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAc    >  2:216457/1‑100 (MQ=255)                                                                                              gTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAccc  >  1:1041701/1‑100 (MQ=255)                                                                                              gTTCCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAAccc  >  1:1033905/1‑100 (MQ=255)                                                                                                 |                                                                                                 AAATATTTAGATGAGGATGTGGCAAACGTTTGCTTTTCTCTTTCGGTCAGCTAAAATGCCCGCTTTGCTACCACGGGATTGTTTTCGATGTCTGTTTCCACCCTCGAGTCAGAAAATGCGCAACCGGTTGCGCAGACTCAAAACAGCGAACTGATTTACCGTCTTGAAGATCGTCCGCCGCTTCCTCAAACCC  >  NC_000913/3828858‑3829050

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