Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,040,794	T→A	10.9%	intergenic (+19/+15)	PP_5292 → / ← PP_5293	catabolite repression control protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,040,794	0	T	A	10.9%	173.7 / 10.6	64	intergenic (+19/+15)	PP_5292/PP_5293	catabolite repression control protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/30);  new base A (3/4);  total (30/34)
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.66e-01

GTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAA  >  NC_002947/6040649‑6040940                                                                                                                                                   |                                                                                                                                                     gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCAtt                                                                                                                                                     >  2:36757/1‑147 (MQ=255) gTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCAtt                                                                                                                                                     <  1:36757/147‑1 (MQ=255)     gACTACCAGATCCTCACCCCGGGGCTGCGCCGCTTCGTGCGCAACGCTAGCCTGCCGCGTCAGCCACGCTTCACCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTTCATTGCAGc                                                                                                                                                <  6:315613/148‑1 (MQ=255)     gACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTATAAGGCCATTGGGGCTGCAATGCAGcc                                                                                                                                               <  7:260077/149‑1 (MQ=255)      aCTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGccc                                                                                                                                              <  1:24256/149‑1 (MQ=255)      aCTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc                                                                                                                                              <  1:1336/149‑1 (MQ=255)      aCTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGccc                                                                                                                                              >  7:274532/1‑149 (MQ=255)        tACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGcccc                                                                                                                                             >  6:10815/1‑148 (MQ=255)           cAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATg                                                                                                                                         <  5:10815/149‑1 (MQ=255)           cAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATg                                                                                                                                         <  8:274532/149‑1 (MQ=255)                 cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCa                                                                                                                                            >  1:343329/1‑140 (MQ=255)                 cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCa                                                                                                                                            <  2:343329/140‑1 (MQ=255)                 cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTAtt                                                                                                                                   <  3:56040/149‑1 (MQ=255)                 cTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCAATGCAGCCCCAATggccttt                                                                                                                                   <  3:174176/149‑6 (MQ=255)                  tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      >  5:180399/1‑145 (MQ=255)                  tCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      <  6:180399/145‑1 (MQ=255)                       ccGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCg                                                                                                                             <  1:279402/149‑1 (MQ=255)                       ccGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCg                                                                                                                             <  5:226071/149‑1 (MQ=255)                                             gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCg                                                                                                             <  7:79428/142‑1 (MQ=255)                                             gcAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGAAGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCg                                                                                                             >  8:79428/1‑142 (MQ=255)                                              cAACGCTCGCCTGCCGCGTCAGCGACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          >  1:344544/1‑113 (MQ=255)                                              cAACGCTCGCCTGCCGCGTCAGCGACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          <  2:344544/113‑1 (MQ=255)                                              cAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCCGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGa                                                                                                     >  6:262873/1‑149 (MQ=255)                                                 cGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATa‑a                                                                                                 <  7:249279/149‑1 (MQ=255)                                                      gCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGc                                                                                                       <  2:23021/139‑1 (MQ=255)                                                      gCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGc                                                                                                       >  1:23021/1‑139 (MQ=255)                                                      gCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCtt                                                                                            >  5:44261/1‑149 (MQ=255)                                                        cTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGgc‑cgc                                                                                                                       >  2:103038/1‑121 (MQ=255)                                                        cTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGgc‑cgc                                                                                                                       <  1:103038/121‑1 (MQ=255)                                                                  aGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCGTTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTc                                                                                      <  2:294576/143‑1 (MQ=255)                                                                  aGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCGTTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTc                                                                                      >  1:294576/1‑143 (MQ=255)                                                                        gCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCAc                                                                          >  2:135349/1‑149 (MQ=255)                                                                         cTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCgg                                                                                                            >  8:193578/1‑115 (MQ=255)                                                                         cTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCgg                                                                                                            <  7:193578/115‑1 (MQ=255)                                                                          ttCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGCGTGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGatcatc                                                                                                              >  1:68018/1‑113 (MQ=255)                                                                          ttCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGatcatc                                                                                                              <  2:68018/112‑1 (MQ=255)                                                                            ctcCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCAccgcc                                                                      >  2:62338/1‑149 (MQ=255)                                                                               ccAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACcgccgc                                                                    >  4:258761/1‑148 (MQ=255)                                                                               ccAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACcgccgc                                                                    <  6:44261/148‑1 (MQ=255)                                                                                        aCCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATc                                                          <  3:61954/149‑1 (MQ=255)                                                                                           gCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAat                                                       <  7:224792/149‑1 (MQ=255)                                                                                               aTCGTCGACTATGACTGGACGC‑TGACCATGTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATa‑a                                                                                                 >  2:28579/1‑103 (MQ=255)                                                                                               aTCGTCGACTATGACTGGACGC‑TGACCATGTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTCATCGGC‑CGCCAGATCATCGGGTAGCGATa‑a                                                                                                 <  1:28579/103‑1 (MQ=255)                                                                                                  gtcgACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATatcatc                                                 <  1:135349/148‑1 (MQ=255)                                                                                                        tATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCa                                          >  4:139908/1‑149 (MQ=255)                                                                                                            aCTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGt                                      <  1:62338/149‑1 (MQ=255)                                                                                                              tGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGt                                      >  1:1698/1‑147 (MQ=255)                                                                                                                gACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAg                                   >  2:99648/1‑148 (MQ=255)                                                                                                                gACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGc                                  <  3:258761/149‑1 (MQ=255)                                                                                                                   gC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCgatgatga                                                            >  8:254280/1‑120 (MQ=255)                                                                                                                   gC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCgatgatga                                                            <  7:254280/120‑1 (MQ=255)                                                                                                                            tCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCaa                                                        >  1:110440/1‑116 (MQ=255)                                                                                                                            tCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCaa                                                        <  2:110440/116‑1 (MQ=255)                                                                                                                             cTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCa                                          <  6:211394/129‑1 (MQ=255)                                                                                                                             cTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCa                                          >  5:211394/1‑129 (MQ=255)                                                                                                                                aataaccATTGGGGCTGCAATGCAGCCCCAATGGTTATTTGATCGGCACGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCAc                    >  1:252649/6‑149 (MQ=255)                                                                                                                                   ggCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCaaa                          >  4:100620/1‑139 (MQ=255)                                                                                                                                   ggCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCaaa                          <  3:100620/139‑1 (MQ=255)                                                                                                                                    gCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATATGGGCTTGCCCTCATTCGCCCGCACCGCCGCgatgat                                                              >  1:104780/1‑103 (MQ=255)                                                                                                                                      cATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCg             <  1:99648/149‑1 (MQ=255)                                                                                                                                        ttGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGatc                                                                                                                 <  6:281793/48‑1 (MQ=255)                                                                                                                                        ttGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGatc                                                                                                                 >  5:281793/1‑48 (MQ=255)                                                                                                                                           gggCTGCAATGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATaaa        >  4:280192/1‑149 (MQ=255)                                                                                                                                                gCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATAATCGGGTAGCGATA‑AGCCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTCAGCATGGCAAACAGCACCAGCACGATAAACACg    >  3:61718/1‑148 (MQ=255)                                                                                                                                                 cATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGaa  <  3:139908/149‑1 (MQ=255)                                                                                                                                                  aTTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGaa  >  3:123541/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GTTCGACTACCAGATCCTCACCCCGGGCCTGCGCCGCTTCGTGCGCAACGCTCGCCTGCCGCGTCAGCCACGCTTCTCCCAGCATGCACCGCTGATCGTCGACTATGACTGGACGC‑TGACCATCTAAAAGGCCATTGGGGCTGCATTGCAGCCCCAATGGTTATTTGATCGGC‑CGCCAGATCATCGGGTAGCGATA‑AGGCTTGCCCTCATTCGCCCGCACCGCCGCGATGATGATCAATATCATCACTGCCACCGTGAGCATGGCAAACAGCACCAGCCCGATAAACACGAA  >  NC_002947/6040649‑6040940

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: ‑