Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,566,544	C→A	S13I (AGC→ATC)	glgP ←	glycogen phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,566,544	0	C	A	98.9%	325.5 / ‑3.4	93	S13I (AGC→ATC)	glgP	glycogen phosphorylase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (51/41);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

GTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGC  >  NC_000913/3566448‑3566638                                                                                                  |                                                                                               gTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATaa                                                                                             >  1:1036167/1‑100 (MQ=255) gTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATaa                                                                                             >  1:2168278/1‑100 (MQ=255) gTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATaa                                                                                             <  2:337984/100‑1 (MQ=255)  ttGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAg                                                                                            <  2:1550081/100‑1 (MQ=255)     cGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCg                                                                                          <  1:1628851/99‑1 (MQ=255)      gTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTg                                                                                        >  1:242894/1‑100 (MQ=255)       ttCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTgg                                                                                       <  2:303935/100‑1 (MQ=255)       ttCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTgg                                                                                       >  2:479863/1‑100 (MQ=255)        tCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTggg                                                                                      <  2:150040/100‑1 (MQ=255)        tCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTggg                                                                                      <  1:2090589/100‑1 (MQ=255)        tCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTggg                                                                                      <  2:340888/100‑1 (MQ=255)         cAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGc                                                                                     >  2:1405026/1‑100 (MQ=255)         cAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGc                                                                                     >  1:1851739/1‑100 (MQ=255)           gCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGa                                                                                   >  2:2058829/1‑100 (MQ=255)             cATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATg                                                                                 <  2:37988/100‑1 (MQ=255)             cATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATg                                                                                 >  2:1250650/1‑100 (MQ=255)              aTTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGa                                                                                >  2:19692/1‑100 (MQ=255)                 cATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAata                                                                             >  2:981109/1‑100 (MQ=255)                   tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                           <  2:1179408/100‑1 (MQ=255)                   tGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATg                                                                           <  2:702436/100‑1 (MQ=255)                    gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                          >  1:920268/1‑100 (MQ=255)                    gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                          >  1:919849/1‑100 (MQ=255)                    gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                          >  2:1806037/1‑100 (MQ=255)                    gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                          >  2:404811/1‑100 (MQ=255)                    gTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGt                                                                          >  1:1866926/1‑100 (MQ=255)                      ttATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTATACGATAAGCGTGGGCGATGAATATGTaa                                                                        >  2:91775/1‑100 (MQ=255)                         ttGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACg                                                                     >  1:191985/1‑100 (MQ=255)                           ggCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGa                                                                   >  1:1708879/1‑100 (MQ=255)                           ggCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGa                                                                   >  1:1709098/1‑100 (MQ=255)                            gCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAg                                                                  <  1:1922941/100‑1 (MQ=255)                             cgacgaCCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGc                                                                 >  1:1111199/1‑100 (MQ=255)                               acgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAt                                                               <  1:1811402/100‑1 (MQ=255)                                cgacCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCAtt                                                              >  2:108456/1‑100 (MQ=255)                                  acCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCa                                                            <  1:2058844/100‑1 (MQ=255)                                  acCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATACAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGc                                                                 >  2:726294/1‑95 (MQ=255)                                    cGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCata                                                          <  1:284207/100‑1 (MQ=255)                                    cGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCata                                                          <  1:1876109/100‑1 (MQ=255)                                       gTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAg                                                       >  2:1602788/1‑100 (MQ=255)                                       gTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAg                                                       <  1:1466815/100‑1 (MQ=255)                                        tCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAgg                                                      <  1:356818/100‑1 (MQ=255)                                          cTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCg                                                    >  2:1546947/1‑100 (MQ=255)                                              ccAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTc                                                >  1:108197/1‑100 (MQ=255)                                               cAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                               <  2:1851727/100‑1 (MQ=255)                                               cAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTcc                                               >  2:1557965/1‑100 (MQ=255)                                                  tCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGa                                            >  1:759492/1‑100 (MQ=255)                                                    gTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaa                                          <  2:533080/100‑1 (MQ=255)                                                    gTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaaa                                          <  1:1900417/100‑1 (MQ=255)                                                      aaaCATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                                        <  1:122015/100‑1 (MQ=255)                                                      aaaCATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                                        >  2:1048950/1‑100 (MQ=255)                                                      aaaCATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAc                                        >  1:699133/1‑100 (MQ=255)                                                      aaaCATCAG‑TTGTAAAGGGAAAGATTGCTTAAGAGCTTCTACAATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGaa                                           <  2:1111192/97‑1 (MQ=255)                                                       aaCATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACt                                       >  1:1767290/1‑100 (MQ=255)                                                        aCATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTa                                      >  2:792810/1‑100 (MQ=255)                                                         cATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                                     >  1:127868/1‑100 (MQ=255)                                                         cATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAt                                     >  2:203643/1‑100 (MQ=255)                                                           tCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTAttt                                   <  1:1944890/100‑1 (MQ=255)                                                            cAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                                  >  1:1537908/1‑100 (MQ=255)                                                            cAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTc                                  <  2:1146706/100‑1 (MQ=255)                                                              gCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                                >  1:1400619/1‑100 (MQ=255)                                                              gCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCaa                                >  1:2133926/1‑100 (MQ=255)                                                                ttGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc                              >  2:1773383/1‑100 (MQ=255)                                                                ttGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc                              <  1:521563/100‑1 (MQ=255)                                                                ttGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGc                              <  1:232231/100‑1 (MQ=255)                                                                 tGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCg                             >  2:647365/1‑100 (MQ=255)                                                                   t‑AAGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGAt                           <  2:1161028/100‑1 (MQ=255)                                                                     aaGCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATa                          <  2:191982/100‑1 (MQ=255)                                                                       gCGATAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGt                        <  1:203646/100‑1 (MQ=255)                                                                           tAGAGTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAg                    >  1:859465/1‑100 (MQ=255)                                                                             gagTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCt                  >  1:924251/1‑100 (MQ=255)                                                                              agTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTc                 >  1:1331629/1‑100 (MQ=255)                                                                               gTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                >  2:89769/1‑100 (MQ=255)                                                                               gTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                >  2:52706/1‑100 (MQ=255)                                                                               gTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                <  2:310550/100‑1 (MQ=255)                                                                               gTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCa                >  2:1165363/1‑100 (MQ=255)                                                                               gTGCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAAATATTTAAAGCGATAGTAAAGCTCa                >  1:979170/1‑100 (MQ=255)                                                                                 gCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg              >  2:2052467/1‑100 (MQ=255)                                                                                 gCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg              <  1:255572/100‑1 (MQ=255)                                                                                 gCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg              <  2:242891/100‑1 (MQ=255)                                                                                 gCTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACg              >  1:247262/1‑100 (MQ=255)                                                                                  cTTAAGAGCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACgg             <  1:108459/100‑1 (MQ=255)                                                                                      tccgatcttACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg         <  2:857091/92‑1 (MQ=255)                                                                                      agagCTTCTACGATAAGCGTGGGCGATGAATATGTAACCGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg         <  1:91777/100‑1 (MQ=255)                                                                                      agagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg         <  1:1250657/100‑1 (MQ=255)                                                                                      agagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg         >  2:1555255/1‑100 (MQ=255)                                                                                      agagCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACg         <  2:1331620/100‑1 (MQ=255)                                                                                        agCTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGAc       <  2:759485/100‑1 (MQ=255)                                                                                          cTTCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACtt     >  2:1152684/1‑100 (MQ=255)                                                                                           ttCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTc    <  2:919842/100‑1 (MQ=255)                                                                                           ttCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTc    <  2:920261/100‑1 (MQ=255)                                                                                            tCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcg   >  1:556752/1‑100 (MQ=255)                                                                                            tCTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTc    <  2:2168263/99‑1 (MQ=255)                                                                                             cTACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTcgc  >  1:1753861/1‑100 (MQ=255)                                                                                              tACGATAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGagatcgg  >  1:857098/1‑93 (MQ=255)                                                                                                  |                                                                                               GTTGCGTTCAGCCATTCATGTTTATTGGCGACGACCGGGTCCTTTCCAATCGTAAACATCAGCTTGT‑AAGCGATAGAGTGCTTAAGAGCTTCTACGCTAAGCGTGGGCGATGAATATGTAAACGGAGCATTCATATAGGCGTTTCCTGAAAACTATTTCAAGCGATAGTAAAGCTCACGGTACGACTTCGC  >  NC_000913/3566448‑3566638

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