Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,123,420	A→C	6.6%	noncoding (381/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,420	0	A	C	6.6%	189.0 / 4.5	61	noncoding (381/398 nt)	PP_mr13	cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base C (2/2);  total (30/31)
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.44e-01

CACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAA  >  NC_002947/1123273‑1123567                                                                                                                                                    |                                                                                                                                                    cacaGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCa                                                                                                                                                     <  5:323532/148‑1 (MQ=255)  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                                                                                                                                 <  8:202046/149‑1 (MQ=255)                    aCATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                 <  4:244352/149‑1 (MQ=255)                    aCATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCaaa                                                                                                                                 <  4:242060/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                                                                                                                  >  1:141577/1‑133 (MQ=255)                                                   cGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATc                                                                                                                  <  2:141577/133‑1 (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                    >  5:271673/1‑145 (MQ=255)                                                                                                                                     aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    >  1:161719/1‑145 (MQ=255)                                                                                                                                     aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  2:161719/145‑1 (MQ=255)                                                                                                                                     aaaaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  6:271673/145‑1 (MQ=255)                                                                                                                                      aaaaaaGCCCCTGCACTTGGCAGGGGCTTTTGAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCtt               >  1:62797/1‑149 (MQ=255)                                                                                                                                       aaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              <  7:214454/149‑1 (MQ=255)                                                                                                                                       aaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              <  8:15918/149‑1 (MQ=255)                                                                                                                                       aaaaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              <  3:79650/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)                                                                                                                                                    |                                                                                                                                                    CACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAA  >  NC_002947/1123273‑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: ‑