Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,123,421	C→A	6.4%	noncoding (380/398 nt)	PP_mr13 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,123,421	0	C	A	6.4%	197.5 / 3.7	63	noncoding (380/398 nt)	PP_mr13	cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base A (2/2);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-01

ACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAA  >  NC_002947/1123274‑1123567                                                                                                                                                    |                                                                                                                                                   acaGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGACCCTGCACt                                                                                                                                                   >  7:175253/1‑149 (MQ=255)  aaGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACtt                                                                                                                                                  <  5:301933/148‑1 (MQ=255)       aaaaGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTTGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAg                                                                                                                                             <  5:241116/149‑1 (MQ=255)                aGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtttt                                                                                                                                     >  6:303331/1‑148 (MQ=255)                   aCATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                 <  4:244352/149‑1 (MQ=255)                   aCATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                 <  4:242060/149‑1 (MQ=255)                   aCATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                 <  8:202046/149‑1 (MQ=255)                   aCATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCa                                                                                                                                   >  8:214454/1‑147 (MQ=255)                     aTCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCCAAGTGCAGGGGCttttttttgc                                                                                                                               <  2:118882/149‑7 (MQ=39)                        gCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCTGGGGCTTTTCAAAGTCTg                                                                                                                            <  6:70973/149‑1 (MQ=255)                         cTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTgg                                                                                                                           <  3:66512/149‑1 (MQ=255)                          tGCCAGTCGCCGGCCAAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGATTGGACATGCAGGAAGCTCCAAAAGAAGAACCCACAATTCTACGTGCAAAAAAAAGCCACTTAACTTAGCAGGGGCTTTTCAAAGTCTg                                                                                                                            >  4:66643/1‑147 (MQ=255)                            ccAGGCGCCGGCCCAGGACGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGCCTTGGACATTCAGGACCCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGgcgc                                                                                                                        <  4:315833/149‑1 (MQ=255)                                                  cGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATc                                                                                                                  <  2:141577/133‑1 (MQ=255)                                                  cGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATc                                                                                                                  >  1:141577/1‑133 (MQ=255)                                                  cGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCtt                                                                                                  >  5:238028/1‑149 (MQ=255)                                                    gTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTg                                                                                                 >  5:253172/1‑148 (MQ=255)                                                        aaaaTCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGaa                                                                                                      >  7:65719/1‑139 (MQ=255)                                                        aaaaTCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGaa                                                                                                      <  8:65719/139‑1 (MQ=255)                                                                ccTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGc                                                                                    <  8:134782/149‑1 (MQ=255)                                                                         tgagtAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCt                                                                                                                             <  7:314453/99‑1 (MQ=255)                                                                         tgagtAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCt                                                                                                                             >  8:314453/1‑99 (MQ=255)                                                                         tgagtAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTAGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGc                                                                           <  5:303331/149‑1 (MQ=255)                                                                           agtAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCt                                                                         >  7:15918/1‑149 (MQ=255)                                                                                  gCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTg                                                                  >  6:259137/1‑149 (MQ=255)                                                                                  gCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTg                                                                  >  5:129895/1‑149 (MQ=255)                                                                                      gtaCATGCAGGACGCTCAAATAGACGAGCCCACTATTCTACGGTCAAAAAAATGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGTCGCTGGATCAGACAGGTTGAACCTTTTCATCCTGCATGCCTTTCTGGCCTTGGACTGCtt                                                               <  3:66643/146‑1 (MQ=255)                                                                                         cgatctggACGCTCCAAAAGACGATCCCAATATTCTACTGTCAAAAAAAAGCCCCTGCACTTGTCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGaa                                                           <  3:259529/143‑1 (MQ=255)                                                                                         cATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGaa                                                           <  6:253172/149‑1 (MQ=255)                                                                                         cATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGaa                                                           <  6:129895/149‑1 (MQ=255)                                                                                                gACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACt                                                    >  7:171665/1‑149 (MQ=255)                                                                                                caCGCTCAAAAAGACGAGCCCACTATTATAAGGGCAAAAAAAAAACCCTCCACTTGCAACGGGCTTTTAAAAGACTGGAGCTACATCAGACAGGTTTAAACTTCACACACTCCAACCCTTTCTGGCATTGGACTGCTTCGAAgagatc                                                     >  4:259529/2‑143 (MQ=255)                                                                                                     tCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTg                                               >  1:216597/1‑149 (MQ=255)                                                                                                      ccAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTgg                                              >  3:178404/1‑149 (MQ=255)                                                                                                       cAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGc                                             >  6:250084/1‑149 (MQ=255)                                                                                                       cAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGc                                             >  2:331969/1‑149 (MQ=255)                                                                                                        aaaaGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGcc                                            >  4:79650/1‑149 (MQ=255)                                                                                                                  cccACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAgg                                  <  5:250084/149‑1 (MQ=255)                                                                                                                  cccACTATACTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAgg                                  >  1:135380/1‑149 (MQ=255)                                                                                                                            tACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGcc                                                                          >  4:281139/1‑99 (MQ=255)                                                                                                                            tACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGcc                                                                          <  3:281139/99‑1 (MQ=255)                                                                                                                               gggCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGt                     >  4:268796/1‑149 (MQ=255)                                                                                                                                ggCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  1:331969/149‑1 (MQ=255)                                                                                                                                 gCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTg                                                                       >  8:329547/1‑97 (MQ=255)                                                                                                                                 gCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTg                                                                       <  7:329547/97‑1 (MQ=255)                                                                                                                                   tttgaaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  4:349670/142‑1 (MQ=39)                                                                                                                                   tttgaaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    >  3:349670/5‑146 (MQ=39)                                                                                                                                    ttgaaaaGCCCCTGCCAAGTGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCt                >  8:247961/4‑149 (MQ=39)                                                                                                                                    aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    >  1:161719/1‑145 (MQ=255)                                                                                                                                    aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  2:161719/145‑1 (MQ=255)                                                                                                                                    aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  6:271673/145‑1 (MQ=255)                                                                                                                                    aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    >  5:271673/1‑145 (MQ=255)                                                                                                                                     aaaaaaGCCCCTGCACTTGGCAGGGGCTTTTGAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCtt               >  1:62797/1‑149 (MQ=255)                                                                                                                                      aaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              <  7:214454/149‑1 (MQ=255)                                                                                                                                      aaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              <  3:79650/149‑1 (MQ=255)                                                                                                                                      aaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              <  8:15918/149‑1 (MQ=255)                                                                                                                                             cccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAgcg       <  6:238028/149‑1 (MQ=255)                                                                                                                                              ccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTTCCTTCGATAgcgc      <  2:135380/149‑1 (MQ=255)                                                                                                                                               cTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCAGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGcgcg     >  8:76536/1‑149 (MQ=255)                                                                                                                                                 gCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCgg    >  3:68004/1‑148 (MQ=255)                                                                                                                                                  cACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGaa  <  4:178404/149‑1 (MQ=255)                                                                                                                                                   aCTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTc                                     >  7:337530/1‑113 (MQ=255)                                                                                                                                                   aCTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTc                                     <  8:337530/113‑1 (MQ=255)                                                                                                                                                    cTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGctc                               >  7:82257/1‑118 (MQ=255)                                                                                                                                                    cTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGctc                               <  8:82257/118‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   ACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAA  >  NC_002947/1123274‑1123567

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