Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	143,604	C→G	7.5%	S18R (AGC→AGG)	PP_0136 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	143,604	0	C	G	7.5%	208.7 / 6.6	67	S18R (AGC→AGG)	PP_0136	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/27);  new base G (2/3);  total (37/30)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

GCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCGGCAGTTCTCA  >  NC_002947/143456‑143750                                                                                                                                                      |                                                                                                                                                     gcgGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAgc                                                                                                                                                      >  8:267179/1‑149 (MQ=255)       gctgctACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCgcggcg                                                                                                                                                <  7:267179/149‑1 (MQ=255)         tgctACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCgcgg                                                                                                                                                  <  6:273075/145‑1 (MQ=255)         tgctACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCgcgg                                                                                                                                                  >  5:273075/1‑145 (MQ=255)           ctACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCGCTGCCGCCCAGCAGCAGGGCGGCGACCa                                                                                                                                            <  2:45167/149‑1 (MQ=38)           ctACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCggcggc                                                                                                                                              <  3:131872/147‑1 (MQ=255)              cAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATgg                                                                                                                                         >  6:156488/1‑149 (MQ=255)               aGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCCGCCATggc                                                                                                                                        <  1:313983/149‑1 (MQ=255)                 gTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATggcgg                                                                                                                                      >  4:87670/1‑149 (MQ=255)                    tCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTGTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGcc                                                                                                                                    <  5:170153/148‑1 (MQ=255)                      gCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCaa                                                                                                                                  <  8:13056/148‑1 (MQ=255)                       ccAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATg                                                                                                                                >  4:220287/1‑149 (MQ=255)                       ccAAGATCCTGTTCATTGTCTTCCTGGGGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATg                                                                                                                                >  8:31603/1‑149 (MQ=255)                             tCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCGCTGCCGCCCAGCAGCAGGgcggcg                                                                                                                                                >  7:135923/1‑127 (MQ=38)                             tCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCGCTGCCGCCCAGCAGCAGGgcggcg                                                                                                                                                <  8:135923/127‑1 (MQ=38)                                 gTTCATTGTCTTCATGGTGTTTTTTGTAGCATCCTTATTCTTTGGACGCTGACTAGGTTTAATATGAGCAGGAAGCTTCGTTCAGCAGTGGTC‑GCCGCCATGATGCTTGGAGGCAGCGAGGCGGCAATGGCGGCAAATGACGGCCAGGt                                                                                                                      >  5:164065/1‑149 (MQ=255)                                  ttCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTCTGCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAg                                                                                                                        <  1:247728/147‑1 (MQ=255)                                     aTTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCgg                                                                                                                  >  5:10787/1‑149 (MQ=255)                                      ttGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCggg                                                                                                                 >  3:101656/1‑149 (MQ=255)                                          cTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGa                                                                                                             >  1:307361/1‑149 (MQ=255)                                                    ttgtttgtAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCCGTCATTGGCCGCCatg                                                                                                                       <  1:64896/129‑3 (MQ=255)                                                    ttgtttgtAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCCGTCATTGGCCGCCatg                                                                                                                       >  2:64896/1‑127 (MQ=255)                                                       tttgtAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTgctcgg‑‑c                                                                                              >  6:249563/1‑149 (MQ=255)                                                       tttgtAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTgctcgg                                                                                                 <  4:101656/148‑1 (MQ=255)                                                                  cttcttcttTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCg                                                                                     >  6:16772/1‑147 (MQ=255)                                                                   ttcttcttTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAg                                                                                  >  3:88276/1‑149 (MQ=255)                                                                         tttGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCgg                                                                            <  3:220287/149‑1 (MQ=255)                                                                          ttGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCgg                                                                            >  8:249718/1‑148 (MQ=255)                                                                          ttGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCgg                                                                            <  6:10787/148‑1 (MQ=255)                                                                                gcgGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCt                                                                     >  7:147045/1‑149 (MQ=255)                                                                                  ggACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTgg                                                                   >  7:201397/1‑149 (MQ=255)                                                                                    aCGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCa                                                                 <  5:156488/149‑1 (MQ=255)                                                                                             aaGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTgct                                                                                                    <  7:147425/107‑1 (MQ=255)                                                                                             aaGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTgct                                                                                                    >  8:147425/1‑107 (MQ=255)                                                                                             aaGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCt                                                        <  8:201397/149‑1 (MQ=255)                                                                                              aGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCAGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCta                                                       >  1:82287/1‑148 (MQ=255)                                                                                                   aGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAgg                                                  >  8:242650/1‑149 (MQ=255)                                                                                                    gCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGc                                                 >  2:97518/1‑149 (MQ=255)                                                                                                      aGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGa                                               <  8:147045/149‑1 (MQ=255)                                                                                                      aGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGa                                               <  6:124460/149‑1 (MQ=255)                                                                                                             ttCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCa                                                                 >  5:249183/1‑124 (MQ=255)                                                                                                             ttCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCa                                                                 <  6:249183/124‑1 (MQ=255)                                                                                                             ttCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGa                                                              >  2:76998/1‑127 (MQ=255)                                                                                                             ttCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGa                                                              <  1:76998/127‑1 (MQ=255)                                                                                                                gTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTg                                     >  5:13688/1‑149 (MQ=255)                                                                                                                   cAGCCTTGGTC‑GCCTCCCTGCTGCTGGGCGGCAGCTCTGCGGCAATGGCGGCCAATGACGGCCAGGTTCGGGTAGACCAACTGCTCTG‑‑CTCGGACCCGGAGTACCGGTAAACCTGTCAGGACACCCTCTAAGGCGAGGAACGCCTGcc                                   <  6:164065/148‑1 (MQ=255)                                                                                                                   cAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGccc                                  >  8:96044/1‑149 (MQ=255)                                                                                                                     gCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGa                                >  1:163726/1‑149 (MQ=255)                                                                                                                     gCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGa                                >  5:196681/1‑149 (MQ=255)                                                                                                                       cGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCg                                                      >  3:136473/1‑125 (MQ=255)                                                                                                                       cGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCg                                                      <  4:136473/125‑1 (MQ=255)                                                                                                                        gTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTg                             >  1:53075/1‑149 (MQ=255)                                                                                                                           gcc‑gccgcGCTGCCGCCCAGCAGCAGGGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTggg                           >  3:205193/3‑148 (MQ=38)                                                                                                                           gcc‑gccgcGCTGCCGCCCAGCAGCAGGGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTggg                           <  4:205193/146‑1 (MQ=38)                                                                                                                           gTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGGCCCCCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTgg                            <  1:97518/149‑1 (MQ=255)                                                                                                                             c‑gccgcCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCa                                                                 <  3:67028/108‑1 (MQ=255)                                                                                                                             c‑gccgcCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCa                                                                 >  4:67028/1‑108 (MQ=255)                                                                                                                               gccgccCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGggtggt                       <  7:249718/149‑1 (MQ=255)                                                                                                                               gccgccCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGggtggt                       <  7:242650/149‑1 (MQ=255)                                                                                                                                ccgccCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGggtggt                       >  3:277201/1‑148 (MQ=255)                                                                                                                                     ctgctgctgGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTTGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAAAGCCTGCCAGATTGGGTGGTCAATCt                 >  8:237723/1‑149 (MQ=255)                                                                                                                                         tgctgGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACCCCTGCCCGCTTGGGTGGTCAATCTGAcc             >  8:237136/1‑149 (MQ=255)                                                                                                                                            tgGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCgg           >  4:230000/1‑148 (MQ=255)                                                                                                                                             gggCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCgg           >  8:44136/1‑147 (MQ=255)                                                                                                                                              ggcggcAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGAcc             >  8:108647/1‑144 (MQ=255)                                                                                                                                              ggcggcAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGAcc             <  7:108647/144‑1 (MQ=255)                                                                                                                                                   cAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCGGCAGTtctc   <  2:163726/149‑1 (MQ=255)                                                                                                                                                    aGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCGGCAGTTCTca  <  3:230000/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCGGCAGTTCTCA  >  NC_002947/143456‑143750

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