Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	143,586	C→G	8.2%	A12A (GCC→GCG)	PP_0136 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	143,586	0	C	G	8.2%	181.6 / 6.9	61	A12A (GCC→GCG)	PP_0136	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/23);  new base G (2/3);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

TTCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGT  >  NC_002947/143444‑143729                                                                                                                                                |                                                                                                                                                  ttCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCctgctg                                                                                                                                             >  4:46469/1‑149 (MQ=255) ttCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCctgctg                                                                                                                                             >  2:247728/1‑149 (MQ=255) ttCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCctgc                                                                                                                                               <  7:114687/147‑1 (MQ=255) ttCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTACGTTCAGCCGTGGTC‑GCCGCCctgc                                                                                                                                               >  8:114687/1‑147 (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                                                                                                  >  2:64896/1‑127 (MQ=255)                                                                ttgtttgtAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCCGTCATTGGCCGCCatg                                                                                                  <  1:64896/129‑3 (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                                            <  6:249183/124‑1 (MQ=255)                                                                                                                         ttCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCa                                            >  5:249183/1‑124 (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                                 <  4:136473/125‑1 (MQ=255)                                                                                                                                   cGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCg                                 >  3:136473/1‑125 (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)                                                                                                                                                |                                                                                                                                                  TTCGGTGGTATCGCGGGCGCTGCTACAGGTATCGCCAAGATCCTGTTCATTGTCTTCCTGGTGTTGTTTGTAGCCTCCTTCTTCTTTGGACGCGGACGAGGTTGAAGATGAGCAGGAAGCTTCGTTCAGCCGTGGTC‑GCCGCCCTGCTGCTGGGCGGCAGCGCGGCGGCCATGGCGGCCAATGACGGCCAGGTTCGGGTCGACCAACTGCTCGG‑‑CTCGGACCCGGAGTACCGGGAAACCTGGCAGGACACCCTCGAAGGCGAGGAACGCCTGCCCGATTGGGTGGT  >  NC_002947/143444‑143729

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