Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,016,329	A→T	5.3%	I345N (ATC→AAC)	PP_4421 ←	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,016,329	0	A	T	5.3%	236.9 / 2.7	76	I345N (ATC→AAC)	PP_4421	aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base T (2/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GGATCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTAGG  >  NC_002947/5016182‑5016476                                                                                                                                                    |                                                                                                                                                     ggATCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGAt                                                                                                                                                     >  3:165048/1‑149 (MQ=255)   aTCCGCTACGAACTCAATGGCTGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTATCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGt                                                                                                                                                   <  8:312854/149‑1 (MQ=255)    tCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATTGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTc                                                                                                                                                  <  4:19060/149‑1 (MQ=255)    tCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTc                                                                                                                                                  >  6:335974/1‑149 (MQ=255)       gCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCgcg                                                                                                                                               >  6:403415/1‑147 (MQ=255)            gAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAg                                                                                                                                          >  6:76072/1‑149 (MQ=255)              aCTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGAGAGGAGAAGGTCTCATGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGATTGTTCTCACCCTCGACGATGTCGAGCAAAGCg                                                                                                                                        >  8:139202/1‑149 (MQ=255)                  aaTGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGTCGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGttt                                                                                                                                     <  5:335974/148‑1 (MQ=255)                  aaTGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGttt                                                                                                                                     <  1:180037/148‑1 (MQ=255)                  aaTGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGtt                                                                                                                                      <  3:169630/147‑1 (MQ=255)                  aaTGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTtg                                                                                                                                    >  8:289068/1‑149 (MQ=255)                   aTGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTtgc                                                                                                                                   <  5:131095/149‑1 (MQ=255)                       cggcCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGcc                                                                                                                               <  6:313371/149‑1 (MQ=255)                            atcatGTCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCGGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGa                                                                                                                           >  3:147078/1‑148 (MQ=255)                             tcatGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGaca                                                                                                                         >  3:237074/1‑149 (MQ=255)                              catGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGacac                                                                                                                        >  5:192549/1‑149 (MQ=255)                                 gCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGa                                                                                                                     >  6:41240/1‑149 (MQ=255)                                   ccACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGa                                                                                                                     >  7:120601/1‑147 (MQ=255)                                     acacCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAgg                                                                                                                 <  5:76072/149‑1 (MQ=255)                                        ccacGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTg                                                                                                              <  1:316905/149‑1 (MQ=255)                                          acGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGc                                                                                                            <  8:120601/149‑1 (MQ=255)                                            gCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGccc                                                                                                          <  7:289068/149‑1 (MQ=255)                                               cTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGAGACACCGATAGGGTGGCCCGAg                                                                                                       >  4:355379/1‑149 (MQ=255)                                                ttCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAg                                                                                                       >  6:303223/1‑148 (MQ=255)                                                  cGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCAGAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGt                                                                                                      <  4:200109/147‑1 (MQ=255)                                                     cgacgaTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGgtg                                                                                                 <  6:192549/149‑1 (MQ=255)                                                     cgacgaTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAgg                                                                                                                                              <  3:374278/104‑2 (MQ=255)                                                     cgacgaTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAgg                                                                                                                                              >  7:90244/1‑103 (MQ=255)                                                     cgacgaTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAgg                                                                                                                                              <  8:90244/104‑2 (MQ=255)                                                     cgacgaTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACATCGTCGAgg                                                                                                                                              >  4:374278/1‑103 (MQ=255)                                                      gacgaTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGgtg                                                                                                 >  2:5028/1‑148 (MQ=255)                                                         gaTTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGc                                                                                             >  8:226717/1‑149 (MQ=255)                                                              gCAGCGAGGAGAAGGTCTCCAGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTg                                                                                         <  4:147078/148‑1 (MQ=255)                                                                  cGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGgaga                                                                                    <  1:342147/149‑1 (MQ=255)                                                                   gaggagAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGaa                                                                                   <  1:273476/149‑1 (MQ=255)                                                                        gaAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCAc                                                                              >  1:172331/1‑149 (MQ=255)                                                                        gaAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCAc                                                                              >  2:292997/1‑149 (MQ=255)                                                                            gTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTAGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AAc                                                                         >  4:167076/1‑149 (MQ=255)                                                                               tcCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGt                                                                                          <  8:113001/130‑1 (MQ=255)                                                                               tcCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGt                                                                                          >  7:113001/1‑130 (MQ=255)                                                                                ccTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGa                                                                     >  5:286953/1‑149 (MQ=255)                                                                                 cTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGAt                                                                    <  5:41240/149‑1 (MQ=255)                                                                                    cATCGACTTTTGGAAATACTGACCAACCTAGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAg                                                                 >  1:236747/1‑149 (MQ=255)                                                                                        gACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGc                                                                                <  7:254625/131‑1 (MQ=255)                                                                                        gACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGc                                                                                >  8:254625/1‑131 (MQ=255)                                                                                             ttGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAAGGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTtc                                                          <  2:324848/147‑1 (MQ=255)                                                                                                gAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAg                                                                                              <  8:273614/109‑1 (MQ=255)                                                                                                gAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAg                                                                                              >  7:273614/1‑109 (MQ=255)                                                                                                 aaaTACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCAt                                                    <  2:140853/149‑1 (MQ=255)                                                                                                       tGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCt                                               >  4:197983/1‑148 (MQ=255)                                                                                                           cAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATac                                          <  7:226717/149‑1 (MQ=255)                                                                                                            aaCCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATaca                                         >  6:351689/1‑149 (MQ=255)                                                                                                             aCCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCt                                               >  3:31252/1‑142 (MQ=255)                                                                                                             aCCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCt                                               <  4:31252/142‑1 (MQ=255)                                                                                                              ccTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTg                                                                                         >  6:198950/1‑100 (MQ=255)                                                                                                              ccTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTg                                                                                         <  5:198950/100‑1 (MQ=255)                                                                                                              ccTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACAcc                                       >  2:120170/1‑149 (MQ=255)                                                                                                                ttGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCt                                      >  1:32336/1‑148 (MQ=255)                                                                                                                  gttgaagTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCAGAGAGTTTGCTGACGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTtc                                    <  7:139202/142‑1 (MQ=255)                                                                                                                     gACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCa                                 >  8:167302/1‑148 (MQ=255)                                                                                                                     gACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGGTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCa                                 <  7:167302/148‑1 (MQ=255)                                                                                                                       cGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACca                              <  4:237074/149‑1 (MQ=255)                                                                                                                        gTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACcaa                             <  5:303223/149‑1 (MQ=255)                                                                                                                             cTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATc                        >  8:34270/1‑149 (MQ=255)                                                                                                                             cTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATc                        >  3:359836/1‑149 (MQ=255)                                                                                                                             cTGGAATGTTCTCACCCTCGACGAGTTCGAGCACAGCGTTTGCTGCCGCGACACCGATATGGTGGCCCGAGTAGGTGTGGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATc                        <  5:318776/149‑1 (MQ=255)                                                                                                                                gAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCAATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCg                       >  1:101163/1‑147 (MQ=255)                                                                                                                                        tcACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAg             >  8:20637/1‑149 (MQ=255)                                                                                                                                        tcACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAg             >  7:36271/1‑149 (MQ=255)                                                                                                                                          aCCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGTGTGGCCCGAGTAGGTGTATCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACa              <  2:255861/146‑1 (MQ=255)                                                                                                                                          aCCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCAAC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCg           <  2:32336/149‑1 (MQ=255)                                                                                                                                           cccTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTATCCGTGGGAGAAGGCACCTTACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCg           <  1:120170/149‑1 (MQ=255)                                                                                                                                            ccTCGACGATGTCGAGCACAGCGTTTGATGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATacac                                        >  6:208038/1‑118 (MQ=255)                                                                                                                                            ccTCGACGATGTCGAGCACAGCGTTTGATGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATacac                                        <  5:208038/118‑1 (MQ=255)                                                                                                                                             cTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTg        >  8:295278/1‑149 (MQ=255)                                                                                                                                              tCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGc       >  6:305075/1‑149 (MQ=255)                                                                                                                                                gacgaTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGc       <  8:36271/147‑1 (MQ=255)                                                                                                                                                gacgaTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCg      >  7:244336/1‑148 (MQ=255)                                                                                                                                                gacgaTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGt     <  2:236747/149‑1 (MQ=255)                                                                                                                                                 acgaTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTa    <  5:305075/149‑1 (MQ=255)                                                                                                                                                 acgaTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTa    <  4:165048/149‑1 (MQ=255)                                                                                                                                                    aTGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTAgg  >  4:168231/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GGATCCGCTACGAACTCAATGGCGGCCATCATGCCCACACCACGCACTTCGCCGACGATTGGCAGCGAGGAGAAGGTCTCCTGCATCGACTTTTGGAAATACTGACCAACCTTGTCGACGTGCGCTGGAATGTTCTCACCCTCGACGATGTCGAGCACAGCGTTTGCTGCCGCGACACCGATAGGGTGGCCCGAGTAGGTGTAGCCGTGGGAGAAGGCACC‑AACGCGATCAGCGCCTTCCTCCATCACCTTATACACCTTCTCACCAACAATCGAAGCCGACAGCGGTGCGTAGG  >  NC_002947/5016182‑5016476

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