Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,123,421	C→A	6.5%	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.5%	219.6 / 3.5	71	noncoding (380/398 nt)	PP_mr13	cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base C (27/39);  new base A (3/2);  total (30/41)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CACAGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGAA  >  NC_002947/1123273‑1123567                                                                                                                                                     |                                                                                                                                                   cacaGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCAc                                                                                                                                                    >  8:46639/1‑149 (MQ=255)  acaGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACt                                                                                                                                                   >  3:228726/1‑149 (MQ=255)    aGTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTg                                                                                                                                                 >  5:122577/1‑149 (MQ=255)     gTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACGGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGCACATGCAGGACGCTCCAAAAGACGAGCCGACTATTCTACGGGCAAAAAAAAGCCCCAGCACTAgg                                                                                                                                                >  5:247899/1‑149 (MQ=255)     gTAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTgg                                                                                                                                                <  5:266614/149‑1 (MQ=255)      tAAAAAGTCATAGAACATCGCTGCCAGGCGCCGGCCCACGCCGACCGGTGGAAAATCCCACTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAACCCCCTCCCCTTgc                                                                                                                                                >  6:61349/1‑147 (MQ=255)        aaaaGTCATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCa                                                                                                                                              <  3:23811/148‑1 (MQ=255)            gtactAGAACCTCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACCAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGc                                                                                                                                         <  6:213305/145‑1 (MQ=255)             tCATAGAACATCGCTGCCAGGCGCCGCCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGAAATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGCCAAAAAAAATCCCCTGCACTGGGCAGTGGCt                                                                                                                                        >  8:268461/1‑149 (MQ=255)              cATAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCt                                                                                                                                        <  6:23595/148‑1 (MQ=255)               aTAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCttt                                                                                                                                      <  6:122577/149‑1 (MQ=255)               aTAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCttt                                                                                                                                      >  1:379757/1‑149 (MQ=255)                tAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGATCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtttt                                                                                                                                     <  8:152027/149‑1 (MQ=255)                tAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtttt                                                                                                                                     >  6:69046/1‑149 (MQ=255)                tAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCttt                                                                                                                                      >  7:154140/1‑148 (MQ=255)                tAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCttt                                                                                                                                      <  8:154140/148‑1 (MQ=255)                tAGAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTACCAAGTGCAGGGGCtttt                                                                                                                                     >  4:117436/1‑149 (MQ=255)                  gAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCCAAGTGCAGGGGCttttt                                                                                                                                    <  1:257892/148‑2 (MQ=255)                  gAACATCGCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTc                                                                                                                                    <  7:46639/148‑1 (MQ=255)                         gCTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCt                                                                                                                             <  7:268461/148‑1 (MQ=255)                          cTGCCAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTgg                                                                                                                           <  1:261714/149‑1 (MQ=255)                             ccAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCGCCAAAAGACGAGCCCACTAATCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGgc                                                                                                                          <  6:297544/147‑1 (MQ=255)                              cAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtttt                                                                                                                                     <  7:174456/135‑1 (MQ=255)                              cAGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCtttt                                                                                                                                     >  8:174456/1‑135 (MQ=255)                               aGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACtt                                                                                                                                                  >  4:186759/1‑121 (MQ=255)                               aGGCGCCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACtt                                                                                                                                                  <  3:186759/121‑1 (MQ=255)                                   gcCGGCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGAt                                                                                                                   >  7:265745/1‑148 (MQ=255)                                       gCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCt                                                                                                                                        <  8:103921/123‑1 (MQ=255)                                       gCCCAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCt                                                                                                                                        >  7:103921/1‑123 (MQ=255)                                         ccAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGTCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAg                                                                                                            >  5:233922/1‑149 (MQ=255)                                         ccAGGCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAg                                                                                                            >  1:351976/1‑149 (MQ=255)                                            ggCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGtt                                                                                                         <  5:69046/149‑1 (MQ=255)                                             gCCGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGtt                                                                                                         >  1:70337/1‑148 (MQ=255)                                              ccGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGa                                                                                                       <  7:105339/149‑1 (MQ=255)                                               cGACCGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGaa                                                                                                      <  2:379757/149‑1 (MQ=255)                                                gtccGGTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCTAAAGCCGAGCCCACTATTCTACGGGCTAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAAc                                                                                                     <  5:61349/147‑1 (MQ=255)                                                     gTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCCAAGTGCAGGGGCttttttt                                                                                                                                  <  4:262582/115‑4 (MQ=38)                                                     gTGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCCAAGTGCAGGGGCttttttt                                                                                                                                  >  3:262582/1‑112 (MQ=38)                                                      tGGAAAATCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGt                                                                                                >  4:271126/1‑148 (MQ=255)                                                          aaaTCCCCCTGGGGAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGcc                                                                                           <  3:35637/149‑1 (MQ=255)                                                                     gggAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGc                                                                                    <  1:263403/145‑1 (MQ=255)                                                                     gggAGTGAGTAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGc                                                                                    >  2:263403/1‑145 (MQ=255)                                                                             gtAGAGGCTTGGAAATGCAGGACGCTCCAAAAGACGAGGCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCtt                                                                        <  6:247899/149‑1 (MQ=255)                                                                              tAGAGGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTg                                                                       <  2:178476/149‑1 (MQ=255)                                                                                gagGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTgg                                                                      <  5:22222/148‑1 (MQ=255)                                                                                gagGCTTGGACATGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGa                                                                     <  6:85930/149‑1 (MQ=255)                                                                                          cATGCAGGACGCTCCAAAAGACCCGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGaa                                                           <  4:228726/149‑1 (MQ=255)                                                                                           aTGCAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAg                                                          >  8:146576/1‑149 (MQ=255)                                                                                              cAGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGt                                                        <  8:265745/148‑1 (MQ=255)                                                                                               aGGACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGa                                                      <  3:271126/149‑1 (MQ=255)                                                                                                 gACGCTCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACt                                                    >  8:272464/1‑149 (MQ=255)                                                                                                      tCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTg                                               >  8:102343/1‑149 (MQ=255)                                                                                                      tCCAAAAGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTg                                               >  5:303096/1‑149 (MQ=255)                                                                                                         aaaaGACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGcc                                            >  6:135470/1‑149 (MQ=255)                                                                                                             gACGAGCCCACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCcttct                                        <  2:287386/149‑1 (MQ=255)                                                                                                                   cccACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCa                                    <  2:179732/147‑1 (MQ=255)                                                                                                                   cccACTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCa                                    >  1:179732/1‑147 (MQ=255)                                                                                                                      aCTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTg                                                                                         <  6:291137/91‑1 (MQ=255)                                                                                                                      aCTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGctc                               <  7:102343/149‑1 (MQ=255)                                                                                                                      aCTATTCTACGGGCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGctc                               <  7:272464/149‑1 (MQ=255)                                                                                                                        tATTCTACGGGCAAAAAAAAGCCCCTGAACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCATTTCTGGCCTTGGACTGCTTCGAAGTTGACTTGCTGGCCTTCTTTCAGGCTCtt                             >  2:230364/1‑149 (MQ=255)                                                                                                                                 ggCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  5:135470/149‑1 (MQ=255)                                                                                                                                 ggCAAAAAAAAGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGtt                    <  6:135958/149‑1 (MQ=255)                                                                                                                                    aaaaaaaaTCCCATGAAATTGGCAGGGGCTTTTAAAAGTATGGCGCAGGATCAGACAGCTTGAACCTTATCAGCCTCCATGACTATCTGGCCATCTACTCCTTCGAAGGTGACTTGCTGGACTTCTTTCAGGCTCTTGAACCCGTTGcc                 >  4:349356/1‑149 (MQ=255)                                                                                                                                         aaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGa            <  2:70337/149‑1 (MQ=255)                                                                                                                                         aaaGCCCCTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGa            <  6:233922/149‑1 (MQ=255)                                                                                                                                               ccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCtt               >  4:54281/1‑140 (MQ=255)                                                                                                                                               ccTGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCtt               <  3:54281/140‑1 (MQ=255)                                                                                                                                                 tGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              >  7:349211/1‑139 (MQ=255)                                                                                                                                                 tGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTc              <  8:349211/139‑1 (MQ=255)                                                                                                                                                 tGCACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGcgcg     >  8:267594/1‑148 (MQ=255)                                                                                                                                                   cACTTGGCAGGGGCTTTTCAAAGTCTGGCGCTGGATCAGACAGGTTGAACCTTGTCAGCCTGCATGCCTTTCTGGCCTTGGACTGCTTCGAAGGTGACTTGCTGGCCTTCTTTCAGGCTCTTGAAGCCGTTGCCTTCGATAGCGCGGaa  <  7:146576/149‑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: ‑