Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,853,771	T→C	6.8%	intergenic (+54/+28)	PP_3401 → / ← PP_3402	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,853,771	0	T	C	6.8%	210.3 / 3.5	73	intergenic (+54/+28)	PP_3401/PP_3402	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (38/29);  new base C (2/3);  total (41/32)
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.87e-01

CTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAGGGTCAC  >  NC_002947/3853624‑3853918                                                                                                                                                    |                                                                                                                                                    cTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTg                                                                                                                                                    >  7:184654/1‑149 (MQ=255) cTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTg                                                                                                                                                    <  4:126371/149‑1 (MQ=255)  ttCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTgg                                                                                                                                                   >  3:27652/1‑149 (MQ=255)        tGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGACCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGAc                                                                                                                                             >  7:185604/1‑149 (MQ=255)          ggATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGAc                                                                                                                                             <  4:319957/147‑1 (MQ=255)          ggATGGGCGGAG‑TATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGTCAGCCCAGGCTGGCCgg                                                                                                                                          <  2:323877/149‑1 (MQ=255)           gATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGTCAGCCCAGGCTGGCCgg                                                                                                                                          <  1:162310/149‑1 (MQ=255)           gATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCgg                                                                                                                                          >  1:174365/1‑149 (MQ=255)            aTGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCGGGGCTGCCCggg                                                                                                                                         >  3:84410/1‑149 (MQ=255)                         tCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGgcgc                                                                                                                            >  2:284483/1‑149 (MQ=255)                               cGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAg                                                                                                                      <  2:174365/149‑1 (MQ=255)                                gTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAgg                                                                                                                     <  7:79745/149‑1 (MQ=255)                                gTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAgg                                                                                                                     >  5:215225/1‑149 (MQ=255)                                 tGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGGCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGc                                                                                                                    <  8:74308/149‑1 (MQ=255)                                     tCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATc                                                                                                                >  5:40502/1‑149 (MQ=255)                                        cgAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCa                                                                                                                                      >  8:197122/1‑124 (MQ=255)                                        cgAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCa                                                                                                                                      <  7:197122/124‑1 (MQ=255)                                         gagaAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCt                                                                                                            >  6:386022/1‑149 (MQ=255)                                         gagaAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCt                                                                                                            <  4:27652/149‑1 (MQ=255)                                          agaAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAg                                                                                                                      <  3:103969/138‑1 (MQ=255)                                          agaAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAg                                                                                                                      >  4:103969/1‑138 (MQ=255)                                             aGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGa                                                                                                         <  6:35595/148‑1 (MQ=255)                                              gggAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGt                                                                                                       >  2:119220/1‑149 (MQ=255)                                               ggAAAGTGCGGCGACGGTGGCCTGAACCTGCTACTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTg                                                                                                      >  1:241924/1‑149 (MQ=255)                                                 aaaGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGACCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGgc                                                                                                    <  8:185604/149‑1 (MQ=255)                                                        ggcgACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGt                                                                                             >  8:203476/1‑149 (MQ=255)                                                                tGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCa                                                                                     <  1:49743/149‑1 (MQ=255)                                                                  gCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCac                                                                                    <  8:227927/148‑1 (MQ=255)                                                                   ccTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCt                                                                                         >  6:317683/1‑142 (MQ=255)                                                                   ccTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCt                                                                                         <  5:317683/142‑1 (MQ=255)                                                                    cTGAACCTGCTGCTGGAACTGCAATAAGCAGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCAAGACAAGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGCTGGATATCaccc                                                                                  >  4:60668/1‑146 (MQ=255)                                                                       aaCCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCacgcaac                                                                               <  5:386022/148‑1 (MQ=255)                                                                        aCCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCacgcaacgc                                                                             <  1:119220/149‑1 (MQ=255)                                                                         ccTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGTCAGCCCAGGCTGGCCgg                                                                                                                                          >  2:181852/1‑87 (MQ=39)                                                                          ctgctgctgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCAt                                                                           >  6:332152/1‑149 (MQ=255)                                                                                ctgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgct                                                                     >  4:143883/1‑149 (MQ=255)                                                                                ctgGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgct                                                                     <  5:86285/149‑1 (MQ=255)                                                                                  ggAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgctgc                                                                   >  7:367041/1‑149 (MQ=255)                                                                                   gAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAgctgct                                                                  >  4:179000/1‑149 (MQ=255)                                                                                     aCTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTc                                                                 >  7:276272/1‑148 (MQ=255)                                                                                         cAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTc                                                                 <  2:261781/144‑1 (MQ=255)                                                                                         cAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTc                                                                 >  1:261781/1‑144 (MQ=255)                                                                                           aTAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAgg                                                          <  6:40502/149‑1 (MQ=255)                                                                                             aaGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGtt                                                        >  7:181989/1‑149 (MQ=255)                                                                                               gCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGa                                                      <  8:181989/149‑1 (MQ=255)                                                                                                  gATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGAc                                                     >  4:29743/1‑147 (MQ=255)                                                                                                       gcggcgGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCa                                              >  4:119889/1‑149 (MQ=255)                                                                                                               cGTGGCCTCGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCTCTCAGGCAATCGCGTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGGAGCTGCTCCACCAGGTTGACCAGGTCACGACAACt                                      <  3:60668/149‑1 (MQ=255)                                                                                                               cGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACt                                      <  6:215225/149‑1 (MQ=255)                                                                                                                gcagCCACGAGCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTg                                     >  7:146669/4‑149 (MQ=255)                                                                                                                gTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTg                                     <  1:179987/149‑1 (MQ=255)                                                                                                                   gCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCGCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGtt                                   <  8:184654/148‑1 (MQ=255)                                                                                                                     cACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAg                                <  6:11274/149‑1 (MQ=255)                                                                                                                     cACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAg                                >  5:11274/1‑149 (MQ=255)                                                                                                                        gCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTc                             >  1:292725/1‑149 (MQ=255)                                                                                                                        gCTAGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACACAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTc                             >  7:353941/1‑149 (MQ=255)                                                                                                                            gTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTg                          <  2:343559/148‑1 (MQ=255)                                                                                                                              ggCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATc                                                                                      <  3:244944/86‑1 (MQ=255)                                                                                                                              ggCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATc                                                                                      >  4:244944/1‑86 (MQ=255)                                                                                                                                 tgctTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAgcgc                    >  1:348323/1‑149 (MQ=255)                                                                                                                                         cccGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACACAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGt            <  8:353941/149‑1 (MQ=255)                                                                                                                                          ccGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCacgca                                                                                 <  5:17085/79‑1 (MQ=255)                                                                                                                                          ccGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCacgca                                                                                 >  6:17085/1‑79 (MQ=255)                                                                                                                                          ccGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGtt           >  5:338829/1‑149 (MQ=255)                                                                                                                                          ccGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGtt           >  2:385331/1‑149 (MQ=255)                                                                                                                                           cGGTCAGCCCAGGCTGGCCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAg                                                  <  2:296197/109‑1 (MQ=255)                                                                                                                                           cGGTCAGCCCAGGCTGGCCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAg                                                  >  1:296197/1‑109 (MQ=255)                                                                                                                                               cAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAggg      >  6:161381/1‑149 (MQ=255)                                                                                                                                                  ccTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACt                                      >  4:287185/1‑114 (MQ=255)                                                                                                                                                  ccTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACt                                      <  3:287185/114‑1 (MQ=255)                                                                                                                                                  ccTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAGGGTCa   >  3:323374/1‑149 (MQ=255)                                                                                                                                                  ccTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAGGGTCa   >  5:981/1‑149 (MQ=255)                                                                                                                                                   cTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAGGGTCAc  >  5:327682/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CTTCGTCTGGGATGGGCGGAGTTATCACGGCGTGGTTCGCGAGAAGGGAAAGTGCGGCGACGGTGGCCTGAACCTGCTGCTGGAACTGCAATAAGCCGATCAGCGGCGGCCGTGGCCACGCTCGTGGCTGCTTATACCCGGCCAGCCTGGGCTGACCGGGGCACGACAGGCGCGCTCAGGCAATCGCCTTGAGTGGCGCCGGGTGGCTATCACGCAACGCATGCAGCTGCTCCACCAGGTTGACCAGGTCACGACAACTGTTGAGGTCTTGCAGCGCGACCCCGTTGAGGGTCAC  >  NC_002947/3853624‑3853918

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