Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,634,458	C→A	6.2%	intergenic (+34/+349)	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,458	0	C	A	6.2%	202.3 / 3.6	64	intergenic (+34/+349)	putA/PP_5718	proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base A (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.18e-01

TGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTC  >  NC_002947/5634315‑5634603                                                                                                                                                 |                                                                                                                                                  tGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGctc                                                                                                                                              >  1:367775/1‑149 (MQ=255) tGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGctc                                                                                                                                              <  3:121445/149‑1 (MQ=255)  gTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGctct                                                                                                                                             >  7:18339/1‑149 (MQ=255)          gggATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCttttt                                                                                                                                      >  6:33610/1‑148 (MQ=255)           ggATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTa                                                                                                                                    >  7:340866/1‑149 (MQ=255)             aTCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttct                                                                                                                                       >  5:94718/1‑142 (MQ=255)             aTCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttct                                                                                                                                       <  6:94718/144‑3 (MQ=255)              tCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGc                                                                                                                                 <  5:33610/149‑1 (MQ=255)               cACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTAt                                                                                                                                   >  4:45916/1‑146 (MQ=255)               cACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTAt                                                                                                                                   <  3:45916/146‑1 (MQ=255)                aCCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAATGAGAGCCTCGGCTCTCCTTTTTTATGCGt                                                                                                                               >  6:112542/1‑149 (MQ=255)                  cAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTg                                                                                                                              >  7:22737/1‑148 (MQ=255)                   aGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCa                                                                                                                            >  3:85617/1‑149 (MQ=255)                    gATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGGTCTCCTTTTTTATGCGTGCAt                                                                                                                           >  3:10487/1‑149 (MQ=255)                          gcTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTATCGCCAGCCTGATGCTTTATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCt                                                                                                                      <  1:298955/149‑1 (MQ=255)                            tGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTg                                                                                                                     >  6:130555/1‑147 (MQ=255)                              gAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGaa                                                                                                                   <  4:85617/147‑1 (MQ=255)                              gAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAg                                                                                                                  <  8:18339/148‑1 (MQ=255)                              gAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCGCTCCTTTTTTATGCGTGCATGTGCTGCAGc                                                                                                                 >  8:238727/1‑149 (MQ=255)                                gcgcCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAg                                                                                                               >  3:66977/1‑149 (MQ=255)                                       gTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTg                                                                                                         <  8:340866/148‑1 (MQ=255)                                        tgaTCGAGCGCGCGGTGAGCGTGAACACGGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGa                                                                                                       <  5:366224/147‑1 (MQ=255)                                          aTCGAGCGCGTGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGAc                                                                                                      >  8:288747/1‑148 (MQ=255)                                              cgcgcgcGGGGAGCGTGAACACTGCAGCGGCGGGTTGTAACGCCAGCCTGATG‑AGCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAACCAGGGGCTGGACTAACa                                                                                                 <  4:10487/148‑1 (MQ=255)                                               gcgcgcGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAg                                                                                                >  6:9470/1‑149 (MQ=255)                                                 gcgcgGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGt                                                                                              <  6:312586/149‑1 (MQ=255)                                                   gcgGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttc                                                                                                                                        <  7:52604/105‑2 (MQ=255)                                                   gcgGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCGAGGCTCTCCttc                                                                                                                                        >  8:52604/1‑104 (MQ=255)                                                             tGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAAc                                                                                  <  2:367775/149‑1 (MQ=255)                                                                       aGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGc                                                                        <  8:22737/149‑1 (MQ=255)                                                                                 ggtAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGAc                                                              >  8:212738/1‑149 (MQ=255)                                                                                   tAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACtt                                                            >  2:267959/1‑149 (MQ=255)                                                                                   tAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACtt                                                            <  2:53614/149‑1 (MQ=255)                                                                                       gCCAGCCTGATG‑ACCTTTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGTCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                        <  3:267758/149‑1 (MQ=255)                                                                                       gCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                        >  8:251111/1‑149 (MQ=255)                                                                                       gCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgt                                                        <  7:212738/149‑1 (MQ=255)                                                                                        ccAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGgtg                                                       >  8:177158/1‑149 (MQ=255)                                                                                               gatg‑aCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtg                                                >  4:203091/1‑149 (MQ=255)                                                                                                 tg‑aCCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                              <  7:288747/149‑1 (MQ=255)                                                                                                        ttGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTAcc                                        <  1:267959/149‑1 (MQ=255)                                                                                                           gCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGt                                                                 >  8:175953/1‑121 (MQ=255)                                                                                                           gCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGt                                                                 <  7:175953/121‑1 (MQ=255)                                                                                                             tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGc                                                                        <  7:293988/112‑1 (MQ=255)                                                                                                             tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGc                                                                        >  8:293988/1‑112 (MQ=255)                                                                                                             tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                              <  7:162659/138‑1 (MQ=255)                                                                                                             tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCtgtt                                              >  8:162659/1‑138 (MQ=255)                                                                                                             tAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGtttg                                   >  4:301436/1‑148 (MQ=255)                                                                                                                   cTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCt                             >  8:174162/1‑149 (MQ=255)                                                                                                                    tGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTg                            >  6:228628/1‑149 (MQ=255)                                                                                                                        ggTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGaaa                        >  4:15729/1‑149 (MQ=255)                                                                                                                             ttgAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATtctc                   <  5:9470/149‑1 (MQ=255)                                                                                                                               gaagaagaaGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCg                  <  7:251111/148‑1 (MQ=255)                                                                                                                                aagaagaagGCGAGCCTCGGCTCTCCGTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGGCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCgg                 <  7:238727/148‑1 (MQ=255)                                                                                                                                aagaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGc                >  8:182469/1‑149 (MQ=255)                                                                                                                                  gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTg                            >  7:122266/1‑135 (MQ=255)                                                                                                                                  gaagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTg                            <  8:122266/135‑1 (MQ=255)                                                                                                                                   aagaagGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGc             <  4:66977/149‑1 (MQ=255)                                                                                                                                     gaagGAGAGCCTCGGCTCTCCTTTTTTTTGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCt           >  6:267324/1‑149 (MQ=255)                                                                                                                                         gagaGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGc       <  5:112542/149‑1 (MQ=255)                                                                                                                                             gCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGt   <  5:228628/149‑1 (MQ=255)                                                                                                                                              ccTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTc  >  2:31527/1‑149 (MQ=255)                                                                                                                                              ccTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTc  >  2:40794/1‑149 (MQ=255)                                                                                                                                                tCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGt   <  2:363027/146‑1 (MQ=255)                                                                                                                                                tCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGt   >  1:363027/1‑146 (MQ=255)                                                                                                                                                 |                                                                                                                                                  TGTCCAGCGGGGATCACCAGATTGCGCTGGAGCGCCTGGTCATCGAGCGCGCGGTGAGCGTGAACACTGCAGCGGCGGGTGGTAACGCCAGCCTGATG‑ACCATTGGCTAATGGCTGAGGGTTGTTGAAGAAGAAGGAGAGCCTCGGCTCTCCTTTTTTATGCGTGCATGTCCTGAAGCAGGGGCTGGACTAACAGGTGCTCACCAGAACGTGGTCATGCCAACGGTGACTTGGGTGTCTTCTGTTGTTACCGTTTATCCCCTGGAAATTCTCGGCTGCCTGGGCCGGTC  >  NC_002947/5634315‑5634603

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