Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,544	G→C	10.1%	P13A (CCC→GCC)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,544	0	G	C	10.1%	189.0 / 12.7	70	P13A (CCC→GCC)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base C (4/3);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.85e-01

CCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGC  >  NC_002947/5518397‑5518689                                                                                                                                                     |                                                                                                                                                  ccAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCggg                                                                                                                                                   >  2:325347/1‑148 (MQ=255)     ttCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGAAGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGc                                                                                                                                              >  2:183823/1‑149 (MQ=255)       cGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGcc                                                                                                                                             >  2:97693/1‑148 (MQ=255)        gAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGc                                                                                                                                           <  8:226594/149‑1 (MQ=255)         aaCACCAGCAGACCGGCCTGCAGGCTTGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                          <  3:24458/149‑1 (MQ=255)         aaCACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCg                                                                                                                                            >  1:51582/1‑147 (MQ=255)           cacCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGt                                                                                                                                        <  1:325347/149‑1 (MQ=255)            acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTc                                                                                                                                       >  5:101953/1‑149 (MQ=255)            acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGa                                                                                                                                                  >  1:155475/1‑138 (MQ=255)            acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGa                                                                                                                                                  <  2:155475/138‑1 (MQ=255)              cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTg                                                                                                                                     <  6:101953/149‑1 (MQ=255)              cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGCCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGAGTATACAGCCTCAACGCCACACAGCGCGAATGCCGCaaaat                                                                                                                                      >  1:119794/1‑145 (MQ=255)               agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGc                                                                                                                                    <  5:72283/149‑1 (MQ=255)               agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGc                                                                                                                                    >  8:28448/1‑149 (MQ=255)                 cagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                                  <  7:28448/149‑1 (MQ=255)                        tccAGCAGGGTGGTGATCAACTTCAGGGTCAAAGCGAGCAGCAGGGCGGCGATGACGATGATCAGGTTGGCGAGGGTGATCTGTCTGTGGGTAAAGG‑TCATTCAGCCTCGACGGCAGACAGCTGTAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                           <  5:361345/148‑1 (MQ=255)                          cagcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggcg                                                                                                                         <  1:45229/149‑1 (MQ=255)                                  tGGGCAGCAACTTCAGGTGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGCCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                           <  5:182186/139‑1 (MQ=255)                                  tGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGCCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                           >  6:182186/1‑139 (MQ=255)                                   gggCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCa                                                                                                                >  4:297250/1‑149 (MQ=255)                                      cagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGAGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                             >  4:293446/1‑149 (MQ=255)                                       agcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                            <  4:82951/149‑1 (MQ=255)                                       agcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                            >  4:204075/1‑149 (MQ=255)                                          aaCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGccgtc                                                                                                                                <  2:272953/126‑4 (MQ=38)                                          aaCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGccgtc                                                                                                                                >  1:272953/1‑123 (MQ=38)                                           aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                   <  5:296231/138‑1 (MQ=255)                                           aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                   >  6:296231/1‑138 (MQ=255)                                               cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCa                                                                                                        <  5:93106/145‑1 (MQ=255)                                               cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCa                                                                                                        >  6:93106/1‑145 (MQ=255)                                                   ggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGTAGGTCAttt                                                                                                >  5:285884/1‑149 (MQ=255)                                                    gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                               <  1:183823/149‑1 (MQ=255)                                                    gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                               >  1:149683/1‑149 (MQ=255)                                                       aaGCCAGCAGCAGCCCGGCCATTACGATGATCAGGCTGGCGAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGAGGGCAGAGAGCGGTAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                            <  1:97693/149‑1 (MQ=255)                                                       aaGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                            <  3:11807/149‑1 (MQ=255)                                                         gCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGTCGGCAGACAGCGGGAAGGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttctt                                                                                          <  7:157570/149‑1 (MQ=255)                                                           cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        <  1:8709/149‑1 (MQ=255)                                                           cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        <  2:51582/149‑1 (MQ=255)                                                             gcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGTAGGTCATTTCTTCTTCAgg                                                                                      <  6:285884/149‑1 (MQ=255)                                                             gcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAgg                                                                                      <  3:297250/149‑1 (MQ=255)                                                                 cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttctt                                                                                             >  1:34218/1‑138 (MQ=255)                                                                 cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttctt                                                                                             <  2:34218/138‑1 (MQ=255)                                                                   gcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCg                                                                                >  7:146843/1‑149 (MQ=255)                                                                    cccGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGc                                                                               >  8:70816/1‑149 (MQ=255)                                                                       ggCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCa                                                                            <  2:149683/149‑1 (MQ=255)                                                                        gCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCAt                                                                           >  7:92053/1‑149 (MQ=255)                                                                                  tgatCTGGCTGGCCAGGGTGTTCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCTGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTg                                                                 <  2:119794/149‑1 (MQ=255)                                                                                   gatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                >  3:76914/1‑149 (MQ=255)                                                                                    atCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                <  8:92053/148‑1 (MQ=255)                                                                                     tCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgcgc                                                              <  8:146843/149‑1 (MQ=255)                                                                                                aGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCtt                                                   >  1:173368/1‑149 (MQ=255)                                                                                                   gTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAg                                                >  8:118308/1‑149 (MQ=255)                                                                                                     gATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGc                                                                     <  7:210884/126‑1 (MQ=255)                                                                                                     gATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGc                                                                     >  8:210884/1‑126 (MQ=255)                                                                                                        cTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        >  8:274989/1‑104 (MQ=255)                                                                                                        cTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                        <  7:274989/104‑1 (MQ=255)                                                                                                         tGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCttt                                          >  3:20985/1‑149 (MQ=255)                                                                                                             ggggggTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtg                                      <  7:118308/149‑1 (MQ=255)                                                                                                               ggggTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                    <  4:76914/149‑1 (MQ=255)                                                                                                                  gTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGc                                 <  2:173368/149‑1 (MQ=255)                                                                                                                       gg‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAg                            >  8:88393/1‑149 (MQ=255)                                                                                                                               cAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGc                      <  4:20985/148‑1 (MQ=255)                                                                                                                                aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                                  <  3:102302/135‑1 (MQ=255)                                                                                                                                aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                                  >  4:102302/1‑135 (MQ=255)                                                                                                                                       ctggCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCa             >  8:225392/3‑149 (MQ=255)                                                                                                                                          gCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGGGGTGGCGTCAGCATTGGAGGTATTCaata          >  3:26475/1‑146 (MQ=255)                                                                                                                                            aGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtgg                                        >  2:126337/1‑117 (MQ=255)                                                                                                                                            aGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtgg                                        <  1:126337/117‑1 (MQ=255)                                                                                                                                             gACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttttt       >  6:95453/1‑149 (MQ=255)                                                                                                                                                aGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCtt    <  3:204075/149‑1 (MQ=255)                                                                                                                                                  cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGc  >  7:16252/1‑149 (MQ=255)                                                                                                                                                  cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGc  >  8:288163/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGC  >  NC_002947/5518397‑5518689

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