Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,474,862	C→G	G214G (GGC→GGG)	ydbC →	putative NAD(P)‑binding oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,474,862	0	C	G	98.5%	243.9 / ‑2.6	65	G214G (GGC→GGG)	ydbC	putative NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (41/23);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.73e-01

CCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGA  >  NC_000913/1474716‑1475008                                                                                                                                                   |                                                                                                                                                   ccGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGttt                                                                                                                                                 >  1:662607/1‑150 (MQ=255) ccGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTa                                                                                                                                                <  1:620188/151‑1 (MQ=255)     cGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACAcc                                                                                                                                            <  2:950033/151‑1 (MQ=255)     cGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACAcc                                                                                                                                            >  1:130835/1‑151 (MQ=255)        aGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCt                                                                                                                                         >  1:364247/1‑151 (MQ=255)           ttGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGc                                                                                                                                       >  1:872728/1‑150 (MQ=255)               agagGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCg                                                                                                                                  >  2:486706/1‑151 (MQ=255)                     gcgcAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCcac                                                                                                                             >  1:585787/1‑150 (MQ=255)                     gcgcAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCcac                                                                                                                             >  2:895444/1‑150 (MQ=255)                      cgcAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTATGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCacac                                                                                                                           >  1:928281/1‑151 (MQ=255)                           gATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTcc                                                                                                                      >  2:656108/1‑151 (MQ=255)                            aTTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCg                                                                                                                     <  1:117709/151‑1 (MQ=255)                            aTTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCg                                                                                                                     >  2:338005/1‑151 (MQ=255)                               gCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGAt                                                                                                                   <  1:457002/150‑1 (MQ=255)                               gCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATg                                                                                                                  >  2:1149462/1‑151 (MQ=255)                                ccGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGt                                                                                                                 >  1:322150/1‑151 (MQ=255)                                 cGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGtt                                                                                                                >  1:474980/1‑151 (MQ=255)                                  gAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTtg                                                                                                               >  2:904878/1‑151 (MQ=255)                                  gAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTtg                                                                                                               >  1:48503/1‑151 (MQ=255)                                  gAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTtg                                                                                                               <  2:759893/151‑1 (MQ=255)                                    aaTTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTtgct                                                                                                             >  1:952874/1‑151 (MQ=255)                                     aTTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTtgctg                                                                                                            <  2:952930/151‑1 (MQ=255)                                       tGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTtgctgc                                                                                                           <  2:334143/150‑1 (MQ=255)                                         tCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCg                                                                                                          >  2:310469/1‑149 (MQ=255)                                            gtgtgCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAg                                                                                                        <  2:36867/148‑1 (MQ=255)                                            gtgtgCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCt                                                                                                     >  1:982227/1‑151 (MQ=255)                                             tgtgCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTg                                                                                                    >  1:42605/1‑151 (MQ=255)                                               tgCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTggg                                                                                                  >  1:866034/1‑151 (MQ=255)                                                   aaaCGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGc                                                                                               >  2:316623/1‑150 (MQ=255)                                                     aCGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAac                                                                                            >  2:92669/1‑151 (MQ=255)                                                       gAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAacac                                                                                          >  2:295519/1‑151 (MQ=255)                                                          tACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCa                                                                                        >  2:512707/1‑150 (MQ=255)                                                           acaacaTCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAAt                                                                                      >  1:473240/1‑151 (MQ=255)                                                           acaacaTCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAAt                                                                                      >  2:276853/1‑151 (MQ=255)                                                                aTCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAgg                                                                                 <  2:792909/151‑1 (MQ=255)                                                                    cgcACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGgc                                                                             >  1:1162618/1‑151 (MQ=255)                                                                      cACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGgcgc                                                                           <  2:386447/151‑1 (MQ=255)                                                                      cACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGgcgc                                                                           <  1:312600/151‑1 (MQ=255)                                                                        ccGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTg                                                                         <  2:149933/151‑1 (MQ=255)                                                                          gTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGc                                                                       <  1:1148637/151‑1 (MQ=255)                                                                          gTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGc                                                                       <  2:761360/151‑1 (MQ=255)                                                                           tGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCg                                                                      >  2:535098/1‑151 (MQ=255)                                                                            gCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGAGGCGCTGGCg                                                                      >  1:219729/1‑150 (MQ=255)                                                                               gatgatGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTg                                                                    >  1:753249/1‑149 (MQ=255)                                                                                  gatgCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGt                                                               >  1:1210576/1‑151 (MQ=255)                                                                                     gCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTAc                                                            <  2:1187237/151‑1 (MQ=255)                                                                                     gCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTAc                                                            >  1:450166/1‑151 (MQ=255)                                                                                              gatGCTTTGGCCCACGATGGTATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCAc                                                   <  1:380928/151‑1 (MQ=255)                                                                                                      ggCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATAtt                                           <  1:499335/151‑1 (MQ=255)                                                                                                         ccACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTg                                        >  1:1110474/1‑151 (MQ=255)                                                                                                               tGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGAt                                   <  1:929409/150‑1 (MQ=255)                                                                                                               tGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATc                                  >  2:566669/1‑151 (MQ=255)                                                                                                                     tGCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAggg                            >  2:472499/1‑151 (MQ=255)                                                                                                                      gCCTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAggg                            >  2:869527/1‑150 (MQ=255)                                                                                                                        cTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACg                         >  2:151992/1‑151 (MQ=255)                                                                                                                        cTACGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACg                         >  1:435226/1‑151 (MQ=255)                                                                                                                           cGTGCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGt                        >  1:678002/1‑149 (MQ=255)                                                                                                                            gTGCCGTTCTTCCCGCTCGGGGGGTTTAAACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGAGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTc                       >  1:158691/1‑149 (MQ=255)                                                                                                                              gCCGTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTc                    <  1:328674/150‑1 (MQ=255)                                                                                                                                 gTTCTTCCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCg                   <  1:1101766/148‑1 (MQ=255)                                                                                                                                   tcttcCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTgc              >  1:187344/1‑151 (MQ=255)                                                                                                                                   tcttcCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTg               <  1:956406/150‑1 (MQ=255)                                                                                                                                      tcCCGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTgcgc            <  1:491010/150‑1 (MQ=255)                                                                                                                                         cGCTCGGGGGGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCAttt        <  1:486692/151‑1 (MQ=255)                                                                                                                                              ggggggTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACggg   >  2:477010/1‑151 (MQ=255)                                                                                                                                               ggggGTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGga  <  1:895389/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCGACGCAGGTTGCAGAGGCGCGCAAGATTGCCGAAATTGTCTGTGTGCAAAACGAATACAACATCGCGCACCGTGCTGATGATGCAATGATTGATGCTTTGGCCCACGATGGCATTGCCTACGTGCCGTTCTTCCCGCTCGGGGGCTTTACACCGCTGCAATCGTCCACACTTTCCGATGTTGCTGCGAGCCTGGGTGCAACACCAATGCAGGTGGCGCTGGCGTGGCTGTTACAGCGTTCACCGAATATTTTGCTGATCCCAGGGACGTCTTCGGTTGCGCATTTACGGGA  >  NC_000913/1474716‑1475008

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