Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,574,852	A→C	5.7%	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	5.7%	208.5 / 3.5	70	D10A (GAC→GCC)	PP_2258	sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (34/32);  new base C (2/2);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.77e-01

CCCAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAAC  >  NC_002947/2574707‑2574992                                                                                                                                                  |                                                                                                                                              cccAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGAcgc                                                                                                                                            >  4:174535/1‑149 (MQ=255)  ccAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGAcgcg                                                                                                                                           >  8:7202/1‑149 (MQ=255)   cAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGTCCACGGCcgcgc                                                                                                                                          >  8:262127/1‑148 (MQ=255)   cAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCgg                                                                                                                                          >  1:174304/1‑149 (MQ=255)       acaGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGTCCACGGCcgcgcgata                                                                                                                                      <  7:262127/149‑6 (MQ=255)        caGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGCGTCCACGGCcgcg                                                                                                                                     <  6:105075/149‑1 (MQ=255)              ccATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCttt                                                                                                                               >  7:40177/1‑149 (MQ=255)              ccATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCttt                                                                                                                               >  6:149250/1‑149 (MQ=255)               cATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCtttt                                                                                                                              <  2:235828/149‑1 (MQ=255)                aTAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTc                                                                                                                             <  5:120137/149‑1 (MQ=255)                aTAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTc                                                                                                                             <  3:226288/149‑1 (MQ=255)                 tAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCg                                                                                                                            >  6:107950/1‑149 (MQ=255)                 tAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCg                                                                                                                            <  5:149250/149‑1 (MQ=255)                     gcCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGaaa                                                                                                                        >  8:99966/1‑149 (MQ=255)                        aaCGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCg                                                                                                                     >  7:348281/1‑149 (MQ=255)                          cgccgccgcCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTTCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGGTGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGAc                                                                                                                   <  6:302551/149‑1 (MQ=255)                           gccgccgccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGAcc                                                                                                                  >  4:132338/1‑149 (MQ=255)                             cgccgccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAAc                                                                                                                       >  3:356336/1‑142 (MQ=255)                             cgccgccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAAc                                                                                                                       <  4:356336/142‑1 (MQ=255)                              gccgcgTCGCGCGGTACCGTACTGGATCACGTGTGGGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCTGCCGTGGACGCGGCCGCGATCTTTTCTGAAACCGACCTGc                                                                                                               <  8:168086/149‑1 (MQ=255)                              gccgccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGc                                                                                                               <  8:209098/149‑1 (MQ=255)                              gccgccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGc                                                                                                               >  8:337328/1‑149 (MQ=255)                                cgccTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGt                                                                                                             <  5:324223/149‑1 (MQ=255)                                            tACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCAcc                                                                                                 >  8:240289/1‑149 (MQ=255)                                            tACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCAcc                                                                                                 <  8:268970/149‑1 (MQ=255)                                              ccGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAAACGACCTGCGTGGGAACATCACCTa                                                                                               >  5:116574/1‑149 (MQ=255)                                                gAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACg                                                                                             >  2:115946/1‑149 (MQ=255)                                                 aaCTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGt                                                                                            >  7:101154/1‑149 (MQ=255)                                                      gAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAAt                                                                                        >  8:12201/1‑148 (MQ=255)                                                      gAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTCCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAAt                                                                                        >  4:258518/1‑148 (MQ=255)                                                         cacGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATcag                                                                                     >  4:362435/1‑148 (MQ=255)                                                            gtgtGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGt                                                                                 >  7:349454/1‑149 (MQ=255)                                                               tgCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCt                                                                              <  7:7202/149‑1 (MQ=255)                                                                  gcCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTg‑ca                                                                          <  5:107950/149‑1 (MQ=255)                                                                  gcCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTg‑ca                                                                          <  7:99966/149‑1 (MQ=255)                                                                    cGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTg‑cagc                                                                        <  2:174304/149‑1 (MQ=255)                                                                           cgcgcCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCgg                                                                 <  4:344250/149‑1 (MQ=255)                                                                            gcgcCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCggg                                                                <  7:12201/149‑1 (MQ=255)                                                                              gcCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGc                                                               >  8:132789/1‑148 (MQ=255)                                                                                ccTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTa                                                             >  2:226327/1‑148 (MQ=255)                                                                                       gAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTg                                                     <  3:132338/149‑1 (MQ=255)                                                                                        aGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGa                                                    <  8:101154/149‑1 (MQ=255)                                                                                         gcgcTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAg                                                   >  5:10483/1‑149 (MQ=255)                                                                                         gcgcTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAg                                                   <  7:240289/149‑1 (MQ=255)                                                                                                    gCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTggg                                        <  3:362435/149‑1 (MQ=255)                                                                                                          gCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCaa                                  >  4:280263/1‑149 (MQ=255)                                                                                                             ggATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGATACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAacca                               <  8:349454/149‑1 (MQ=255)                                                                                                              gATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAaccac                              >  7:345194/1‑149 (MQ=255)                                                                                                               aTGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAaccacc                             <  7:337328/149‑1 (MQ=255)                                                                                                                  cTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGc                           <  7:132789/148‑1 (MQ=255)                                                                                                                   tttATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAaccacc                             >  8:206487/1‑145 (MQ=255)                                                                                                                   tttATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAaccacc                             <  7:206487/145‑1 (MQ=255)                                                                                                                   tttATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCAt                         >  7:81881/1‑149 (MQ=255)                                                                                                                   tttATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCAt                         >  3:365287/1‑149 (MQ=255)                                                                                                                      aTGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATc                        <  4:45613/147‑1 (MQ=255)                                                                                                                      aTGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATc                        >  3:45613/1‑147 (MQ=255)                                                                                                                       tGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCt                      >  3:357243/1‑148 (MQ=255)                                                                                                                         gatgatAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCaa                   >  2:28538/1‑149 (MQ=255)                                                                                                                            gatAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAAc                                               <  2:182981/118‑1 (MQ=255)                                                                                                                            gatAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAAc                                               >  1:182981/1‑118 (MQ=255)                                                                                                                              tAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCg              <  3:280263/149‑1 (MQ=255)                                                                                                                                  tATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCt          >  3:279669/1‑149 (MQ=255)                                                                                                                                     atcgcgGCCGCGTCCACGGCCGCGCTCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCa       >  3:222894/3‑149 (MQ=255)                                                                                                                                      gcgcgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTGTTAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCa       <  1:115946/149‑1 (MQ=255)                                                                                                                                      gcgcgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGt                                                      <  6:291705/101‑1 (MQ=255)                                                                                                                                      gcgcgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGt                                                      >  5:291705/1‑101 (MQ=255)                                                                                                                                      gcgcgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCAt      >  7:260503/1‑149 (MQ=255)                                                                                                                                      gcgcgGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCAt      <  3:174535/149‑1 (MQ=255)                                                                                                                                          ggCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAAc  <  1:28538/149‑1 (MQ=255)                                                                                                                                          ggCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAAc  <  8:40177/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                              CCCAGCACAGCGCCCATAGCGCCAACGCCGCCGCCTCGCGCGGTACCGAACTGGAACACGTGTGCGCCGGATTGCGCGCCCTGGTGGAGCGCTTCAACCGCTAAGGCAGGATGCTTTATGGATGATAACTATCGCGCGGCCGTGGACGCGGCCGCGATCTTTTCGGAAACCGACCTGCGTGGGAACATCACCTACGTCAATCAGCAGTTCTG‑CAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAAC  >  NC_002947/2574707‑2574992

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