Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,574,852	A→C	6.6%	D10A (GAC→GCC)	PP_2258 →	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,574,852	0	A	C	6.6%	204.4 / 2.5	68	D10A (GAC→GCC)	PP_2258	sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (26/37);  new base C (3/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCCAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAA  >  NC_002947/2574706‑2574991                                                                                                                                                   |                                                                                                                                            gCCCAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGAcg                                                                                                                                           >  6:237453/1‑149 (MQ=255)  cccAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGAcgc                                                                                                                                          >  4:283928/1‑149 (MQ=255)    cAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCgg                                                                                                                                        <  6:61896/149‑1 (MQ=255)    cAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCgg                                                                                                                                        <  5:202650/149‑1 (MQ=255)     aGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGTAAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGc                                                                                                                                       <  1:80636/149‑1 (MQ=255)        acaGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGCGTCCACGGCcgc                                                                                                                                    <  7:161006/149‑1 (MQ=255)         caGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCcgc                                                                                                                                    >  7:297847/1‑148 (MQ=255)         caGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCcgc                                                                                                                                    >  2:16761/1‑148 (MQ=255)              cccATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCtt                                                                                                                              >  4:43950/1‑149 (MQ=255)                cATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGGGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGCCAGGATGCTTTATGGATGAGAACTCTCGCGCGGCCGTGGACGCGGCCGCGATCtttt                                                                                                                            <  1:176062/149‑1 (MQ=255)                 aTAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTc                                                                                                                           <  6:155493/149‑1 (MQ=255)                     cgcCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGaa                                                                                                                       <  1:196341/149‑1 (MQ=255)                      gcCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGaaa                                                                                                                      >  8:253530/1‑149 (MQ=255)                       ccAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGAGCTTTTCGGaaaa                                                                                                                     >  4:6755/1‑148 (MQ=255)                           cgccgccgcCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGAc                                                                                                                 <  5:240374/149‑1 (MQ=255)                                ccgccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCg                                                                                                            <  1:210404/149‑1 (MQ=255)                                  gccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTg                                                                                                          <  8:53174/149‑1 (MQ=255)                                       gcgcgGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGa                                                                                                       <  3:129045/147‑1 (MQ=255)                                       gcgcgGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGa                                                                                                       <  8:297847/147‑1 (MQ=255)                                          cgGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACAt                                                                                                   >  8:16028/1‑148 (MQ=255)                                             tACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCATTTAGGAAACCCACCTGCGTGGGAACAacac                                                                                                >  2:324599/1‑148 (MQ=255)                                                cGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTAc                                                                                            >  7:171558/1‑149 (MQ=255)                                                 gAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACg                                                                                           >  5:253932/1‑149 (MQ=255)                                                  aaCTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGt                                                                                          <  5:237453/149‑1 (MQ=255)                                                    cTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTc                                                                                         >  8:142472/1‑148 (MQ=255)                                                       gAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAAt                                                                                      <  6:179824/148‑1 (MQ=255)                                                       gAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAAt                                                                                      <  7:253530/148‑1 (MQ=255)                                                        aaCACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATca                                                                                    <  3:77310/149‑1 (MQ=255)                                                        aaCACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATca                                                                                    <  1:218702/149‑1 (MQ=255)                                                         acacGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATcag                                                                                   <  1:3083/149‑1 (MQ=255)                                                           acGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATcagca                                                                                 <  7:142472/149‑1 (MQ=255)                                                                  cgcCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACTCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTgc                                                                          <  7:245826/149‑1 (MQ=255)                                                                   gcCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTgca                                                                         <  5:262855/149‑1 (MQ=255)                                                                   gcCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTgca                                                                         <  3:283928/149‑1 (MQ=255)                                                                   gcCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTg                                                                           >  2:94487/1‑147 (MQ=255)                                                                    cccgATTGCGCGACCTGGTGGACCCCAACAACCGCTAACGCAGGATGCTTTATCGATGATAACTATCGCTCGGCATTGGACGCGGCCCCGATCTTTTCCTAAACCGACCTGCTTGGGAACATCACAAACGTCAATCAGCAGTTCTgcag                                                                        >  2:296240/4‑149 (MQ=255)                                                                    ccGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTgcag                                                                        <  1:324599/149‑1 (MQ=255)                                                                         ttGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTgca                                                                         <  5:55916/143‑1 (MQ=255)                                                                         ttGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTgca                                                                         >  6:55916/1‑143 (MQ=255)                                                                              cgcCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGc                                                              <  3:43950/149‑1 (MQ=255)                                                                               gcCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCt                                                             <  8:64980/149‑1 (MQ=255)                                                                                cccAGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGaa                                                                                                      <  7:76249/107‑1 (MQ=255)                                                                                cccAGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGaa                                                                                                      >  8:76249/1‑107 (MQ=255)                                                                                         aGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGt                                                     >  7:137635/1‑147 (MQ=255)                                                                                          gcgcTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAg                                                  <  4:190631/149‑1 (MQ=255)                                                                                          gcgcTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAg                                                  <  6:253932/149‑1 (MQ=255)                                                                                                  aCCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGACAATCAGCAGTTCTgca                                                                         >  3:307945/1‑118 (MQ=255)                                                                                                  aCCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGACAATCAGCAGTTCTgca                                                                         <  4:307945/118‑1 (MQ=255)                                                                                                     gCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTggg                                       <  3:101975/149‑1 (MQ=255)                                                                                                      cTAAGGCAGTATGCTTTATGGATGAAAACTATCGCGCGGCCGTGGACGCGCCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGAAGCATCTCTGGCTACAGCCGTGAGGAACTGCTgggg                                      >  2:26239/1‑149 (MQ=255)                                                                                                      cTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTgggg                                      <  1:16761/149‑1 (MQ=255)                                                                                                               gATGCTTTATGGATGATAACTATCGTGCGGCAGTGGACGCGGCCGCGAAATTTTAGGTAAACTACCATACTGGGAACATCACCTACCTCAATCATCAGTACAGCAGCATCTCGCGCTACAGCCGTGAGGAACTGCTGGGTGCCAacgac                             >  6:80651/1‑149 (MQ=255)                                                                                                                    tttATGGATGATAACTATCGCGCGGCCGTGGACGCGGCAGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCAt                        >  7:116425/1‑149 (MQ=255)                                                                                                                      tATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATc                       <  8:171558/148‑1 (MQ=255)                                                                                                                          gatgatAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCaa                  >  3:86432/1‑149 (MQ=255)                                                                                                                           atgatAACTATCGCGCGGCCGTGGACGCGGCCGCCATCTTTTCGGAAACCGACCTGCGAGGGAACATCACCTACCTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGTGGAACTGCTGGGGGCCAACCACCGCATCCTCAAt                 >  6:263497/1‑149 (MQ=255)                                                                                                                            tgatAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGACCTGCTGGGGGCCa                                  >  5:182299/1‑131 (MQ=255)                                                                                                                            tgatAACTAGCGCGCGGCCGGGGCCGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCa                                  <  6:182299/131‑1 (MQ=255)                                                                                                                               tAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCg             <  8:116425/149‑1 (MQ=255)                                                                                                                                 aCTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCggg           <  4:306070/149‑1 (MQ=255)                                                                                                                                 aCTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCggg           >  3:306070/1‑149 (MQ=255)                                                                                                                                   tATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCt         >  3:314573/1‑149 (MQ=255)                                                                                                                                     gatcgcgGCCGCGTCCACGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGc       >  6:209647/4‑149 (MQ=255)                                                                                                                                       tcgcgGCCGCGTCCACGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCAt     >  7:100329/2‑149 (MQ=255)                                                                                                                                       tcgcgGCCGCGTCCACGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCAt     >  3:290718/2‑149 (MQ=255)                                                                                                                                       gcgcgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCTGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCa      <  1:26239/148‑1 (MQ=255)                                                                                                                                          cgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGaa  <  4:86432/149‑1 (MQ=255)                                                                                                                                          cgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGaa  <  7:16028/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                            GCCCAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAA  >  NC_002947/2574706‑2574991

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