Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,630,665	T→A	5.8%	L429Q (CTG→CAG)	PP_3200 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,630,665	0	T	A	5.8%	208.6 / 2.4	69	L429Q (CTG→CAG)	PP_3200	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/27);  new base A (2/2);  total (40/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 = 8.06e-01

GGTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGATGTAC  >  NC_002947/3630520‑3630813                                                                                                                                                  |                                                                                                                                                      ggTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGa                                                                                                                                                    <  2:177400/149‑1 (MQ=255)  gTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAg                                                                                                                                                   <  2:277988/149‑1 (MQ=255)      gTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTAGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCg                                                                                                                                                >  2:248413/1‑148 (MQ=255)          tCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACt                                                                                                                                            >  4:189972/1‑148 (MQ=255)            aaCCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCa                                                                                                                                         >  7:304751/1‑149 (MQ=255)              cccGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTTCTGGCTTGCCAGCCTGGAGCCGAACTTCATc                                                                                                                                       <  2:326950/149‑1 (MQ=255)                   tGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGc                                                                                                                                  >  2:227055/1‑149 (MQ=255)                     gAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCAAGCCTGGAGCCGACATTCATCGCTGcc                                                                                                                                 >  6:281363/1‑148 (MQ=255)                      aaCCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGa                                                                                                                               <  1:227055/149‑1 (MQ=255)                        ccAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGt                                                                                                                             >  1:245797/1‑149 (MQ=255)                        ccAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGt                                                                                                                             >  1:250964/1‑149 (MQ=255)                                  gCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTc                                                                                                                            >  2:201959/1‑140 (MQ=255)                                  gCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTc                                                                                                                            <  1:201959/140‑1 (MQ=255)                                    cTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCg                                                                                                                >  8:318135/1‑149 (MQ=255)                                           aaccaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGaaa                                                                                                          >  1:351607/1‑148 (MQ=255)                                           aaccaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGaaa                                                                                                          <  2:351607/148‑1 (MQ=255)                                           aaccaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAAt                                                                                                         >  1:355941/1‑149 (MQ=255)                                            accaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCa                                                                                                                                         <  5:58160/117‑1 (MQ=255)                                            accaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCa                                                                                                                                         >  6:58160/1‑117 (MQ=255)                                             ccaGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCg                                                                                                       >  1:256548/1‑149 (MQ=255)                                                  aTGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTg                                                                                                   >  3:225563/1‑148 (MQ=255)                                                         ggCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCGTCGGCGAAATCGCCTGGAACCgg                                                                                            <  3:189972/148‑1 (MQ=255)                                                             gTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGtgat                                                                                       <  2:355941/149‑1 (MQ=255)                                                                   gaggcgaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGTCTTTAAAGCCTGGAGCCGAACTTAATCGCTGCCGATTCGA‑GCTTCATCGGAGAAATCGCCTGGAACCGGGTGATGAGCGt                                                                                 >  4:158421/1‑149 (MQ=255)                                                                   gaggcgaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGACTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGt                                                                                 <  5:281363/149‑1 (MQ=255)                                                                       cgaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGc                                                                             <  8:304751/149‑1 (MQ=255)                                                                        gaTCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTg                                                                                                   >  6:220532/1‑126 (MQ=255)                                                                          tCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGc                                                                             >  8:150350/1‑146 (MQ=255)                                                                            gACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGaa                                                                        <  1:248413/149‑1 (MQ=255)                                                                             aCAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCAGCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAAc                                                                       >  5:365089/1‑149 (MQ=255)                                                                               aaTGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGc                                                                     >  8:279766/1‑149 (MQ=255)                                                                                 tGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTACCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                            <  4:290956/124‑1 (MQ=255)                                                                                 tGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTACCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCgg                                                                                            >  3:290956/1‑124 (MQ=255)                                                                                          cccGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAAt                                                                                                         >  1:173911/1‑102 (MQ=255)                                                                                          cccGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAAt                                                                                                         <  2:173911/102‑1 (MQ=255)                                                                                          cccGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGtg                                                                                         <  4:167587/118‑1 (MQ=255)                                                                                          cccGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGtg                                                                                         >  3:167587/1‑118 (MQ=255)                                                                                          cccGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGtt                                                           >  5:193533/1‑148 (MQ=255)                                                                                             ggTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGAcc                                                       <  2:245797/149‑1 (MQ=255)                                                                                             ggTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGAcc                                                       <  2:250964/149‑1 (MQ=255)                                                                                                  aTGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAAc                                                  >  2:109069/1‑149 (MQ=255)                                                                                                   tGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAAc                                                  <  6:365089/148‑1 (MQ=255)                                                                                                     ccACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGcc                                               >  3:128771/1‑149 (MQ=255)                                                                                                      cACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGa                                                                               >  7:57927/1‑116 (MQ=255)                                                                                                      cACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGa                                                                               <  8:57927/116‑1 (MQ=255)                                                                                                      cACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGcc                                               >  1:203136/1‑148 (MQ=255)                                                                                                       aCCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGtt                                                           <  8:330593/135‑1 (MQ=255)                                                                                                       aCCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGtt                                                           >  7:330593/1‑135 (MQ=255)                                                                                                                aCGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGc                                    <  6:193533/149‑1 (MQ=255)                                                                                                                      cACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGa                                                                                                            >  5:242173/1‑71 (MQ=255)                                                                                                                      cACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGa                                                                                                            <  6:242173/71‑1 (MQ=255)                                                                                                                         gCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCt                           >  1:303876/1‑149 (MQ=255)                                                                                                                         gCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCt                           >  7:174858/1‑149 (MQ=255)                                                                                                                           gttcggctcCAGGCTGGCAAGCCAGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAAACGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTgc                         >  6:15352/7‑149 (MQ=255)                                                                                                                           cAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTgc                         >  4:36672/1‑149 (MQ=255)                                                                                                                                tctcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTg                    <  4:128771/149‑1 (MQ=255)                                                                                                                                tctcCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTg                    <  7:150350/149‑1 (MQ=255)                                                                                                                                gctcCAGGCTGGCAAGCCAGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTg                    >  2:223064/2‑149 (MQ=255)                                                                                                                                    caggCTGGCAAGCCAGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGgt                  <  1:223064/145‑1 (MQ=255)                                                                                                                                       gCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGtt               <  4:225563/147‑1 (MQ=255)                                                                                                                                       gCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTc              >  6:348146/1‑148 (MQ=255)                                                                                                                                        ctggcaaGCCAGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGa            <  5:15352/143‑1 (MQ=255)                                                                                                                                        cTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCTAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGa            <  3:158421/149‑1 (MQ=255)                                                                                                                                         ttGCCAGCCTGGAGCCGGACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAg           >  3:155031/1‑149 (MQ=255)                                                                                                                                             cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGa       >  8:85459/1‑149 (MQ=255)                                                                                                                                             cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGa       >  4:48696/1‑149 (MQ=255)                                                                                                                                             cAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGa       <  2:256548/149‑1 (MQ=255)                                                                                                                                              aGCCTGGAGCCGAAATTCATCGCTGCCGAGTCGA‑GCTTACTCGGCGAAATCGCGTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGAGCCTAACGCAGACCGTGACGCCACCATCCTGCGACTGGGGTTCGAGcctag      >  6:356659/1‑146 (MQ=255)                                                                                                                                                ccTGGAGCCGAACTTCATCGCTGCCGAGTCGATGCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGATg     >  1:42795/1‑149 (MQ=255)                                                                                                                                                  tGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGtt               <  4:104341/136‑1 (MQ=255)                                                                                                                                                  tGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGtt               >  3:104341/1‑136 (MQ=255)                                                                                                                                                  tGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGATGTAc  >  1:208947/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GGTGAGTTGTCAACCCGCTGGAACCAGCCGTTGGCCTCGACCAACCAGCATGCGTTGGCCGTGGGCGAGGCGATCGACAATGACGATAACCCGGTTTATGCCACCGGCCGCACGCTGCACGCCAACGTCTCCTGGCTTGCCAGCCTGGAGCCGAACTTCATCGCTGCCGAGTCGA‑GCTTCCTCGGCGAAATCGCCTGGAACCGGGTGATGAGCGTGAGCAAGAACGCCGATGCCGTTGACCCTAACGCCGACCGTGACGCCACCAGCCTGCGCCTGGTGTTCGAGCCGATGTAC  >  NC_002947/3630520‑3630813

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