Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,176,698	G→A	7.8%	R460H (CGT→CAT) ‡	PP_3676 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,176,698	0	G	A	7.8%	180.9 / 3.9	64	R460H (CGT→CAT) ‡	PP_3676	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.65e-01

AAATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGT  >  NC_002947/4176551‑4176836                                                                                                                                                    |                                                                                                                                           aaaTACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCg                                                                                                                                            <  6:3713/148‑1 (MQ=255)     aCCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAAATATTTCTAACATCGACATCAAAGGGAACGGCGAGTTTCTCAAATTTCATAGGTAAACTTAATTCTACCAAGGGTGTTACGCTACATCAAATTTACCGTGAt                                                                                                                                        >  5:296316/1‑148 (MQ=255)          gagCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAAttt                                                                                                                                  <  4:232839/149‑1 (MQ=255)            gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                <  5:226180/149‑1 (MQ=255)            gcgcGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                <  7:121968/149‑1 (MQ=255)                ggCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            <  8:299149/149‑5 (MQ=255)                    aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            <  6:242601/145‑5 (MQ=255)                    aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            >  5:242601/1‑141 (MQ=255)                    aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            >  1:153427/1‑141 (MQ=255)                    aCTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTATCACGGTAAATTTGatgta                                                                                                                            <  2:153427/145‑5 (MQ=255)                      tGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTgg                                                                                                                       >  4:55313/1‑148 (MQ=255)                        cAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCGAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCAt                                                                                                                    >  7:69774/1‑149 (MQ=255)                           gACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCAc                                                                                                                 >  5:39312/1‑149 (MQ=255)                             ccGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCAccc                                                                                                               <  1:228730/149‑1 (MQ=255)                                  actCAGGATGGCGTAAATTTTTTCTATCTTCGACATCAACGGGAACGGCGCGTTTCTCATATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATATATTTGAGGCCGCTGGCATCACCCGCGGa                                                                                                          <  2:89221/147‑1 (MQ=255)                                  aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                <  1:286170/127‑1 (MQ=255)                                  aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGa                                                                                                                                >  2:286170/1‑127 (MQ=255)                                  aaTCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGa                                                                                                          <  2:54857/149‑1 (MQ=255)                                         aTGGCGTAAATTATTTCTAACATCGACATCAAAGGGTACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCa                                                                                                   <  6:296316/149‑1 (MQ=255)                                           ggCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCAc                                                                                                  >  3:44001/1‑148 (MQ=255)                                             cGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                               <  7:241739/149‑1 (MQ=255)                                                         cTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTc                                                                                    >  4:110895/1‑148 (MQ=255)                                                           aaCATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTTATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGAt                                                                                 <  6:39312/149‑1 (MQ=255)                                                               tCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGca                                                                              >  8:169400/1‑148 (MQ=255)                                                                  aCATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcacaca                                                                          <  7:169400/149‑1 (MQ=255)                                                                     tCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAggg                                                                       >  2:262141/1‑149 (MQ=255)                                                                       aaaGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGAt                                                                                 <  2:319555/137‑1 (MQ=255)                                                                       aaaGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGAt                                                                                 >  1:319555/1‑137 (MQ=255)                                                                         aGGGAACGGCGCGTTTCTCAAATTTCATAGTTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCAtt                                                                   <  4:44001/149‑1 (MQ=255)                                                                         aGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                               >  4:170065/1‑121 (MQ=255)                                                                         aGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCt                                                                                               <  3:170065/121‑1 (MQ=255)                                                                           ggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                  >  5:244385/1‑148 (MQ=255)                                                                           ggAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTc                                                                  >  1:255111/1‑148 (MQ=255)                                                                                gtcgcgTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAt                                                            <  6:281903/147‑1 (MQ=255)                                                                                   gcgTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAtttt                                                         >  5:165724/1‑149 (MQ=255)                                                                                         ctaaaaTTTCATAGGTAAACTTAATTCTTAAAAGGGTGTTAAGCTACATCAAATTTAACGTGATAAATTTTAGGCCGCTGGCATTAACCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTATAGAATTTTTACgg                                                    >  2:67402/4‑148 (MQ=255)                                                                                          tcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGaa                                                  >  3:73351/1‑149 (MQ=255)                                                                                          tcAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGaa                                                  >  5:162816/1‑149 (MQ=255)                                                                                               tttCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAt                                                            <  4:235965/134‑1 (MQ=255)                                                                                               tttCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAAt                                                            >  3:235965/1‑134 (MQ=255)                                                                                                     aGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCATACAGGGCATTCTAGAATTTTTACGGAACGGCAGCAtt                                        >  2:28362/1‑148 (MQ=255)                                                                                                      ggTAAACTTAATTCTTCGAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTc                                       <  8:69774/148‑1 (MQ=255)                                                                                                       gTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAActg                                               >  8:14383/1‑137 (MQ=255)                                                                                                       gTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACgg                                               <  7:14383/139‑1 (MQ=255)                                                                                                          aaCTTAATTCTTCCAATGTTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATTACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGTGCATTCTATAATTTTTACGGAACGGCAGCATTCAAaga                                  <  1:67402/149‑1 (MQ=255)                                                                                                          aaCTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCATACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAag                                   <  1:28362/148‑1 (MQ=255)                                                                                                          aaCTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCaaa                                    >  6:193315/1‑147 (MQ=255)                                                                                                                    ttcCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAAc                                                 <  3:156720/124‑1 (MQ=255)                                                                                                                    ttcCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAAc                                                 >  4:156720/1‑124 (MQ=255)                                                                                                                      ccAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGAc                         <  5:193315/146‑1 (MQ=255)                                                                                                                       cAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcacac                                                                           <  6:294885/95‑1 (MQ=255)                                                                                                                       cAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGcacac                                                                           >  5:294885/1‑95 (MQ=255)                                                                                                                        aaGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACt                    >  4:239693/1‑149 (MQ=255)                                                                                                                         aGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTg                   <  3:239693/149‑1 (MQ=255)                                                                                                                              gtTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGt                                                                                        <  3:266853/75‑1 (MQ=255)                                                                                                                              gtTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGt                                                                                        >  4:266853/1‑75 (MQ=255)                                                                                                                                    cTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATg        <  3:208041/149‑1 (MQ=255)                                                                                                                                        aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGcc      >  7:228645/1‑147 (MQ=255)                                                                                                                                        aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCaaa    >  4:207296/1‑149 (MQ=255)                                                                                                                                        aTCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCaa    <  3:110895/149‑1 (MQ=255)                                                                                                                                         tCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATTCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAg   <  6:244385/149‑1 (MQ=255)                                                                                                                                         tCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAg   >  6:114597/1‑149 (MQ=255)                                                                                                                                          cAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGt  >  6:125948/1‑149 (MQ=255)                                                                                                                                             aTTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCTCGGAGTGCTCACCCTATTCCGTATTCGATTCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGc       <  1:52472/141‑1 (MQ=255)                                                                                                                                             aTTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGc       >  2:52472/1‑141 (MQ=255)                                                                                                                                                    |                                                                                                                                           AAATACCAAGAGCGCGGCTACTGCAAGACCGAAAATCAGGATGGCGTAAATTATTTCTAACATCGACATCAAAGGGAACGGCGCGTTTCTCAAATTTCATAGGTAAACTTAATTCTTCCAAGGGTGTTACGCTACATCAAATTTACCGTGATAAATTTGAGGCCGCTGGCATCACCCGCGGAGTGCTCACCCTATTCCGTATTCGATGCACACAGGGCATTCTAGAATTTTTACGGAACGGCAGCATTCAAAGAGTAGCTGACTTACTGGGAAATACATGCCAAGT  >  NC_002947/4176551‑4176836

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