Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,804,063	C→G	6.0%	intergenic (‑104/‑40)	PP_3361 ← / → PP_3362	aminotransferase/AsnC family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,804,063	0	C	G	6.0%	209.5 / 3.9	69	intergenic (‑104/‑40)	PP_3361/PP_3362	aminotransferase/AsnC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base G (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

TCGCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAA  >  NC_002947/3803916‑3804205                                                                                                                                                    |                                                                                                                                                tCGCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCg                                                                                                                                               >  1:925013/1‑150 (MQ=255)   gCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑Ac                                                                                                                                            >  1:253339/1‑150 (MQ=255)   gCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑Ac                                                                                                                                            <  2:1618573/150‑1 (MQ=255)   gCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCACAACAGGCATTATGCCGGTTCCG‑Ac                                                                                                                                            <  1:283954/150‑1 (MQ=255)         cGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATTTCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCAt                                                                                                                                      >  2:2300848/1‑150 (MQ=255)         cGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCAt                                                                                                                                      >  1:2162395/1‑150 (MQ=255)            aaaCCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAAt                                                                                                                                   <  1:966406/150‑1 (MQ=255)             aaCCGAGGTGCCGGGTGGGAAAGGTCTGTCATCTCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATg                                                                                                                                  >  1:788476/1‑150 (MQ=255)              aCCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCTCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGa                                                                                                                                 >  1:502046/1‑150 (MQ=255)                       ccGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTCGGA‑ACCGGCATAATgcctgtt                                                                                                                            >  1:2244177/1‑140 (MQ=255)                       ccGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTCGGA‑ACCGGCATAATgcctgtt                                                                                                                            <  2:2244177/146‑7 (MQ=255)                         gggTGGGAAAGGTCTGTCATATCCTCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑a                                                                                                                                             <  1:1985390/127‑1 (MQ=255)                         gggTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑a                                                                                                                                             >  2:1985390/1‑127 (MQ=255)                          ggTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTAt                                                                                                                     >  1:1510918/1‑150 (MQ=255)                              gaaaaGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAg                                                                                                                 <  2:654258/148‑1 (MQ=255)                               gAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAgg                                                                                                                <  2:502046/150‑1 (MQ=255)                                    gccTGTCATATCATCGGGTTCCGGGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACCCCCTGCATACGGACGGGCGCAACATGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGAGCTg                                                                                                           <  2:1602769/148‑1 (MQ=255)                                    gTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTg                                                                                                           >  2:1312048/1‑150 (MQ=255)                                       tGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCaa                                                                                                        >  2:308828/1‑150 (MQ=255)                                          cATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGc                                                                                                     <  1:180791/150‑1 (MQ=255)                                            tatCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGt                                                                                                   <  2:925013/150‑1 (MQ=255)                                              tcatcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAgggg                                                                                                              >  1:1791107/1‑137 (MQ=255)                                              tcatcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAgggg                                                                                                              <  2:1791107/137‑1 (MQ=255)                                              tcatcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGa                                                                                                 >  2:1954129/1‑150 (MQ=255)                                                atcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGAcc                                                                                               >  2:416706/1‑150 (MQ=255)                                                 tcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATCATGACGGCAGGCTTATCGAGGGGCTGCGATGCGTGACCt                                                                                              >  2:724688/1‑150 (MQ=255)                                                 tcaGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCAACCACC‑GGCGGCATAATGACGGCAAGCTTAGCAAAGGGCTGCAATGCGTGCCCt                                                                                              >  1:1970026/1‑150 (MQ=255)                                                         cGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTa                                                                                      <  1:2373361/150‑1 (MQ=255)                                                           aGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTaaa                                                                                    <  1:1312048/150‑1 (MQ=255)                                                           aGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTaaa                                                                                    <  1:1954129/150‑1 (MQ=255)                                                                   gATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGa                                                                                                 <  1:34502/129‑1 (MQ=255)                                                                   gATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGa                                                                                                 >  2:34502/1‑129 (MQ=255)                                                                   gATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCga                                                                            >  1:306810/1‑150 (MQ=255)                                                                   gATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGAGACGCGTGACCTGGACGCTAAAGATCGCga                                                                            >  2:2007677/1‑150 (MQ=255)                                                                        tttttCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAg                                                                                                                 <  2:2294073/108‑1 (MQ=255)                                                                        tttttCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAg                                                                                                                 >  1:2294073/1‑108 (MQ=255)                                                                        tttttCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGG‑TCGGAACCGGCATAAt                                                                                                                                   >  1:2096768/1‑90 (MQ=38)                                                                        tttttCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGG‑TCGGAACCGGCATAAt                                                                                                                                   <  2:2096768/90‑1 (MQ=38)                                                                        tttttCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGACATTATGCCCGTTTCG‑ACCGGCATAATGATGGCATGCTTATCTAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTc                                                                       <  2:270627/150‑1 (MQ=255)                                                                             cAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAg                                                                  <  1:416706/150‑1 (MQ=255)                                                                              aaTAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAgggg                                                                                                              >  2:2011235/1‑105 (MQ=255)                                                                              aaTAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAgggg                                                                                                              <  1:2011235/105‑1 (MQ=255)                                                                              aaTAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAgg                                                                 >  1:2347396/1‑150 (MQ=255)                                                                                   cTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACtt                                                            >  2:214234/1‑150 (MQ=255)                                                                                      gTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGa                                                                   >  1:12408/1‑140 (MQ=255)                                                                                      gTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGa                                                                   <  2:12408/140‑1 (MQ=255)                                                                                              ccgccgAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGAcgc                                                 >  1:1826422/1‑150 (MQ=255)                                                                                                      aCGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGt                                                                                                   >  2:1204998/1‑92 (MQ=255)                                                                                                      aCGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGt                                                                                                   <  1:1204998/92‑1 (MQ=255)                                                                                                      aCGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTg                                         >  2:694512/1‑150 (MQ=255)                                                                                                         gACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCtt                                                                     >  2:790965/1‑119 (MQ=255)                                                                                                         gACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCtt                                                                     <  1:790965/119‑1 (MQ=255)                                                                                                           cacCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTc                                    <  2:253339/150‑1 (MQ=255)                                                                                                           cacCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTc                                    >  2:1742246/1‑150 (MQ=255)                                                                                                            acCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGc                                                                                                     >  2:798197/1‑84 (MQ=255)                                                                                                            acCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGc                                                                                                     <  1:798197/84‑1 (MQ=255)                                                                                                             cccTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATATCGACGTTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCaa                                  <  2:1380782/150‑1 (MQ=255)                                                                                                                   aTACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGNGtt                            <  2:1510918/150‑1 (MQ=255)                                                                                                                    tACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTa                           >  2:2058107/1‑150 (MQ=255)                                                                                                                        gTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgccg                       >  1:277813/1‑150 (MQ=255)                                                                                                                         tCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAgccgg                      >  1:1524066/1‑150 (MQ=255)                                                                                                                               ggaaccggCATAATGCCTGTTGCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCAt                >  1:49814/7‑150 (MQ=255)                                                                                                                                 aaCAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAg              <  1:308828/150‑1 (MQ=255)                                                                                                                                 aaCAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAg              >  1:69123/1‑150 (MQ=255)                                                                                                                                  aCAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAgg             <  1:694512/150‑1 (MQ=255)                                                                                                                                         ttATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGTGCGCAAGGCGTTAGCCGGCCGCATAGGCCTTTcc      <  1:2300848/150‑1 (MQ=255)                                                                                                                                         ttATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTcc      >  2:945715/1‑150 (MQ=255)                                                                                                                                             gCCTGTTGCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTACCGaa  <  2:49814/150‑1 (MQ=255)                                                                                                                                             gCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGaa  <  2:69123/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                TCGCTTGCCGAAAACCGAGGTGCCGGGTGGGAAAGGTCTGTCATATCATCAGGTTCCGAGCGCGCTGATTGTTTTTCAATAGCTTGTCATGCGCCGCCGAAACGGACACCCTGCATACGGTCGGGCGCAACAGGCATTATGCCGGTTCCG‑ACCGGCATAATGACGGCAAGCTTATCAAGGGGCTGCAATGCGTGACCTGGACGCTAAAGATCGCGAGATTCTTGAGGTACTTTCGAAGGACGCGCGTATTGCGCTCAAGGCGTTAGCCGGCCGCATAGGCCTTTCCCGAA  >  NC_002947/3803916‑3804205

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