Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,632,137	T→A	11.7%	P283P (CCT→CCA) ‡	lepA →	elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,632,137	0	T	A	11.7%	196.7 / 12.6	67	P283P (CCT→CCA) ‡	lepA	elongation factor 4
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base A (5/3);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAA  >  NC_002947/1631998‑1632285                                                                                                                                            |                                                                                                                                                      tCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAgg                                                                                                                                             <  2:507190/150‑1 (MQ=255)                  ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgt                                                                                                                                            >  1:1072537/1‑134 (MQ=255)                  ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGgt                                                                                                                                            <  2:1072537/134‑1 (MQ=255)                  ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                            >  1:939989/1‑150 (MQ=255)                  ttCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaag                                                                                                                            >  1:203693/1‑150 (MQ=255)                      ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                             <  1:874568/145‑1 (MQ=255)                      ccccGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCaa                                                                                                                             >  2:874568/1‑145 (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                                                                                                                                               >  2:698879/1‑103 (MQ=255)                                              aTCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCTCAGGGGTCGa                                                                                                                                               <  1:698879/103‑1 (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                                                                                    >  2:6293/1‑123 (MQ=255)                                                                                    aTCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTc                                                                                    <  1:6293/123‑1 (MQ=255)                                                                                      cAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGa                                                       <  2:294361/150‑1 (MQ=255)                                                                                      cAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGa                                                       <  2:345137/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                                                                                                                        <  1:357357/57‑1 (MQ=255)                                                                                                                   acacCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAagaaga                                                                                                                        >  2:357357/1‑57 (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:1293460/1‑150 (MQ=255)                                                                                                                                  ccTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACa           <  1:791672/150‑1 (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)                                                                                                                                           cTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGaa                                           >  1:68745/1‑109 (MQ=255)                                                                                                                                           cTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGaa                                           <  2:68745/109‑1 (MQ=255)                                                                                                                                           cTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGaaa  >  1:782500/1‑150 (MQ=255)                                                                                                                                            |                                                                                                                                                      TCGACAGCGTCGGTGTCTTCACCCCGAAGCACACCCAGACCGCTGATCTGAAAGCCGGCGAAGTAGGCTTCATCATCGCCAGCATCAAGGACATCCACGGTGCCCCGGTTGGTGACACCCTGACCCTGTCCTCGACCCCTGAGGTCGAGGTGTTGCCAGGCTTCAAGAAGATTCAGCCGCAG‑GTCTACGCCGGCCTGTTCCCGGTCAGCTCCGACGACTTCGAGGATTTCCGCGACGCTCTGCAGAAGCTGACCCTGAACGACTCCTCGCTGCAGTACATGCCGGAAA  >  NC_002947/1631998‑1632285

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