Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,835,622	G→A	32.8%	R36R (CGG→CGA)	ynjA →	carboxymuconolactone decarboxylase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,835,622	0	G	A	32.8%	107.5 / 62.9	82	R36R (CGG→CGA)	ynjA	carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base A (10/17);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

ACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGC  >  NC_000913/1835523‑1835721                                                                                                    |                                                                                                    aCCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCgg                                                                                                     >  2:257178/1‑100 (MQ=255)  ccgccgCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGa                                                                                                    <  1:612092/100‑1 (MQ=255)  ccgccgCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGa                                                                                                    >  1:1882703/1‑100 (MQ=255)          tAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGttt                                                                                            <  2:508113/100‑1 (MQ=255)          tAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGttt                                                                                            <  2:2052101/100‑1 (MQ=255)           aGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGttta                                                                                           <  2:1186362/100‑1 (MQ=255)            gCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTATGGCGAGGTGATAAGCACCATTCGCAGGTGGGGGCGGATCCCGtttat                                                                                          >  2:1587867/1‑100 (MQ=255)               aaaTTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGtttatctt                                                                                       <  1:763402/100‑1 (MQ=255)                aaTTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttt                                                                                      <  2:936446/100‑1 (MQ=255)                   tCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttat                                                                                   <  2:1822612/100‑1 (MQ=255)                     cTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttatct                                                                                 >  2:456938/1‑100 (MQ=255)                     cTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGtttatcttttatct                                                                                 <  2:2181875/100‑1 (MQ=255)                     cTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGtttatcttttatct                                                                                 >  2:784048/1‑100 (MQ=255)                          aTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGgtgt                                                                            <  1:1066039/100‑1 (MQ=255)                          aTTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGgtgt                                                                            >  1:1830957/1‑100 (MQ=255)                           ttttACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTc                                                                           <  2:2075513/100‑1 (MQ=255)                            tttACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCg                                                                          <  2:2333900/100‑1 (MQ=255)                              tACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAt                                                                        <  1:1090657/100‑1 (MQ=255)                              tACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAt                                                                        <  1:33204/100‑1 (MQ=255)                                cGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgt                                                                      >  2:357837/1‑100 (MQ=255)                                cGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGAtgt                                                                      <  2:15697/100‑1 (MQ=255)                                 gTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtt                                                                     <  1:578482/100‑1 (MQ=255)                                 gTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtt                                                                     >  2:771878/1‑100 (MQ=255)                                   ccACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttg                                                                   >  2:691765/1‑100 (MQ=255)                                   ccACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGAtgtttg                                                                   <  2:1875251/100‑1 (MQ=255)                                    cacaGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGAtgtttgt                                                                  <  2:2303935/100‑1 (MQ=255)                                      caGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTttgttg                                                                <  1:129701/100‑1 (MQ=255)                                      caGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTttgttg                                                                <  1:2056211/100‑1 (MQ=255)                                      caGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCCGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTttgttg                                                                >  1:1882613/1‑100 (MQ=255)                                       aGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtgg                                                               >  1:9193/1‑100 (MQ=255)                                        ggCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtggc                                                              >  1:2135848/1‑100 (MQ=255)                                          cGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCTCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtggctg                                                            <  2:2221229/100‑1 (MQ=255)                                          cGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCTATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtggctg                                                            <  2:1173741/100‑1 (MQ=255)                                          cGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTtggctg                                                            >  1:1860742/1‑100 (MQ=255)                                             ggCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTggctggct                                                         >  1:2044662/1‑100 (MQ=255)                                             ggCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTggctggct                                                         >  2:280272/1‑100 (MQ=255)                                              gCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGgctggctg                                                        <  2:416979/100‑1 (MQ=255)                                                  cATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAgc                                                    >  2:2273080/1‑100 (MQ=255)                                                  cATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAgc                                                    <  1:2110963/100‑1 (MQ=255)                                                   aTCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAgcg                                                   <  2:2341302/100‑1 (MQ=255)                                                      gtcgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCaa                                                >  2:62678/1‑100 (MQ=255)                                                        cgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAAc                                              <  1:940525/100‑1 (MQ=255)                                                        cgCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAAc                                              <  2:1564940/100‑1 (MQ=255)                                                        cgCCATTACGGCGAGGTGCGAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCGGGCTGGAGCGCAAAc                                              <  1:1021895/100‑1 (MQ=255)                                                          ccATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGc                                            >  1:385433/1‑100 (MQ=255)                                                           cATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCt                                           >  1:2036691/1‑100 (MQ=255)                                                            aTTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTc                                          <  2:1958369/100‑1 (MQ=255)                                                            aTTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTc                                          <  1:997405/100‑1 (MQ=255)                                                              tACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCAc                                        <  1:2017275/100‑1 (MQ=255)                                                              tACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCAc                                        >  1:1263402/1‑100 (MQ=255)                                                               aCGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCAcc                                       >  2:309675/1‑100 (MQ=255)                                                                  gCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCt                                    >  1:225345/1‑100 (MQ=255)                                                                   cGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTc                                   <  1:740492/100‑1 (MQ=255)                                                                    gAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCg                                  <  1:485593/100‑1 (MQ=255)                                                                      ggTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGAt                                <  2:1697485/100‑1 (MQ=255)                                                                        tGCTAAGCCCCATTCGCTGGTGGGGTCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATc                               <  2:2039901/99‑1 (MQ=255)                                                                        tGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATcc                              <  1:1986559/100‑1 (MQ=255)                                                                        tGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATcc                              >  2:1507903/1‑100 (MQ=255)                                                                        tGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATcc                              <  1:784051/100‑1 (MQ=255)                                                                        tGCTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTCGAGCGCAAACGCTCACCGCTCGATcc                              <  1:1433591/100‑1 (MQ=255)                                                                         gCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCg                             >  2:283195/1‑100 (MQ=255)                                                                         gCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCg                             >  1:1725750/1‑100 (MQ=255)                                                                         gCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCg                             >  1:810967/1‑100 (MQ=255)                                                                         gCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCg                             >  2:1476776/1‑100 (MQ=255)                                                                         gCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCg                             >  2:2221357/1‑100 (MQ=255)                                                                          cTAAGCCCCATTCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCgg                            <  1:1848919/100‑1 (MQ=255)                                                                           tAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGgt                           <  1:2221376/100‑1 (MQ=255)                                                                            aaGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGgta                          >  1:1830196/1‑100 (MQ=255)                                                                            aaGCCCCATTCGCCGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTGTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGgt                           <  2:185731/99‑1 (MQ=255)                                                                              gCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGgtagt                        >  1:2311746/1‑100 (MQ=255)                                                                               ccccATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGgtagta                       <  1:384736/100‑1 (MQ=255)                                                                                 ccATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACg                     >  1:418513/1‑100 (MQ=255)                                                                                 ccATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACg                     >  2:1109638/1‑100 (MQ=255)                                                                                     tCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCg                 <  2:2214478/100‑1 (MQ=255)                                                                                     tCGCTGGTGGGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCg                 >  1:932659/1‑100 (MQ=255)                                                                                          ggtggGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGt            >  1:1919016/1‑100 (MQ=255)                                                                                           gtggGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTc           <  2:519636/100‑1 (MQ=255)                                                                                            tggGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCa          <  1:1867879/100‑1 (MQ=255)                                                                                            tggGGGCGAATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCa          >  2:481124/1‑100 (MQ=255)                                                                                             gggggCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAg         <  1:2104206/100‑1 (MQ=255)                                                                                             gggggCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAg         >  2:1641360/1‑100 (MQ=255)                                                                                               gggCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcg       >  2:1886764/1‑100 (MQ=255)                                                                                               gggAGGAGCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACTATCGCTTGTCAgcg       <  2:425597/100‑1 (MQ=255)                                                                                                ggcggcTCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCCCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcgc      >  1:2050222/1‑100 (MQ=255)                                                                                                    gATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcgcgcgc  >  2:1896064/1‑100 (MQ=255)                                                                                                    gATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAgcgcgcgc  <  1:1520/100‑1 (MQ=255)                                                                                                    |                                                                                                    ACCGCCGCTTAGCAAAATTCCTTTAATTTTACGTCCACAGGCGTGGCTGCATCGTCGCCATTACGGCGAGGTGCTAAGCCCCATTCGCTGGTGGGGGCGGATCCCGTTTATCTTTTATCTGGTGTCGATGTTTGTTGGCTGGCTGGAGCGCAAACGCTCACCGCTCGATCCGGTAGTACGATCGCTTGTCAGCGCGCGC  >  NC_000913/1835523‑1835721

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: ‑