Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,409,800	A→G	6.2%	M375T (ATG→ACG)	stcD ←	N‑methylproline demethylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,409,800	0	A	G	6.2%	194.9 / 2.3	64	M375T (ATG→ACG)	stcD	N‑methylproline demethylase
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (2/2);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

CCGCCGGGTTTGTCCGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGG  >  NC_002947/5409659‑5409947                                                                                                                                              |                                                                                                                                                     ccgccgGGTTTGTCCGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGtt                                                                                                                                              >  2:377827/1‑150 (MQ=255)          ttGTCCGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCAAACGTGCCGCCTCCCGGGCAACCAAGAAGAACCAGATAATCACTACGTTTCGCTGTTGCACAGCCTTGTGTATTTCATGAGGCATGGTGGTTTCACGCCCg                                                                                                                                     <  2:1006161/150‑1 (MQ=255)           tGTCCGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCgg                                                                                                                                    >  2:154616/1‑150 (MQ=255)             tCCGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTc                                                                                                                                  >  2:799142/1‑150 (MQ=255)             tCCGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTc                                                                                                                                  >  1:709900/1‑150 (MQ=255)              ccGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCANAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTcg                                                                                                                                 <  1:1499111/150‑1 (MQ=255)                        gTACACCACCACCTCGTGTCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTtgtg                                                                                                                      >  2:887477/1‑150 (MQ=255)                          acacCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGa                                                                                                                    <  2:709900/150‑1 (MQ=255)                            accaccaccTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATg                                                                                                                  <  1:799142/150‑1 (MQ=255)                                       gTGGCCGCACTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCcgc                                                                                                       >  2:191351/1‑150 (MQ=255)                                         ggCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCgcgc                                                                                                     <  1:887477/150‑1 (MQ=255)                                          gCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGcc                                                                                                    >  2:701189/1‑150 (MQ=255)                                            cgcgCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTg                                                                                                  <  1:154616/150‑1 (MQ=255)                                             gcgcTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGc                                                                                                 <  1:191351/150‑1 (MQ=255)                                               gcTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGt                                                                                               >  1:1512929/1‑150 (MQ=255)                                               gcTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGt                                                                                               <  1:266269/150‑1 (MQ=255)                                               gcTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGt                                                                                               <  2:490827/150‑1 (MQ=255)                                                cTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGtt                                                                                              <  2:1512929/150‑1 (MQ=255)                                                  cGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCTTGCGTTGt                                                                                            <  1:144213/150‑1 (MQ=255)                                                     cTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGAt                                                                                                                   <  2:391789/124‑1 (MQ=255)                                                     cTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGAt                                                                                                                   >  1:391789/1‑124 (MQ=255)                                                         cacaCGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCg                                                                                     <  1:223975/150‑1 (MQ=255)                                                          acacGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCgcgc                                                                                                     >  1:1290707/1‑133 (MQ=255)                                                          acacGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCgcgc                                                                                                     <  2:1290707/133‑1 (MQ=255)                                                            acGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGCGGTTTCACGCCCGGTCGCGG‑TGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTc                                                                                  <  1:240291/150‑1 (MQ=255)                                                            acGTCCCGCCGCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCTCTTTTGCACAGCCTTGGGTATTTCATGAGGCACGGGGGCTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTc                                                                                  <  2:803525/150‑1 (MQ=255)                                                                gccgccTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAgg                                                                              <  2:1392517/150‑1 (MQ=255)                                                                   gccTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCaa                                                                           >  2:1165495/1‑150 (MQ=255)                                                                      tCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAg                                                                        >  1:447624/1‑150 (MQ=255)                                                                         aGGCCACCAGGGCCGGCGCCGATAATCGCCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTg                                                                     >  1:457524/1‑150 (MQ=255)                                                                          ggCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTaaa                                                                                         >  2:1548701/1‑129 (MQ=255)                                                                          ggCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTaaa                                                                                         <  1:1548701/129‑1 (MQ=255)                                                                              accaGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGcc                                                                >  2:1231403/1‑150 (MQ=255)                                                                                  gggCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATTGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAgg                                                                              >  2:204240/1‑132 (MQ=255)                                                                                  gggCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATTGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAgg                                                                              <  1:204240/132‑1 (MQ=255)                                                                                        gcgcCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGAAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAggg                                                      >  2:193968/1‑150 (MQ=255)                                                                                          gcCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGGGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTCAATGCGGTCCAGGCAATAGTTGGCGCCTACGCACGGGCg                                                    >  1:639356/1‑150 (MQ=255)                                                                                          gcCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCg                                                    <  1:258373/150‑1 (MQ=255)                                                                                          gcCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCCGATAAGATCGTCG‑GCATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCg                                                    <  2:875494/150‑1 (MQ=255)                                                                                             gATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGat                                                 >  2:643906/1‑150 (MQ=255)                                                                                             gATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGat                                                 <  2:457524/150‑1 (MQ=255)                                                                                                   caccacTTTGCGCTTTTGCACAGCCTGGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTc                                           >  2:63292/1‑150 (MQ=255)                                                                                                      cacTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTc                                        >  2:166812/1‑150 (MQ=255)                                                                                                        cTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCAc                                      <  1:377827/150‑1 (MQ=255)                                                                                                         tttGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCAtg                                                                                                                                                    >  2:1185193/1‑40 (MQ=255)                                                                                                         tttGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCAtg                                                                                                                                                    <  1:1185193/40‑1 (MQ=255)                                                                                                                ttttGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCcgc                                                                                                       <  2:1183022/77‑1 (MQ=255)                                                                                                                ttttGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCcgc                                                                                                       >  1:1183022/1‑77 (MQ=255)                                                                                                                ttttGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCg                              >  1:186179/1‑150 (MQ=255)                                                                                                                 tttGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCga                             <  2:1244630/150‑1 (MQ=255)                                                                                                                  ttGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCgat                            <  1:63292/150‑1 (MQ=255)                                                                                                                      acaGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATc                        >  2:1064412/1‑150 (MQ=255)                                                                                                                          ccTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTgc                    <  1:1064412/150‑1 (MQ=255)                                                                                                                              gggTATTTCATGAGGCGTGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGt                                            <  2:1371325/122‑1 (MQ=255)                                                                                                                              gggTATTTCATGAGGCGTGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGt                                            <  2:1371262/122‑1 (MQ=255)                                                                                                                              gggTATTTCATGAGGCGTGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGt                                            >  1:1371325/1‑122 (MQ=255)                                                                                                                              gggTATTTCATGAGGCGTGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGt                                            >  1:1371262/1‑122 (MQ=255)                                                                                                                                 tATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCtt                      >  2:46732/1‑141 (MQ=255)                                                                                                                                 tATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCtt                      <  1:46732/141‑1 (MQ=255)                                                                                                                                       aTGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTc                                           >  2:1007022/1‑114 (MQ=255)                                                                                                                                       aTGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTc                                           <  1:1007022/114‑1 (MQ=255)                                                                                                                                        tGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCggg      >  2:1407662/1‑150 (MQ=255)                                                                                                                                            gCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGggtcgg  <  1:643906/150‑1 (MQ=255)                                                                                                                                            gCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGggtcgg  >  1:1330700/1‑150 (MQ=255)                                                                                                                                              |                                                                                                                                                     CCGCCGGGTTTGTCCGCCAGCTCGTACACCACCACCTCGTGGCCGCGCTCGCCTGCCACACGTGCCGCCTCCAGGCCACCAGGGCCGGCGCCGATAATCACCACTTTGCGCTTTTGCACAGCCTTGGGTATTTCATGAGGCATGGTGGTTTCACGCCCGGTCGCGG‑CGTTGTGGATGCAGTAGGCCGCGCCTGCGTTGTAAATGCGGTCCAGGCAATAGTTGGCGCCTACGCAAGGGCGGATATCGTCTTCACGGCCTTCGATGATCTTGCGCACGATGTGCGGGTCGG  >  NC_002947/5409659‑5409947

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