Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	143,591	T→C	8.3%	L14P (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,591	0	T	C	8.3%	165.4 / 4.6	60	L14P (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.53e-01

TTCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACC  >  NC_002947/143444‑143739                                                                                                                                                     |                                                                                                                                                       ttCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCctgctg                                                                                                                                                       >  4:46469/1‑149 (MQ=255) ttCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCctgctg                                                                                                                                                       >  2:247728/1‑149 (MQ=255)             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                                                                                         >  8:147425/1‑107 (MQ=255)                                                                                                         aaGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTgct                                                                                         <  7:147425/107‑1 (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                                                   >  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                     >  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)                                                                                                                                                     |                                                                                                                                                       TTCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGTCAATCTGACC  >  NC_002947/143444‑143739

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