Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,256,883	A→C	7.6%	noncoding (432/451 nt)	PP_mr14 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,256,883	0	A	C	7.6%	234.1 / 10.2	79	noncoding (432/451 nt)	PP_mr14	cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (35/37);  new base C (2/4);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.66e-01

CCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTT  >  NC_002947/1256737‑1257031                                                                                                                                                   |                                                                                                                                                     ccATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGaa                                                                                                                                                     <  3:368811/148‑1 (MQ=255)       gCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGt                                                                                                                                              >  3:88571/1‑149 (MQ=255)       gCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGt                                                                                                                                              >  2:271840/1‑149 (MQ=255)        cGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGtt                                                                                                                                             <  6:76403/149‑1 (MQ=255)            gACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttt                                                                                                                                         <  8:284892/149‑1 (MQ=255)            gACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttt                                                                                                                                         <  2:57667/149‑1 (MQ=255)              cTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGttttttta                                                                                                                                       <  7:331774/149‑2 (MQ=255)                 aCCTGCCATCCCGAGTTTGCCCGCACACAGCGACTGCTCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                      <  8:97269/147‑3 (MQ=38)                    tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGc                                                                                                                                  >  7:82436/1‑148 (MQ=255)                      ccATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                <  6:33806/148‑1 (MQ=255)                        aTAGCGTGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGGAGTCCCTGTGGGAGCGGGAGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                <  2:125740/130‑1 (MQ=255)                        aTACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATAACGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                >  1:125740/1‑130 (MQ=255)                           ccGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                          >  6:197035/1‑149 (MQ=255)                            cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGCGGCCAGCGCCTATTTAGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                          >  8:144554/1‑148 (MQ=255)                            cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                          >  8:131019/1‑148 (MQ=255)                            cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGc                                                                                                                          >  3:140810/1‑148 (MQ=255)                            cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCt                                                                                                                         <  3:282948/149‑1 (MQ=255)                                     ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACt                                                                                                                >  8:314125/1‑149 (MQ=255)                                      cGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                               <  2:246187/149‑1 (MQ=255)                                      cGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACt                                                                                                                >  6:16128/1‑148 (MQ=255)                                      cGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACt                                                                                                                >  2:28070/1‑148 (MQ=255)                                      cGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACt                                                                                                                >  3:356541/1‑148 (MQ=255)                                                  cTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGa                                                                                                    >  7:240330/1‑148 (MQ=255)                                                   tGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                  >  3:94647/1‑149 (MQ=255)                                                     cgcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                  <  1:36352/147‑1 (MQ=255)                                                      gcCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtc                                                                                               <  7:314125/149‑1 (MQ=255)                                                       ccGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtct                                                                                              >  1:203062/1‑149 (MQ=255)                                                       ccGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtct                                                                                              >  5:63847/1‑149 (MQ=255)                                                         gTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACttcttc                                                                                            <  3:315680/149‑1 (MQ=255)                                                           aCTTGAAGCCTGTGCAGTCACTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAg                                                                                          >  5:194897/1‑149 (MQ=255)                                                            cTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGc                                                                                         <  7:131019/149‑1 (MQ=255)                                                             ttGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtc                                                                                               <  7:46503/142‑1 (MQ=255)                                                             ttGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtc                                                                                               >  8:46503/1‑142 (MQ=255)                                                                  gCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCAt                                                                                   >  5:158655/1‑149 (MQ=255)                                                                       tgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCtt                                                                              >  7:322101/1‑149 (MQ=255)                                                                         cAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTc                                                                            <  8:240330/149‑1 (MQ=255)                                                                          aGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCt                                                                           <  4:94647/149‑1 (MQ=255)                                                                           gTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                          <  3:158858/149‑1 (MQ=255)                                                                            tCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTg                                                                          <  2:203062/148‑1 (MQ=255)                                                                             cccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGc                                                                        <  6:327758/149‑1 (MQ=255)                                                                               cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtc                                                                                               >  7:320432/1‑124 (MQ=255)                                                                               cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtc                                                                                               <  8:320432/124‑1 (MQ=255)                                                                                  tgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCgg                                                                   >  5:337632/1‑149 (MQ=255)                                                                                  tgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCATTCTGGCCGCgg                                                                   <  8:322101/149‑1 (MQ=255)                                                                                     gAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGAATTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTa                                                                <  1:285305/149‑1 (MQ=255)                                                                                      aGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAg                                                               >  2:243938/1‑149 (MQ=255)                                                                                       gCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGc                                                              <  1:243938/149‑1 (MQ=255)                                                                                        cGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGc                                                              <  3:356916/148‑1 (MQ=255)                                                                                          ggCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgc                                                                     >  4:92190/1‑139 (MQ=255)                                                                                          ggCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgc                                                                     <  3:92190/139‑1 (MQ=255)                                                                                           gCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgac                                                           >  2:8258/1‑148 (MQ=255)                                                                                              aGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGaaa                                                       <  1:8258/149‑1 (MQ=255)                                                                                              aGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGaaa                                                       <  1:255090/149‑1 (MQ=255)                                                                                               gCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAg                                                      >  7:68576/1‑149 (MQ=255)                                                                                                cccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCttt                                                                             <  7:254168/125‑1 (MQ=255)                                                                                                cccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCttt                                                                             >  8:254168/1‑125 (MQ=255)                                                                                                  cgcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                   <  4:356541/149‑1 (MQ=255)                                                                                                  cgcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                   <  1:28070/149‑1 (MQ=255)                                                                                                  cgcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                   <  5:16128/149‑1 (MQ=255)                                                                                                      aaGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGc                                               <  8:82436/149‑1 (MQ=255)                                                                                                         agGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTg                                            <  4:305670/149‑1 (MQ=255)                                                                                                               gcgcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCcttc                                      >  7:311764/1‑149 (MQ=255)                                                                                                                cgcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCcttct                                     <  4:140810/149‑1 (MQ=255)                                                                                                                  ccTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCttt                                   <  7:284012/149‑1 (MQ=255)                                                                                                                     aTTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAg                                >  1:124911/1‑149 (MQ=255)                                                                                                                           aGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGtttt                          >  2:185770/1‑149 (MQ=255)                                                                                                                             cAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGa                        <  6:337632/149‑1 (MQ=255)                                                                                                                                gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGc                     >  2:343302/1‑149 (MQ=255)                                                                                                                                 cATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                    <  6:63847/149‑1 (MQ=255)                                                                                                                                  aTAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCa                                                                                           <  5:175170/77‑1 (MQ=255)                                                                                                                                  aTAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCa                                                                                           >  6:175170/1‑77 (MQ=255)                                                                                                                                  aTAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCggg                                                                  >  4:127423/1‑102 (MQ=255)                                                                                                                                  aTAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCggg                                                                  <  3:127423/102‑1 (MQ=255)                                                                                                                                    aaaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGt                  >  3:68521/1‑148 (MQ=255)                                                                                                                                     aaaaaaCCCGCCTTCTGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGt                  >  3:194179/1‑147 (MQ=255)                                                                                                                                     aaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTTGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACTAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTc                 <  6:158655/148‑1 (MQ=255)                                                                                                                                     aaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTc                 <  4:359351/148‑1 (MQ=255)                                                                                                                                      aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAg               >  8:315903/1‑149 (MQ=255)                                                                                                                                             cGCCAGTAGGCGTGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAATCCGTCAGCTTGGAt        <  2:85866/149‑1 (MQ=255)                                                                                                                                              gCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATg       >  7:192131/1‑149 (MQ=255)                                                                                                                                                   aaGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCttt  >  1:148546/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTT  >  NC_002947/1256737‑1257031

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