Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	233,834	A→T	10.6%	intergenic (+22/+41)	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,834	0	A	T	10.6%	152.3 / 6.3	63	intergenic (+22/+41)	PP_0180/PP_0181	cytochrome c family protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (23/33);  new base T (2/5);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

CTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGC  >  NC_002947/233692‑233980                                                                                                                                               |                                                                                                                                                   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                                                                                                                                         <  5:95229/147‑1 (MQ=255)        gggATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCTGTGGGACCGCGCa                                                                                                                                         <  2:214183/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                                                     >  1:63771/1‑146 (MQ=255)                                                                                             ccGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACg                                                     <  2:63771/146‑1 (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                                                                                       <  1:40646/89‑1 (MQ=255)                                                                                                                    aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCtt                                                                                       >  2:40646/1‑89 (MQ=255)                                                                                                                    aaCTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTa                           <  5:161674/149‑1 (MQ=255)                                                                                                                           tgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTga                                                                                     <  7:138591/84‑1 (MQ=255)                                                                                                                           tgGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCa                    <  7:304296/149‑1 (MQ=255)                                                                                                                               gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgccg                       >  3:47471/1‑142 (MQ=255)                                                                                                                               gCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAgccg                       <  4:47471/142‑1 (MQ=255)                                                                                                                                cttGCGCGGTCCCTGTGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCgg               >  7:242594/3‑149 (MQ=255)                                                                                                                                 cTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGc              <  7:29276/149‑1 (MQ=255)                                                                                                                                         tCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGt      <  2:39161/149‑1 (MQ=255)                                                                                                                                           ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTc     >  1:316488/1‑148 (MQ=255)                                                                                                                                           ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTc     <  2:316488/148‑1 (MQ=255)                                                                                                                                           ccACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCa    <  7:58635/149‑1 (MQ=255)                                                                                                                                            cacaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAg   >  7:93112/1‑149 (MQ=255)                                                                                                                                             acaGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGc  <  8:93112/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   CTGGCTGGGGATTCATGCCGATGCCTACTCGCTGTCGGCGCAAGCGGTGGCCTTGCTGGCCATCGTGTTCCTTTACGGCCGCTCGTGGCTGGCCGAAAAGCGCAGTGCTGCGGCCAACTGACTGGAGCCTGCGCGGTCCCACAGGGACCGCGCAAGGTCTCAGATCTTGCGCAAGCCTGATGCTCAGGCTGGCTGTGGCGCCTTGAGAATGTTCTGTAGCAATTGGACACTTTCAACGGCGTCATGACCCAGGCTGGTCAGGTAGCCGCCATCCGGCTGATCGGTCAGC  >  NC_002947/233692‑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: ‑