Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,816,531	C→G	7.7%	A114G (GCC→GGC)	PP_5092 →	NLP/P60 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,816,531	0	C	G	7.7%	200.2 / 7.2	65	A114G (GCC→GGC)	PP_5092	NLP/P60 protein
Reads supporting (aligned to +/- strand):  ref base C (31/29);  new base G (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

TGTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAC  >  NC_002947/5816385‑5816679                                                                                                                                                   |                                                                                                                                                     tGTCGCCGCCATGTCGCCTAACCAAAGCAGCTATGTGCTGAGCCGTGCGGTAAATTTGCTCGGTACTCCTTATGTTTGGCGCGGCAGCAGCCCGAAAAAATGCTTTTACTGCAGCGGGCTTGTCAAATACGCCTTCAACGACGTCGCCg                                                                                                                                                    <  6:216/149‑1 (MQ=255)   tCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGAc                                                                                                                                                  <  8:207073/149‑1 (MQ=255)         ccATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGGCGACGTCgtt                                                                                                                                            <  5:138489/149‑3 (MQ=255)           aTGTCGCCTAAGCAAAGCAGCGTTGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCt                                                                                                                                          >  5:201364/1‑149 (MQ=255)              tCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGcc                                                                                                                                       <  5:56470/149‑1 (MQ=255)                gCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCcgc                                                                                                                                     >  7:1627/1‑149 (MQ=255)                gCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCcgc                                                                                                                                     >  2:10119/1‑149 (MQ=255)                gCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCcgc                                                                                                                                     >  6:7929/1‑149 (MQ=255)                     agcaacgcaGCGATGTGCTGAGCCGTGCGGTAAATGGGCTCGGTACTCCTAATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCAcc                                                                                                                                <  8:298406/149‑1 (MQ=255)                      gcaaagcaGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCt                                                                                                                               >  5:237804/1‑149 (MQ=255)                        aaagcaGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGGCGACGTCgtt                                                                                                                                            <  6:132496/134‑3 (MQ=255)                        aaagcaGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGGCGACGTCgtt                                                                                                                                            >  5:132496/1‑132 (MQ=255)                        aaagcaGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGGCGACGTCgtt                                                                                                                                            >  3:232770/1‑132 (MQ=255)                        aaagcaGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGGCGACGTCgtt                                                                                                                                            <  4:232770/134‑3 (MQ=255)                         aagcaGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTcc                                                                                                                             >  8:200102/1‑148 (MQ=255)                          agcagaGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCaa                                                                                                                           >  6:116535/1‑149 (MQ=255)                               cGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCa                                                                                                                      >  8:170293/1‑149 (MQ=255)                                     gCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGccc                                                                                                                >  1:55969/1‑149 (MQ=255)                                      cTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCg                                                                                                               >  7:190497/1‑149 (MQ=255)                                         aGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCggg                                                                                                             <  8:190497/148‑1 (MQ=255)                                         aGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCggg                                                                                                             <  3:164599/148‑1 (MQ=255)                                                  gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAACGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGt                                                                                                   <  7:170293/149‑1 (MQ=255)                                                  gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCGTGGCCCCGGGCCATggcat                                                                                                   <  5:116535/149‑4 (MQ=255)                                                  gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCgg                                                                                                              >  1:171776/1‑138 (MQ=255)                                                  gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCgg                                                                                                              <  2:171776/138‑1 (MQ=255)                                                      aTGTGCTCCGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAg                                                                                               <  6:238007/149‑1 (MQ=255)                                                         tgCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGt                                                                                             <  4:9552/148‑1 (MQ=255)                                                          gCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTgg                                                                                           >  2:15459/1‑149 (MQ=255)                                                           cTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTATGACTGCAGCGGGCTGGTCAAATACGCCTACAACGACGTCGCCGACGTCGACCTCCCGCGCACCTCCAATGCCATGCCCCGCGGCCACGGCGTGAAGGTGGc                                                                                          >  6:232787/1‑149 (MQ=255)                                                            tCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGc                                                                                          >  8:245684/1‑148 (MQ=255)                                                                    ccTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGa                                                                                 >  6:78535/1‑149 (MQ=255)                                                                     cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     <  7:68286/128‑1 (MQ=255)                                                                     cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACggg                                                                                                     >  8:68286/1‑128 (MQ=255)                                                                     cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGa                                                                                 >  3:227030/1‑148 (MQ=255)                                                                     cTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                >  3:222817/1‑149 (MQ=255)                                                                      ttATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATc                                                                               >  4:210580/1‑149 (MQ=255)                                                                      ttATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATc                                                                               >  3:7655/1‑149 (MQ=255)                                                                           tttGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGTGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGCCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGaa                                                                           <  7:200102/148‑1 (MQ=255)                                                                                 gcggcAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGc                                                                    <  3:210580/149‑1 (MQ=255)                                                                                    gcagcaGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGAc                                                                 >  5:72256/1‑149 (MQ=255)                                                                                    gcagcaGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGAc                                                                 <  4:222817/149‑1 (MQ=255)                                                                                    gcagcaGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGAc                                                                 <  4:227030/149‑1 (MQ=255)                                                                                            ccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCAACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGAGACCTGAtctt                                                         >  3:59290/1‑149 (MQ=255)                                                                                             cGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtcttc                                                        <  5:7929/149‑1 (MQ=255)                                                                                              gAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCAt                                                                                                                     >  2:153541/1‑87 (MQ=255)                                                                                              gAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCAt                                                                                                                     <  1:153541/87‑1 (MQ=255)                                                                                                aaaaaGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCttcttc                                                     <  5:78535/149‑1 (MQ=255)                                                                                                  aaaGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGTGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCaa                                                   <  3:95466/149‑1 (MQ=255)                                                                                                       cTTTGACTGCAGCGTGCTGGTCAAATACTCCTTCAACGACGTCTCCGACGTCGACCTGCAGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGTCGACCTGATCTTCTTCAACATc                                               <  4:59290/148‑1 (MQ=255)                                                                                                          tGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCGACATCAag                                            <  6:237804/148‑1 (MQ=255)                                                                                                           gACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAaga                                           <  6:201364/148‑1 (MQ=255)                                                                                                            aCTGCAGCGGGCTGGTCAAATACGCCTTCAACGAGGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCTCGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCTACATCAagag                                          <  5:232787/148‑1 (MQ=255)                                                                                                              tGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCg                                       >  1:140534/1‑149 (MQ=255)                                                                                                                cAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGt                                      <  2:55969/148‑1 (MQ=255)                                                                                                                        tGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGaa                                                                           <  4:282235/103‑1 (MQ=255)                                                                                                                        tGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGaa                                                                           >  3:282235/1‑103 (MQ=255)                                                                                                                        tGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGCCCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCCAt                             >  6:236991/1‑149 (MQ=255)                                                                                                                                     cTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACg                         >  7:76291/1‑140 (MQ=255)                                                                                                                                     cTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACg                         <  8:76291/140‑1 (MQ=255)                                                                                                                                     cTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATc                 >  2:23371/1‑148 (MQ=255)                                                                                                                                        cAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc              <  2:140534/148‑1 (MQ=255)                                                                                                                                             acgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTggg         >  4:253743/1‑148 (MQ=255)                                                                                                                                              cgTCGCCGACGTCGACCTGCCGCGCACCTCCAAAGCCATGGCGCGGGGCCACGGGGTGAAGGTGGCCATGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTGAACATCAAGAGGCGGCGGGTCAATCAGGTAGCAATCTACGTGGGt        >  4:142368/1‑148 (MQ=255)                                                                                                                                                tCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAAc     <  7:245684/149‑1 (MQ=255)                                                                                                                                                  gCCGACGTCGACCTGCCGCGCAAATCAAATGCCATGGCCCGGGGCCACGGGGTGAAGGTAGCCAAGGGTGATCTGAAGCCTGGCGACCTGAACTTCTTAAAAATAAAGAGCCGTCGGGTCAATCCCGTACACAACAACCAGGGTAACGa   >  8:202692/1‑149 (MQ=255)                                                                                                                                                   ccGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAc  <  1:10119/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TGTCGCCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGCCATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAC  >  NC_002947/5816385‑5816679

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