Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,391,997	A→G	9.4%	intergenic (+22/‑93)	aspS → / → PP_1214	aspartate‑‑tRNA ligase/transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,391,997	0	A	G	9.4%	177.2 / 7.7	64	intergenic (+22/‑93)	aspS/PP_1214	aspartate‑‑tRNA ligase/transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (2/4);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

GATCGCGTTCCCGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCC  >  NC_002947/1391855‑1392146                                                                                                                                                |                                                                                                                                                      gATCGCGTTCCCGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGACCTGCACATCCG‑ACTGCGCCAACAGCTCAAGCTCGAGTAAGCAGTCCCCGGGCGCATCCATGTGGATgc                                                                                                                                                <  2:356825/150‑1 (MQ=255)        ttCCCGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCATGTGGATGCGcccggg                                                                                                                                         <  1:556778/150‑5 (MQ=255)         tCCCGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGCCGCATCCACATGGATGCGCCTTTGc                                                                                                                                        <  1:1271406/150‑1 (MQ=255)           ccGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGt                                                                                                                                      >  2:113601/1‑150 (MQ=255)           ccGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAACCAGGCCCCGGGCGCATCCACATGGATGCGCCTTTGCGt                                                                                                                                      >  2:1206541/1‑150 (MQ=255)            cGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGtt                                                                                                                                     >  2:719587/1‑150 (MQ=255)            cGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGtt                                                                                                                                     <  2:494399/150‑1 (MQ=255)                     agagCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAg                                                                                                                            <  2:1619905/150‑1 (MQ=255)                            tGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCGAAGGCCCTGCGACAACAGCTCACTCGTGCTGC‑CGATCAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGtttttt                                                                                                                     <  1:1103704/150‑1 (MQ=255)                              cATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTtcc                                                                                                                   <  1:719587/150‑1 (MQ=255)                                   gTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTtcccgtcc                                                                                                              >  2:638060/1‑150 (MQ=255)                                    tCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTtcccgtccc                                                                                                             >  2:1607125/1‑150 (MQ=255)                                        gACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGccc                                                                                                         >  1:1557405/1‑150 (MQ=255)                                         aCCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGcccc                                                                                                        >  2:1415769/1‑150 (MQ=255)                                           ccAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGGCAGTACAGACCAAGGGCCAATAGGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGt                                                                                                      >  1:1019430/1‑150 (MQ=255)                                           ccAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGt                                                                                                      >  1:470006/1‑150 (MQ=255)                                                   cTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGa                                                                                                                                                   >  2:961025/1‑97 (MQ=255)                                                   cTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGa                                                                                                                                                   <  1:961025/97‑1 (MQ=255)                                                      gCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATggggcggg                                                                                           <  1:94441/150‑1 (MQ=255)                                                            tCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGt                                                                                                                                      >  2:113420/1‑101 (MQ=255)                                                            tCGACGCCAAGGCCCTGCGCGAACTGCACACCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGt                                                                                                                                      <  1:113420/101‑1 (MQ=255)                                                                 gCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGa                                                                                >  1:1460447/1‑150 (MQ=255)                                                                 gCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGa                                                                                <  1:390078/150‑1 (MQ=255)                                                                  ccAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGc                                                                                                           <  2:748137/122‑1 (MQ=255)                                                                  ccAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGc                                                                                                           >  1:748137/1‑122 (MQ=255)                                                                    aaGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCATGTGGATGCGcccgg                                                                                                                                          <  1:1455748/89‑4 (MQ=39)                                                                    aaGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCATGTGGATGCGcccgg                                                                                                                                          >  2:1455748/1‑86 (MQ=39)                                                                       gCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCaaa                                                                          <  1:638060/150‑1 (MQ=255)                                                                         ccTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAgg                                                                        >  2:472868/1‑150 (MQ=255)                                                                         ccTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAgg                                                                        >  2:983369/1‑150 (MQ=255)                                                                          cTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAg                                                                                                                            <  1:1140125/97‑1 (MQ=255)                                                                          cTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAg                                                                                                                            >  2:1140125/1‑97 (MQ=255)                                                                          cTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTAATGGGGCGGGGCCGTTTGTGATTCAAAGGa                                                                       >  1:750482/1‑150 (MQ=255)                                                                          cTGCGCGAACTGCACATCCG‑ACTGCGCGAACAAACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGa                                                                       <  1:1591655/150‑1 (MQ=255)                                                                             cgcgAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGTGGGGCCGTTTGTGATTCAAAGGATTc                                                                    <  1:976284/150‑1 (MQ=255)                                                                              gcgAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGTCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCTTGGGGCGGGGCCGTTTGTGATTCAAAGGATTCg                                                                   <  1:1206541/150‑1 (MQ=255)                                                                                             g‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTg                                                    >  1:1151430/1‑150 (MQ=255)                                                                                               aCTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTtcccgtccc                                                                                                             <  1:1459145/92‑1 (MQ=255)                                                                                               aCTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTtcccgtccc                                                                                                             >  2:1459145/1‑92 (MQ=255)                                                                                               aCTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTgg                                                   <  2:1151430/150‑1 (MQ=255)                                                                                               aCTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTgg                                                   >  1:1478565/1‑150 (MQ=255)                                                                                                   cgcgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCAt                                               <  2:1460447/150‑1 (MQ=255)                                                                                                     cgAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTc                                             <  2:1478565/150‑1 (MQ=255)                                                                                                            aCCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTgg                                      >  2:80882/1‑150 (MQ=255)                                                                                                              cgaGGTCGAGTAAGCAGTCCCAGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGc                                    >  2:477219/3‑150 (MQ=255)                                                                                                              cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTc                                                              <  1:1234248/124‑1 (MQ=255)                                                                                                              cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTc                                                              >  2:1234248/1‑124 (MQ=255)                                                                                                              cAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGc                                    <  2:1019430/150‑1 (MQ=255)                                                                                                               aaGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTg                                                    >  1:1015501/1‑133 (MQ=255)                                                                                                               aaGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTg                                                    <  2:1015501/133‑1 (MQ=255)                                                                                                                aGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCg                                   >  2:854911/1‑149 (MQ=255)                                                                                                                aGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCg                                   <  1:854911/149‑1 (MQ=255)                                                                                                                     gAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCTTTTTTTGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATc                             <  2:470006/150‑1 (MQ=255)                                                                                                                       gTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTc                                                                    <  1:378680/109‑1 (MQ=255)                                                                                                                       gTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTc                                                                    >  2:378680/1‑109 (MQ=255)                                                                                                                         aaGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGc                         >  1:1410279/1‑150 (MQ=255)                                                                                                                            cAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTc                                                                             <  2:1203821/95‑1 (MQ=255)                                                                                                                            cAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTc                                                                             >  1:1203821/1‑95 (MQ=255)                                                                                                                                 cccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGttt                                                                                    <  2:1172560/83‑1 (MQ=255)                                                                                                                                 cccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGttt                                                                                    >  1:1172560/1‑83 (MQ=255)                                                                                                                                 cccGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCaaa                 <  1:80882/150‑1 (MQ=255)                                                                                                                                      gcgcATCCATGTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTAt                                                         <  2:841073/105‑1 (MQ=255)                                                                                                                                      gcgcATCCATGTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTAt                                                         >  1:841073/1‑105 (MQ=255)                                                                                                                                       cgcATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAgcgc           <  1:983369/150‑1 (MQ=255)                                                                                                                                             cacaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTtg                                                                                   >  1:400914/1‑72 (MQ=255)                                                                                                                                             cacaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTtg                                                                                   <  2:400914/72‑1 (MQ=255)                                                                                                                                              acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATg                                                                                                  >  2:538661/1‑56 (MQ=255)                                                                                                                                              acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATg                                                                                                  <  1:538661/56‑1 (MQ=255)                                                                                                                                              acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCaaa                                                                          <  2:327569/80‑1 (MQ=255)                                                                                                                                              acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCaaa                                                                          >  1:327569/1‑80 (MQ=255)                                                                                                                                              acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATg    >  2:468347/1‑150 (MQ=255)                                                                                                                                              acaTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATg    >  2:167720/1‑150 (MQ=255)                                                                                                                                                aTGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGcc  <  1:167720/150‑1 (MQ=255)                                                                                                                                                aTGGATGCGCCTTTGCGTTTCGACGCAGGGTTTTTCCCGTCCCGCCCCCTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGcc  >  2:1422388/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      GATCGCGTTCCCGAAAACCCAGAGCGCTGCATGCGTCATGACCCAGGCCCCTGGCATGGTCGACGCCAAGGCCCTGCGCGAACTGCACATCCG‑ACTGCGCGAACAGACCAAGGTCGAGTAAGCAGTCCCCGGGCGCATCCACATGGATGCGCCTTTGCGTTTCGACGCAGGTTTTTTCCCGTCCCGCCCCGTCATGGGGCGGGGCCGTTTGTGATTCAAAGGATTCGGAGTCGTTATGGCTGGTCATTCCAAGTGGGCGAACATCAAGCACCGCAAAGAGCGCCAGGATGCC  >  NC_002947/1391855‑1392146

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