Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,634,457	T→G	6.1%	intergenic (+33/+350)	putA → / ← PP_5718	proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,634,457	0	T	G	6.1%	201.9 / 4.4	66	intergenic (+33/+350)	putA/PP_5718	proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/30);  new base G (2/2);  total (34/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 = 5.68e-01

GCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGG  >  NC_002947/5634313‑5634601                                                                                                                                                 |                                                                                                                                                  gCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGc                                                                                                                                               >  1:212592/1‑149 (MQ=255)     tCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCtctc                                                                                                                                           <  5:134801/149‑1 (MQ=255)         gCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCttt                                                                                                                                       <  2:212592/149‑1 (MQ=255)         gCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttc                                                                                                                                       <  6:102904/149‑2 (MQ=255)          cGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttca                                                                                                                                      <  8:153196/149‑3 (MQ=255)            gggATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCGttttt                                                                                                                                     >  8:110228/1‑148 (MQ=255)              gATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTAt                                                                                                                                  <  1:168172/149‑1 (MQ=255)                   ccAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTg                                                                                                                             >  7:284318/1‑149 (MQ=255)                   ccAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTg                                                                                                                             >  1:568/1‑149 (MQ=255)                    cAGATTGCGCTGGAGCGCCTGGTCAACGAGCGCGCGGTGAGCGTGAACACTACAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGAATGAGGGTTGTTGAAGAAGAAGGAGATCCTAGGCTATCCGTTGTTATGCGTg                                                                                                                             >  7:198330/1‑148 (MQ=255)                           cgcTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCt                                                                                                                     >  2:179581/1‑149 (MQ=255)                              tGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTg                                                                                                                    >  7:137749/1‑147 (MQ=255)                                gAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGTGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGaa                                                                                                                  <  4:139019/147‑1 (MQ=255)                                gAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAg                                                                                                                 >  5:84664/1‑148 (MQ=255)                                   cgcCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgg                                                                                                             <  8:143147/149‑1 (MQ=255)                                      cTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAg                                                                                                              >  4:94930/1‑145 (MQ=255)                                      cTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAg                                                                                                              <  3:94930/145‑1 (MQ=255)                                         gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgggg                                                                                                           <  8:132188/145‑1 (MQ=255)                                         gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAgggg                                                                                                           >  7:132188/1‑145 (MQ=255)                                         gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                        <  8:137749/148‑1 (MQ=255)                                               gAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                      >  3:252506/1‑112 (MQ=255)                                               gAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCtttt                                                                                                                                      <  4:252506/112‑1 (MQ=255)                                                        gTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTCAAGCAGGGGCTGGACTAACA‑GGTGCTCa                                                                                       <  1:179581/149‑1 (MQ=255)                                                              gTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGaa                                                                                 >  1:212406/1‑149 (MQ=255)                                                               tGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAAc                                                                                <  2:568/149‑1 (MQ=255)                                                                      tGCAGCGGCGGGTGGTAATGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCt                                                                                                                     >  3:90133/1‑106 (MQ=255)                                                                      tGCAGCGGCGGGTGGTAATGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCt                                                                                                                     <  4:90133/106‑1 (MQ=255)                                                                          gcggcgGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAAc                                                                                >  3:31136/1‑138 (MQ=255)                                                                          gcggcgGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAAc                                                                                <  4:31136/138‑1 (MQ=255)                                                                                 gtggtAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGGACGTGGTCATGCCAACGGt                                                               >  2:237777/1‑148 (MQ=255)                                                                                   ggtAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTACTTTTTTATGCGTGCATGTCCTGAAGAAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGTTCATTCCAACGGTGaa                                                            >  7:109986/1‑148 (MQ=255)                                                                                      aaCGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTg                                                         >  6:141711/1‑149 (MQ=255)                                                                                       aCGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTgg                                                        <  6:84664/149‑1 (MQ=255)                                                                                        cGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTggg                                                       <  2:275452/149‑1 (MQ=255)                                                                                         gCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                      >  1:51306/1‑149 (MQ=255)                                                                                         gCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                      <  7:110228/149‑1 (MQ=255)                                                                                           cAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgtgt                                                    >  4:266749/1‑149 (MQ=255)                                                                                            aGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                      >  4:213672/1‑99 (MQ=255)                                                                                            aGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                      <  3:213672/99‑1 (MQ=255)                                                                                            aGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGtc                                                   <  3:266749/149‑1 (MQ=255)                                                                                                     aCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACACGGTGCTCACCAGGACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                            <  1:237777/148‑1 (MQ=255)                                                                                                      ccATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTACTGAAGCAGGGGCTGGACTAACA‑GGTGCTAACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgttgt                                          >  7:70952/1‑148 (MQ=255)                                                                                                       cATTGTCTAATGGCTGTTGGTTGTTGAAGAAGAAGGAGAGGCTCGGCTCTACTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGTACTAACA‑TGTGCta                                                                                        <  4:202598/101‑2 (MQ=255)                                                                                                                   gcttcgggTTGATGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCGGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATcccc                            <  8:284318/144‑1 (MQ=255)                                                                                                                    cTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCt                           <  2:51306/149‑1 (MQ=255)                                                                                                                        gggTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑g                                                                                              >  4:284805/1‑78 (MQ=255)                                                                                                                        gggTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑g                                                                                              <  3:284805/78‑1 (MQ=255)                                                                                                                        gggTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGa                        >  7:17172/1‑148 (MQ=255)                                                                                                                         ggTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCaa                                                                   >  3:245476/1‑104 (MQ=255)                                                                                                                         ggTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCaa                                                                   <  4:245476/104‑1 (MQ=255)                                                                                                                         ggTTGTTGAAGAAGAAGGAGAGACTCGGCGATCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGaa                       <  8:109986/148‑1 (MQ=255)                                                                                                                          gttgttGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTg                                                              >  3:260627/1‑108 (MQ=255)                                                                                                                          gttgttGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTg                                                              <  4:260627/108‑1 (MQ=255)                                                                                                                          gttgttGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAAt                     >  8:269521/1‑149 (MQ=255)                                                                                                                              ttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATtctc                 <  2:212406/149‑1 (MQ=255)                                                                                                                               tgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCg                >  7:95367/1‑149 (MQ=255)                                                                                                                                cataaaaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCgg               >  5:119450/8‑149 (MQ=255)                                                                                                                                   gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGtt                                   <  5:256027/126‑1 (MQ=255)                                                                                                                                   gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGtt                                   >  6:256027/1‑126 (MQ=255)                                                                                                                                     aaaagGAGAGCCGAGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGcc          >  5:196811/3‑149 (MQ=255)                                                                                                                                      gaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCa                                                                                       <  8:128951/71‑1 (MQ=255)                                                                                                                                      gaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCa                                                                                       >  7:128951/1‑71 (MQ=255)                                                                                                                                      gaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                            >  7:278648/1‑114 (MQ=255)                                                                                                                                        agGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAAGTCTCGGCTGCCTgg       <  7:269521/149‑1 (MQ=255)                                                                                                                                          gagaGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGc     >  8:12848/1‑149 (MQ=255)                                                                                                                                             agCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCgg  >  7:197015/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GCTGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATGACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACA‑GGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGG  >  NC_002947/5634313‑5634601

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