Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,946,790	T→G	100%	G180G (GGT→GGG)	yebB →	DUF830 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,946,790	0	T	G	94.3%	285.7 / 3.7	86	G180G (GGT→GGG)	yebB	DUF830 family protein
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (35/46);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCG  >  NC_000913/1946698‑1946887                                                                                             |                                                                                                  gCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGCGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgtgt                                                                                            >  2:642696/1‑100 (MQ=255) gCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGgtgt                                                                                            >  1:199930/1‑100 (MQ=255)  cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCCTCAGGGGTTGGTGtt                                                                                           >  1:852859/1‑100 (MQ=255)  cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGtt                                                                                           <  2:1047892/100‑1 (MQ=255)     aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAt                                                                                        <  1:67762/100‑1 (MQ=255)      aCTCGCTTTCTGGAAATTCTGGATCCTGGGTTCTATTCCGTGGGAGCGTTAAACCGTCTCGCCAGCCCGGTTGTGGCATCATCCGGGGATGGTGTTGAt                                                                                        >  2:2293751/1‑99 (MQ=255)      aCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGAt                                                                                        <  2:975092/99‑1 (MQ=255)       cTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTc                                                                                      <  2:2263479/100‑1 (MQ=255)        tCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCa                                                                                     <  1:498191/100‑1 (MQ=255)        tCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCa                                                                                     <  2:1735979/100‑1 (MQ=255)          aCTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATg                                                                                   <  1:32491/100‑1 (MQ=255)           cTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGc                                                                                  <  2:189750/100‑1 (MQ=255)           cTTTCTGGAAATTCTGGTACTTAGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGc                                                                                  >  1:2308578/1‑100 (MQ=255)            tgtCTGGAAATTCTGGTTCCTGGGTTCGATTTCGGGGGATCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATg                                                                                   <  2:1468296/96‑1 (MQ=255)             ttCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCgg                                                                                <  1:2232130/100‑1 (MQ=255)              tCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCgga                                                                               >  2:2232116/1‑100 (MQ=255)               cTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggag                                                                              <  2:2452562/100‑1 (MQ=255)                tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggagg                                                                             >  2:583895/1‑100 (MQ=255)                 ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggaggg                                                                            >  1:1865885/1‑100 (MQ=255)                   aaaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggagggag                                                                          >  1:56006/1‑100 (MQ=255)                    aaTTCTGGTTCCTGTGTTCTATTCCGTGTGAGCGTAAAACCTTCACGCCATCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGt                                                                         <  1:1492554/100‑1 (MQ=255)                    aaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGt                                                                         <  1:285542/100‑1 (MQ=255)                    aaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGt                                                                         <  2:983687/100‑1 (MQ=255)                     aTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTg                                                                        <  2:1436966/100‑1 (MQ=255)                       tCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGa                                                                      >  2:652109/1‑100 (MQ=255)                             tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                                <  2:402844/100‑1 (MQ=255)                             tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                                <  2:650982/100‑1 (MQ=255)                               cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                              >  2:366019/1‑100 (MQ=255)                               cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                              <  2:2071385/100‑1 (MQ=255)                               cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                              >  1:485440/1‑100 (MQ=255)                                tGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCa                                                             <  2:1325470/100‑1 (MQ=255)                                tGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGAATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCa                                                             <  2:852851/100‑1 (MQ=255)                                 gggTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAg                                                            >  2:824790/1‑100 (MQ=255)                                  ggTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGc                                                           <  2:643206/100‑1 (MQ=255)                                   gTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGcc                                                          >  1:1692473/1‑100 (MQ=255)                                    ttCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCg                                                         <  2:1806204/100‑1 (MQ=255)                                      cTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGa                                                       >  1:455440/1‑100 (MQ=255)                                      cTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGa                                                       <  2:92082/100‑1 (MQ=255)                                      cTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGa                                                       <  1:699677/100‑1 (MQ=255)                                      cTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGa                                                       <  2:872546/100‑1 (MQ=255)                                        aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGaan                                                     <  1:2020237/100‑2 (MQ=255)                                        aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                     >  1:1987322/1‑100 (MQ=255)                                        aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                     <  1:2495005/100‑1 (MQ=255)                                         ttCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACt                                                    <  2:349859/100‑1 (MQ=255)                                             gTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGGGTTGATTCATGTGGAGGGAGTGGAAACGCCTCAGCCGGAACTGAcc                                                <  2:431341/100‑1 (MQ=255)                                              tGGGAGCGTAAAACCTTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCg                                               <  2:600554/100‑1 (MQ=255)                                              tGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCg                                               <  2:605360/100‑1 (MQ=255)                                                 gAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAgg                                            >  2:2210850/1‑100 (MQ=255)                                                 gAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAgg                                            >  2:57578/1‑100 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                         >  2:990253/1‑100 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                         <  1:1219211/100‑1 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                         >  1:1665851/1‑100 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                         <  2:1439741/100‑1 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                         <  2:1419316/100‑1 (MQ=255)                                                    cGTAAAACCGTCACGCCAGCCAGATTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                         >  2:1215812/1‑100 (MQ=255)                                                     gTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGt                                        >  1:1344836/1‑100 (MQ=255)                                                      tAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGta                                       >  1:43624/1‑100 (MQ=255)                                                       aaaaCCGTCACGCCAGCCAGTTTGTGGCATCAGCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtat                                      <  2:376230/100‑1 (MQ=255)                                                          aCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAAc                                   <  1:1028006/100‑1 (MQ=255)                                                            cGTCACGCCAGCCAGTTTGCGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGCGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACtt                                 >  2:2479484/1‑100 (MQ=255)                                                             gTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTa                                >  1:978490/1‑100 (MQ=255)                                                             gTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTa                                <  1:1694691/100‑1 (MQ=255)                                                               cACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAAc                              >  1:1148307/1‑100 (MQ=255)                                                               cACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAAc                              >  2:1451127/1‑100 (MQ=255)                                                                  gccagccaGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCa                           >  2:268443/1‑100 (MQ=255)                                                                  gccagccaGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCa                           >  2:1279127/1‑100 (MQ=255)                                                                   ccagccagTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAg                          >  1:1308836/1‑100 (MQ=255)                                                                   ccagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAg                          >  2:2080675/1‑100 (MQ=255)                                                                    cagccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGAAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGt                         <  2:342066/100‑1 (MQ=255)                                                                       ccagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGc                      >  1:2441524/1‑100 (MQ=255)                                                                        cagTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                     <  2:2165731/100‑1 (MQ=255)                                                                        cagTTTGTGGCATCATCCGGGGTTGGTGTTGAGTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                     <  2:2441510/100‑1 (MQ=255)                                                                        cagTTTGTGGCATCACCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGcc                     <  1:2524078/100‑1 (MQ=255)                                                                         agTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGccc                    <  1:269855/100‑1 (MQ=255)                                                                         agTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGccc                    >  2:160780/1‑100 (MQ=255)                                                                         agTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGccc                    <  2:2056451/100‑1 (MQ=255)                                                                          gTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCct                   <  1:2017779/100‑1 (MQ=255)                                                                           tttGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCctc                  <  1:1448010/100‑1 (MQ=255)                                                                               tgGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGc              >  1:1858385/1‑100 (MQ=255)                                                                                ggCATCATCCGTGGTTAGTTTCGATTCATGCGGAGGGAGTGAAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCctc                  <  2:1495571/95‑1 (MQ=255)                                                                                ggCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGcc             >  1:686612/1‑100 (MQ=255)                                                                                ggCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGcc             >  2:2030276/1‑100 (MQ=255)                                                                                    tcatcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGt         <  1:2510057/100‑1 (MQ=255)                                                                                      atcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGGATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc       >  2:2143575/1‑100 (MQ=255)                                                                                       tcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCa      >  2:903418/1‑100 (MQ=255)                                                                                           ggTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGGCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCcg  <  2:2394483/100‑1 (MQ=255)                                                                                             |                                                                                                  GCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAGCCG  >  NC_000913/1946698‑1946887

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