Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,539	A→T	10.0%	A14A (GCT→GCA)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,539	0	A	T	10.0%	191.2 / 10.6	70	A14A (GCT→GCA)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (4/3);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 7.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

GTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTT  >  NC_002947/5518392‑5518687                                                                                                                                                     |                                                                                                                                                     gTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAg                                                                                                                                                     >  2:45229/1‑149 (MQ=255)    gACCAGTTCGAAAACAAGCATACCGGCCATCAGTCTTGGCATCAACTTCAGTGGAAAAGCCAGCATCAGCCCGGCCATGACGATGATCAGGCTGTCCAGGTTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACTGAAGACAGCg                                                                                                                                                   <  6:194151/148‑1 (MQ=255)      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:149683/1‑149 (MQ=255)                                                         gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                             <  1:183823/149‑1 (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                                                                   >  8:210884/1‑126 (MQ=255)                                                                                                          gATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGc                                                                   <  7:210884/126‑1 (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                                >  4:102302/1‑135 (MQ=255)                                                                                                                                     aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                                <  3:102302/135‑1 (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)                                                                                                                                                     |                                                                                                                                                     GTTGACCAGTTCGAACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTT  >  NC_002947/5518392‑5518687

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