Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	143,594	T→C	8.3%	L15P (CTG→CCG) ‡	PP_0136 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	143,594	0	T	C	8.3%	162.4 / 4.7	60	L15P (CTG→CCG) ‡	PP_0136	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/24);  new base C (2/3);  total (33/27)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCGG  >  NC_002947/143456‑143741                                                                                                                                            |                                                                                                                                                      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                                      <  6:124460/149‑1 (MQ=255)                                                                                                      aGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGa                                      <  8:147045/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                                                     <  1:76998/127‑1 (MQ=255)                                                                                                             ttCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGa                                                     >  2:76998/1‑127 (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                       >  5:196681/1‑149 (MQ=255)                                                                                                                     gCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGa                       >  1:163726/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                                                        >  4:67028/1‑108 (MQ=255)                                                                                                                             c‑gccgcCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCa                                                        <  3:67028/108‑1 (MQ=255)                                                                                                                               gccgccCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGggtggt              <  7:242650/149‑1 (MQ=255)                                                                                                                               gccgccCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGggtggt              <  7:249718/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)                                                                                                                                            |                                                                                                                                                      GCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACCGG  >  NC_002947/143456‑143741

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