Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	874,539	G→C	9.9%	E40Q (GAA→CAA) ‡	PP_0759 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	874,539	0	G	C	9.9%	208.2 / 11.8	71	E40Q (GAA→CAA) ‡	PP_0759	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (36/28);  new base C (4/3);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.53e-01

CCACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTG  >  NC_002947/874392‑874686                                                                                                                                                    |                                                                                                                                                    ccACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCg                                                                                                                                                     >  6:242480/1‑148 (MQ=255)     tGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGc                                                                                                                                                <  3:219534/149‑1 (MQ=255)     tGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTAAGGCGTAAGTCGTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGc                                                                                                                                                >  6:407520/1‑149 (MQ=255)       gCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGcc                                                                                                                                               <  7:265841/148‑1 (MQ=255)        cGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCaa                                                                                                                                             <  6:151424/149‑1 (MQ=255)         gACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCaa                                                                                                                                             >  4:214888/1‑148 (MQ=255)         gACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAAc                                                                                                                                            <  5:242480/149‑1 (MQ=255)          aCTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCaa                                                                                                                                             >  4:99667/1‑147 (MQ=255)            tCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGc                                                                                                                                                <  3:245804/142‑1 (MQ=255)            tCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGc                                                                                                                                                >  4:245804/1‑142 (MQ=255)            tCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAg                                                                                                                                          >  7:78191/1‑148 (MQ=255)              ggCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCTGTGCTGTTGGCGTCGAACGCCAACAGTAc                                                                                                                                       <  3:99667/149‑1 (MQ=255)              ggCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTTCGACGCCAACAgc                                                                                                                                         <  5:162966/147‑2 (MQ=255)              ggCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTACTGTTGGCGTTCGACGCCAACAgca                                                                                                                                        <  8:320284/148‑3 (MQ=255)               gCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACa                                                                                                                                      <  5:285602/149‑1 (MQ=255)                  tAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGt                                                                                                                                   >  5:153338/1‑149 (MQ=255)                           gcgcATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGc                                                                                                                          <  3:199270/149‑1 (MQ=255)                           gcgcATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGc                                                                                                                          <  3:214888/149‑1 (MQ=255)                                tgtgGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCg                                                                                                                     <  1:79302/149‑1 (MQ=255)                                    gATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGAcc                                                                                                                  <  6:153338/148‑1 (MQ=255)                                      tAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                               <  4:295481/149‑1 (MQ=255)                                       aGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAAAACCACCAGCGGCCGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGAt                                                                                                               >  7:84394/1‑148 (MQ=255)                                        gggCGTGGTTTCATGCGTAAGGCTTTTGTCATCAGCCGGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGa                                                                                                             <  6:254813/149‑1 (MQ=255)                                        gggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGa                                                                                                             <  6:249099/149‑1 (MQ=255)                                         ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATg                                                                                                                                    >  7:227442/1‑125 (MQ=255)                                         ggAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATg                                                                                                                                    <  8:227442/125‑1 (MQ=255)                                          gAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTAc                                                                                                                                       >  2:101062/1‑121 (MQ=255)                                          gAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTAc                                                                                                                                       <  1:101062/121‑1 (MQ=255)                                           aGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGt                                                                                                          >  5:42831/1‑149 (MQ=255)                                               ttttCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAAc                                                                                                      <  1:336802/149‑1 (MQ=255)                                                       cGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACAAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTa                                                                                              >  5:278074/1‑149 (MQ=255)                                                         tAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTAcg                                                                                            >  3:234883/1‑149 (MQ=255)                                                             tCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCCGCATGCTCTCGACCGATGAGGTCAACCAGATGTAcgcg                                                                                          >  8:131320/1‑147 (MQ=255)                                                              cTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAg                                                                                       >  7:296273/1‑149 (MQ=255)                                                                     catcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTAt                                                                                 >  2:405204/1‑148 (MQ=255)                                                                      atcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGATGTAAACCAGATGTACGCGCAGTCGTAt                                                                                 >  8:181003/1‑147 (MQ=255)                                                                      atcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTAt                                                                                 >  6:7/1‑147 (MQ=255)                                                                      atcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATca                                                                               <  8:78191/149‑1 (MQ=255)                                                                       tcaGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATcag                                                                              >  3:36166/1‑149 (MQ=255)                                                                         aGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATcag                                                                              >  2:239800/1‑147 (MQ=255)                                                                          gCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATcag                                                                              >  3:138064/1‑146 (MQ=255)                                                                           ccTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGa                                                                          <  1:362855/149‑1 (MQ=255)                                                                              gctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGAccc                                                                       >  4:380796/1‑149 (MQ=255)                                                                               ctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCt                                                                      >  4:233796/1‑149 (MQ=255)                                                                               ctgAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCt                                                                      >  2:389916/1‑149 (MQ=255)                                                                                  aGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTgg                                                                    <  6:42831/148‑1 (MQ=255)                                                                                  aGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTgg                                                                    <  3:97273/148‑1 (MQ=255)                                                                                  aGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATAACACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCTACCGATGAGGTCAACCAGATGTACGCGCAGACGTATCAGAAGACCCTggg                                                                   >  4:154396/1‑149 (MQ=255)                                                                                       tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAg                                                              >  4:230925/1‑149 (MQ=255)                                                                                       tGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAg                                                              >  6:63160/1‑149 (MQ=255)                                                                                              ctgctgGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAg                                                       >  5:340942/1‑149 (MQ=255)                                                                                               tgctgGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAg                                                       >  4:225420/1‑148 (MQ=255)                                                                                                     ccGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGgt                                                <  8:296273/149‑1 (MQ=255)                                                                                                      cGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGgtg                                               <  3:233796/149‑1 (MQ=255)                                                                                                                gCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATCAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACaa                                     <  5:407520/149‑1 (MQ=255)                                                                                                                    tGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACACTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTc                                  >  1:133998/1‑148 (MQ=255)                                                                                                                          ccaccaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGcc                           >  8:76849/1‑149 (MQ=255)                                                                                                                           caccaGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGcc                           >  3:25019/1‑148 (MQ=255)                                                                                                                              caGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCaa                                                                                                       <  4:85708/69‑1 (MQ=255)                                                                                                                              caGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCaa                                                                                                       >  3:85708/1‑69 (MQ=255)                                                                                                                              caGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGagg                       >  8:88034/1‑147 (MQ=255)                                                                                                                                    atgtacTGTTGGCGTTCGACGCCAACAGCACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGc                 >  6:223753/6‑149 (MQ=255)                                                                                                                                    atgtacTGTTGGCGTTCGACGCCAACAGCACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGc                 >  1:359520/6‑149 (MQ=255)                                                                                                                                     ggTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCg                <  3:380796/149‑1 (MQ=255)                                                                                                                                       tGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCgtg              >  7:19389/1‑149 (MQ=255)                                                                                                                                       tGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCgtg              >  8:207870/1‑149 (MQ=255)                                                                                                                                         cTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGtgtg            <  5:7/149‑1 (MQ=255)                                                                                                                                          tGTTGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGtgtg            >  4:398303/1‑148 (MQ=255)                                                                                                                                           gTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGtgtg            >  2:235308/1‑147 (MQ=255)                                                                                                                                            ttGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTgg                                                                    >  4:14338/1‑90 (MQ=255)                                                                                                                                            ttGGCGTTCGACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTgg                                                                    <  3:14338/90‑1 (MQ=255)                                                                                                                                                  tCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGTGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTGCAGCGCCGATGCCAAGCGTGTGCAGGCCAtt   <  3:154396/149‑1 (MQ=255)                                                                                                                                                    gAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTg  <  4:138064/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CCACTGGCGACTCGGCCTAGGATACGGCGCATGTGGATAGGGAGTGTTTTCATGCGTAAGTCTTTTGTCATCAGCCTGCTGAGTGCTGGCATCCTGCTGGCCGGCTGCCAGGCGGTGAATACCACCAGCGGCGGTGCTGTTGGCGTCGAACGCCAACAGTACATGTTCAGCATGCTCTCGACCGATGAGGTCAACCAGATGTACGCGCAGTCGTATCAGCAGACCCTGGGTGAAGCATCGAGCAAGGGTGTGCTCGACAAGTCCAGCGCCGATGCCAAGCGTGTGCAGGCCATTG  >  NC_002947/874392‑874686

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