Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,128,232	A→T	10.3%	Q1201H (CAA→CAT)	PP_1887 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,128,232	0	A	T	10.3%	215.8 / 13.2	78	Q1201H (CAA→CAT)	PP_1887	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (34/36);  new base T (4/4);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GCCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGT  >  NC_002947/2128088‑2128373                                                                                                                                                  |                                                                                                                                              gcCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAAt                                                                                                                                             >  3:218930/1‑147 (MQ=255) gcCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGa                                                                                                                                           <  6:27147/149‑1 (MQ=255) gcCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGa                                                                                                                                           <  4:82307/149‑1 (MQ=255)  cccGTGGACCCCTCCGGTCAAACCATTGCGAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAg                                                                                                                                          <  2:50723/149‑1 (MQ=255)   ccGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGc                                                                                                                                         >  2:300624/1‑149 (MQ=255)     gTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACaaa                                                                                                                                              <  8:27975/142‑1 (MQ=255)     gTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACaaa                                                                                                                                              >  7:27975/1‑142 (MQ=255)         aCCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGAt                                                                                                                                    >  1:112373/1‑148 (MQ=255)          ccccccCGGTCAAACCACTGCCAGGCAGATATTCGGCAT‑CGATGGCCTGGACAACATCCTCTCGGCGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATTTGTGCCGGATcg                                                                                                                                  <  8:62714/149‑2 (MQ=255)          ccccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATTTGTGCCGGAt                                                                                                                                    <  5:303969/147‑1 (MQ=255)          ccccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCTCATTTGTGCCGGAt                                                                                                                                    >  6:303969/1‑147 (MQ=255)           cccTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCAt                                                                                                                                 >  8:230466/1‑149 (MQ=255)             cTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGTACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCg                                                                                                                                       <  7:14439/141‑1 (MQ=255)             cTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCg                                                                                                                                       >  8:14439/1‑141 (MQ=255)                cGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCg                                                                                                                             >  4:389852/1‑148 (MQ=255)                cGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGc                                                                                                                            <  2:221115/149‑1 (MQ=255)                     aaaCCATTGCCAGGCAGATATTCGGCTTCCGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTACGCCTGa                                                                                                                        >  1:230489/1‑149 (MQ=255)                      aaCCTTTGCCAGGAAGATTTGTGTCGT‑CGATTGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGAt                                                                                                                       <  6:284106/148‑1 (MQ=255)                      aaCCATTGCCAGGCAGATATTCGGCTT‑CCATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATg                                                                                                                      >  1:230435/1‑149 (MQ=255)                       aCCATTGCCAGGCAGCTATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATgc                                                                                                                     <  8:349656/149‑1 (MQ=255)                       aCCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGACCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATgc                                                                                                                     >  3:81668/1‑149 (MQ=255)                         cATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATgcgc                                                                                                                   <  2:112373/149‑1 (MQ=255)                         cATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATgcgc                                                                                                                   >  2:219742/1‑149 (MQ=255)                               cAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAGTAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTg                                                                                                              >  2:286717/1‑148 (MQ=255)                                     gATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCt                                                                                                       <  7:230466/149‑1 (MQ=255)                                      atatTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGCGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTa                                                                                                      >  6:60587/1‑149 (MQ=255)                                        atTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTAt                                                                                                     >  8:108908/1‑148 (MQ=255)                                          tCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGAc                                                                                                  >  6:97064/1‑149 (MQ=255)                                           cGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAATAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACg                                                                                                 >  7:322504/1‑149 (MQ=255)                                             gCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTc                                                                                               >  8:117339/1‑149 (MQ=255)                                              cTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTc                                                                                               >  6:7374/1‑148 (MQ=255)                                               tt‑CGATGGCCTGGACAACATCCTCTCGTTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACg                                                                                                 >  1:223512/1‑145 (MQ=255)                                                  cGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGAt                                                                                            >  4:19586/1‑148 (MQ=255)                                                  cGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGAt                                                                                            >  4:238508/1‑148 (MQ=255)                                                   gATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCg                                                                                                                             <  3:319368/114‑1 (MQ=255)                                                   gATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCg                                                                                                                             >  4:319368/1‑114 (MQ=255)                                                   gATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCa                                                                                          >  4:98308/1‑149 (MQ=255)                                                    aTGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCAc                                                                                         <  1:219742/149‑1 (MQ=255)                                                    aTGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCAc                                                                                         <  4:292871/149‑1 (MQ=255)                                                    aTGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCAc                                                                                         >  1:395622/1‑149 (MQ=255)                                                      ggCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGa                                                                                       >  7:260799/1‑149 (MQ=255)                                                       gCCTGGACAACATCCTCTCGGTGACCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGAc                                                                                      <  4:81668/149‑1 (MQ=255)                                                        ccTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGAcc                                                                                     <  5:324860/149‑1 (MQ=255)                                                            gACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTg                                                                                   <  2:395622/147‑1 (MQ=255)                                                            gACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTg                                                                                   <  7:375014/147‑1 (MQ=255)                                                                 caTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggag                                                                                >  8:180172/1‑145 (MQ=255)                                                                 caTCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTggag                                                                                <  7:180172/145‑1 (MQ=255)                                                                      tctcGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAgttgtt                                                                        >  2:54140/1‑148 (MQ=255)                                                                           gTGATCACCTACAGACCGGACGGGTCATGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGa                                                                                                                                     >  8:157509/1‑82 (MQ=255)                                                                           gTGATCACCTACAGACCGGACGGGTCATGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGa                                                                                                                                     <  7:157509/82‑1 (MQ=255)                                                                            tGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATc                                                                 <  1:54140/149‑1 (MQ=255)                                                                            tGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATc                                                                 <  3:238508/149‑1 (MQ=255)                                                                               tCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTc                                                              <  5:7374/149‑1 (MQ=255)                                                                                     aCAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATAc                                                        >  5:145634/1‑149 (MQ=255)                                                                                     aCAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATAc                                                        >  2:331848/1‑149 (MQ=255)                                                                                       aGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGa                                                      <  2:230435/149‑1 (MQ=255)                                                                                         aCCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGCCCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATc                                                    >  5:312814/1‑149 (MQ=255)                                                                                         aCCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTc                                                              >  3:86669/1‑139 (MQ=255)                                                                                         aCCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTc                                                              <  4:86669/139‑1 (MQ=255)                                                                                          ccGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAATAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATcc                                                   <  6:212878/149‑1 (MQ=255)                                                                                          ccGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATcc                                                   <  2:223512/149‑1 (MQ=255)                                                                                               cGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGt                                              <  3:98308/149‑1 (MQ=255)                                                                                                 ggTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTg                                            <  1:300624/149‑1 (MQ=255)                                                                                                  gTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTgg                                           <  8:322504/149‑1 (MQ=255)                                                                                                  gTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGAGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTgg                                           <  7:117339/149‑1 (MQ=255)                                                                                                      tGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTAt                                       <  7:108908/149‑1 (MQ=255)                                                                                                            cGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGAACACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGAttt                                 <  4:218930/149‑1 (MQ=255)                                                                                                              gACCTTGTACGAATTCAAGAGTAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGc                               <  1:286717/149‑1 (MQ=255)                                                                                                              gACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTg                                >  7:116068/1‑148 (MQ=255)                                                                                                              gACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGc                               <  3:19586/149‑1 (MQ=255)                                                                                                                ccTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCAt                             <  3:389852/149‑1 (MQ=255)                                                                                                                    gTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTa                           <  5:97064/147‑1 (MQ=255)                                                                                                                        gAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCCCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACgg                     <  1:331848/149‑1 (MQ=255)                                                                                                                              aaGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGATATCGATTTGCATTAAGACGGAAATgg               >  8:191411/1‑149 (MQ=255)                                                                                                                                         tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAAc           >  4:221746/1‑142 (MQ=255)                                                                                                                                         tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAAc           <  3:221746/142‑1 (MQ=255)                                                                                                                                         tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGAt      >  5:25926/1‑147 (MQ=255)                                                                                                                                         tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATc     >  3:12250/1‑148 (MQ=255)                                                                                                                                         tCCGGCTCATTTGTGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATc     <  4:12250/148‑1 (MQ=255)                                                                                                                                            ggCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGt  >  7:267515/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                              GCCCGTGGACCCCTCCGGTCAAACCATTGCCAGGCAGATATTCGGCTT‑CGATGGCCTGGACAACATCCTCTCGGTGATCACCTACAGACCGGACGGGTCGTGGCAGCGGACCTTGTACGAATTCAAGAATAGCGATCCGGCACAAATGAGCCGGATCATTCCGCCTGATGCGCTACTGGAGACCTATGACGTCGATCACGACCTGGAGGAGTTGTTCAAGATCGTCAAATACGATCCGCAGTTGGCTATCGATTTGCATTACGACGGCAATGGCAACCTGATCAGT  >  NC_002947/2128088‑2128373

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