Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	233,833	C→G	9.5%	intergenic (+21/+42)	PP_0180 → / ← PP_0181	cytochrome c family protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	233,833	0	C	G	9.5%	178.8 / 8.4	63	intergenic (+21/+42)	PP_0180/PP_0181	cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/34);  new base G (2/4);  total (25/38)
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.17e-01

AGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGC  >  NC_002947/233685‑233980                                                                                                                                                     |                                                                                                                                                    aGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCcac                                                                                                                                                     <  6:334276/149‑1 (MQ=255) aGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCcac                                                                                                                                                     <  3:352513/149‑1 (MQ=255)        ctggctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGAcc                                                                                                                                              >  6:133893/1‑149 (MQ=255)         tggctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGATAAGCGCAGTGCTGCGGCCCACGGACTGGAGCCTGCGCCGTCCCGGTGGGACcg                                                                                                                                             >  8:231818/1‑149 (MQ=255)           gctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACcg                                                                                                                                             <  2:30846/147‑1 (MQ=255)            ctggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCCGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCgcgc                                                                                                                                          <  5:314032/149‑1 (MQ=255)             tggGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGGGCCACAGGGGCcgccc                                                                                                                                          >  8:185382/1‑146 (MQ=255)               gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCTTTTAGGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTCGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                         <  7:231818/147‑1 (MQ=255)               gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                         <  2:214183/147‑1 (MQ=255)               gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                         <  5:95229/147‑1 (MQ=255)               gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAg                                                                                                                                       >  8:99594/1‑149 (MQ=255)               gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGCAGCCTGCGCGGTCCCACAGGGACCGCGCaa                                                                                                                                        <  1:31229/148‑1 (MQ=255)                      aTGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCa                                                                                                                                 <  4:289375/148‑1 (MQ=255)                              gggtaagcGCTGTCGGCGCATGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGGTCGGGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTgc                                                                                                                         <  6:208631/141‑1 (MQ=255)                               ccTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTgcg                                                                                                                        <  4:145974/148‑1 (MQ=255)                                cTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACACGGGCCGCCCAAGGTCTCACATCATtcgc                                                                                                                       >  6:179054/1‑148 (MQ=255)                                   cTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGc                                                                                                                   <  6:303990/149‑1 (MQ=255)                                      gCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGc                                                                                                                   <  6:176179/146‑1 (MQ=255)                                             gcgcAAGCGTTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCag                                                                                                                                        <  5:92910/118‑2 (MQ=255)                                              cgcAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCGCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAg                                                                                                        <  5:140193/149‑1 (MQ=255)                                                   gCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctgg                                                                                                   <  3:129235/149‑1 (MQ=255)                                                     gggggCCGTGCTGGCCATGGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGGGCGGTCCCACTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAggctggct                                                                                                 <  8:243806/146‑1 (MQ=255)                                                       tGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCtgt                                                                                               <  7:185382/149‑1 (MQ=255)                                                          ccTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGgc                                                                                            <  7:99594/149‑1 (MQ=255)                                                              gCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGcc                                                                                         <  5:133893/148‑1 (MQ=255)                                                               cTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                       <  5:179054/149‑1 (MQ=255)                                                                tGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCAAAGGGACCCCGCCACGTCTCAGATCATGCGCAAGCCTGATGGACAGGCTGGCGGTGGCGccgt                                                                                       >  6:136447/1‑146 (MQ=255)                                                                 ggCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                     <  3:188832/149‑1 (MQ=255)                                                                   ccATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGCTCTTGCGCCAGCCTGATGCTCAGGCTGGGGGTGGCGCCGTgag                                                                                    >  8:29276/1‑148 (MQ=255)                                                                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAgg                                                                                                                                      >  6:184949/1‑97 (MQ=255)                                                                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAgg                                                                                                                                      <  5:184949/97‑1 (MQ=255)                                                                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                       >  2:294890/1‑144 (MQ=255)                                                                    cATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCAGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                       <  1:294890/144‑1 (MQ=255)                                                                         tgtTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTc                                                                             >  8:9830/1‑149 (MQ=255)                                                                         ggggccTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACCGGAGCCTGCGCGGTCCCACACGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGACTGTTc                                                                             <  8:151170/145‑1 (MQ=255)                                                                            tCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCCGTCCCACAGGGACCGCGCAAGGTGGCAGAGAGTTCGCAAGCCTGATGCTCAGGCGGGCGGTGGGGGCTTGGGAATGTTgtgt                                                                          >  3:277976/1‑149 (MQ=255)                                                                             ccTTTACGGCCGCTCGGGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGt                                                                          <  5:136447/148‑1 (MQ=255)                                                                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCCCAGGGACCGCGCAAGGTGGGAGATCTTGCGCAAGCCGGATGCTCAGGCTGGCTGTGGCGCCTTGAGCATGTTCTGTAg                                                                        >  5:292294/1‑149 (MQ=255)                                                                              cTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAg                                                                        >  8:58635/1‑149 (MQ=255)                                                                                     gCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACACGGCCCGCGCAAGGTCTCAGATCTTGCGCAAGCCAGCTGCTCACGCAGGCTGGGGCGCCGTGACAATTTTATGTACCAATTg                                                                  >  5:202709/1‑148 (MQ=255)                                                                                       cGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGtt                                                                              <  5:283806/134‑1 (MQ=255)                                                                                       cGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACACGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCAGATGCTCAGGCTGGCTGTGGCGCCATGAGaaattt                                                                              >  6:283806/1‑130 (MQ=255)                                                                                               gctggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCCGTGCCCCACGGCCCCCGCAAGGTGGGATATATTTGGCAACACAGATGATAAGGCTGGGGGTGGCGCCGGGGGTATGTTGTGGAGCAAATGGAGACgttga                                                        >  4:7504/1‑146 (MQ=255)                                                                                                ctggcCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAAc                                                      >  1:39161/1‑149 (MQ=255)                                                                                                    ccGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACg                                                     <  2:63771/146‑1 (MQ=255)                                                                                                    ccGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACg                                                     >  1:63771/1‑146 (MQ=255)                                                                                                      gAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACACGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGt                                                 >  5:307562/1‑148 (MQ=255)                                                                                                             gcAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGAccc                                         <  6:307562/149‑1 (MQ=255)                                                                                                                   ctgcGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCt                                    <  1:327724/148‑1 (MQ=255)                                                                                                                       ggCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTTCGCAAACCTGATGCTCAGGCAGGCTGTGGCGCCGTGAGAATGTTCTGTATCAATTTGACACGTTAAACGGCGTCATGACCCAAGCTGGt                                 >  6:161674/1‑147 (MQ=255)                                                                                                                          cAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGGGGCGCCGTGAGAATGTTCTGTAGCAATTGGCCACTTTCAACGGCATCATGACCCAGGCTGGTCAGGt                            >  5:169843/1‑149 (MQ=255)                                                                                                                          cAACTGACTGGACCCTGCGCGGTCCCACAGCGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGAAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGt                            <  6:292294/149‑1 (MQ=255)                                                                                                                           aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                       >  2:40646/1‑89 (MQ=255)                                                                                                                           aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                       <  1:40646/89‑1 (MQ=255)                                                                                                                           aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTa                           <  5:161674/149‑1 (MQ=255)                                                                                                                                  tgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                     <  7:138591/84‑1 (MQ=255)                                                                                                                                  tgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCa                    <  7:304296/149‑1 (MQ=255)                                                                                                                                      gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgccg                       <  4:47471/142‑1 (MQ=255)                                                                                                                                      gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgccg                       >  3:47471/1‑142 (MQ=255)                                                                                                                                       cttGCGCGGTCCCTGTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCgg               >  7:242594/3‑149 (MQ=255)                                                                                                                                        cTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGc              <  7:29276/149‑1 (MQ=255)                                                                                                                                                tCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGt      <  2:39161/149‑1 (MQ=255)                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTc     <  2:316488/148‑1 (MQ=255)                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTc     >  1:316488/1‑148 (MQ=255)                                                                                                                                                  ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCa    <  7:58635/149‑1 (MQ=255)                                                                                                                                                   cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg   >  7:93112/1‑149 (MQ=255)                                                                                                                                                    acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc  <  8:93112/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    AGTTCGACTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGC  >  NC_002947/233685‑233980

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