Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,632,135	C→T	12.7%	P283S (CCT→TCT) ‡	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,632,135	0	C	T	12.7%	158.3 / 13.0	63	P283S (CCT→TCT) ‡	lepA	elongation factor 4
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base T (5/3);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 7.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GGTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGA  >  NC_002947/1631987‑1632283                                                                                                                                                     |                                                                                                                                                      ggTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGAcccc                                                                                                                                                      >  2:784555/1‑150 (MQ=255)  gTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGAccctg                                                                                                                                                     <  1:400942/150‑3 (MQ=255)            tCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAgg                                                                                                                                           <  2:507190/150‑1 (MQ=255)                             ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgt                                                                                                                                          >  1:1072537/1‑134 (MQ=255)                             ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgt                                                                                                                                          <  2:1072537/134‑1 (MQ=255)                             ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  1:203693/1‑150 (MQ=255)                             ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                          >  1:939989/1‑150 (MQ=255)                                 ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                           >  2:874568/1‑145 (MQ=255)                                 ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                           <  1:874568/145‑1 (MQ=255)                                  cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAg                                                                                                                                            >  2:503515/1‑127 (MQ=255)                                  cccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGAg                                                                                                                                            <  1:503515/127‑1 (MQ=255)                                             cccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAg‑g                                                                                                         >  1:294361/1‑150 (MQ=255)                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑Gt                                                                                                        >  2:483397/1‑150 (MQ=255)                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACACCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGc                                                                                                                <  2:1170636/143‑1 (MQ=255)                                              ccAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACACCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGc                                                                                                                >  1:1170636/1‑143 (MQ=255)                                                         aTCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGa                                                                                                                                             <  1:698879/103‑1 (MQ=255)                                                         aTCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGa                                                                                                                                             >  2:698879/1‑103 (MQ=255)                                                               aaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGc                                                                                              >  1:429457/1‑143 (MQ=255)                                                               aaGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGc                                                                                              <  2:429457/143‑1 (MQ=255)                                                                 gCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCg                                                                                     >  2:1180/1‑150 (MQ=255)                                                                    ggCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTc                                                                                  >  2:1201135/1‑150 (MQ=255)                                                                          gTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   <  1:1417684/80‑1 (MQ=255)                                                                          gTAGGCTTCATAATAGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGa                                                                                                                                                   >  2:1417684/1‑80 (MQ=255)                                                                           tAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCcg                                                                           >  1:458491/1‑150 (MQ=255)                                                                            aGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCcga                                                                          <  1:1180274/150‑1 (MQ=255)                                                                              gCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTa                                                                                                     >  1:493053/1‑121 (MQ=255)                                                                              gCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTa                                                                                                     <  2:493053/121‑1 (MQ=255)                                                                               cTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                   >  1:503759/1‑91 (MQ=255)                                                                               cTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCa                                                                                                                                   <  2:503759/91‑1 (MQ=255)                                                                                    tcatcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCg                                                                  >  1:1209876/1‑150 (MQ=255)                                                                                     catcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 <  2:1209876/150‑1 (MQ=255)                                                                                     catcGCCAGCATCAAGGACATCCACGGTGCCCCGGTATCTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 <  1:784555/150‑1 (MQ=255)                                                                                       tcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCa                                                                                 >  1:1146520/1‑132 (MQ=255)                                                                                       tcGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCa                                                                                 <  2:1146520/132‑1 (MQ=255)                                                                                               agcAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCcgc                                                       <  2:1223572/148‑1 (MQ=255)                                                                                               aTCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTc                                                                                  <  1:6293/123‑1 (MQ=255)                                                                                               aTCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTc                                                                                  >  2:6293/1‑123 (MQ=255)                                                                                                 cAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGa                                                     <  2:345137/150‑1 (MQ=255)                                                                                                 cAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGa                                                     <  2:294361/150‑1 (MQ=255)                                                                                                  aaGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGAc                                                    <  2:203693/150‑1 (MQ=255)                                                                                                    ggACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGc                                                  >  1:176693/1‑150 (MQ=255)                                                                                                        atacACGGTGCGCGGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTCCCCCGACGCTCTg                                              >  2:707915/4‑150 (MQ=255)                                                                                                        aTCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 >  1:25174/1‑131 (MQ=255)                                                                                                        aTCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGa                                                                 <  2:25174/131‑1 (MQ=255)                                                                                                           cACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCa                                                                                                            >  1:1089611/1‑86 (MQ=255)                                                                                                           cACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCa                                                                                                            <  2:1089611/86‑1 (MQ=255)                                                                                                           cACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGCCGACATAGCGGCTTTCAGCGACTCGCTGCAg                                           >  2:791672/1‑150 (MQ=255)                                                                                                                  ccccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGTTCGAGGTGTTGACCAGTTGCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGa                                    <  1:728706/150‑1 (MQ=255)                                                                                                                    ccGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGAcc                                  >  2:92205/1‑150 (MQ=255)                                                                                                                            tGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGa                          <  2:939989/150‑1 (MQ=255)                                                                                                                              acacCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                      >  2:357357/1‑57 (MQ=255)                                                                                                                              acacCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                      <  1:357357/57‑1 (MQ=255)                                                                                                                              acacCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGAct                        >  1:947071/1‑150 (MQ=255)                                                                                                                               cacCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActc                       <  1:483397/150‑1 (MQ=255)                                                                                                                                 ccctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGc                                             >  1:583760/1‑126 (MQ=255)                                                                                                                                 ccctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGc                                             <  2:583760/126‑1 (MQ=255)                                                                                                                                  cctgaccctgTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGActcctc                    >  1:589054/1‑150 (MQ=255)                                                                                                                                         aacaccTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTc                                                                             <  2:854455/82‑1 (MQ=38)                                                                                                                                         aacaccTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTc                                                                             >  1:854455/5‑86 (MQ=38)                                                                                                                                         aacaccTCGACCTCAGGGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAg             >  1:364315/5‑150 (MQ=255)                                                                                                                                             ccTCGACCTCAGGGGTCGAGGAGTTGAAAGGCTTCAAGAAGATTCAGCC‑CAGTGACTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa         >  1:179771/1‑150 (MQ=255)                                                                                                                                             ccTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa         <  1:791672/150‑1 (MQ=255)                                                                                                                                             ccTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa         >  1:1293460/1‑150 (MQ=255)                                                                                                                                               tCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATg       >  1:705468/1‑150 (MQ=255)                                                                                                                                               tCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATg       >  1:252669/1‑150 (MQ=255)                                                                                                                                                    cccTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGa  >  2:290187/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GGTGCACTTGGTCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGA  >  NC_002947/1631987‑1632283

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