Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,436,515	G→T	6.3%	R136L (CGC→CTC)	gpI →	tail formation protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,436,515	0	G	T	6.3%	198.7 / 3.1	63	R136L (CGC→CTC)	gpI	tail formation protein
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base T (2/2);  total (32/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 = 6.95e-01

CAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTT  >  NC_002947/3436369‑3436660                                                                                                                                                   |                                                                                                                                                  cAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAAGGTCGAAAATTTCTCCGCAGCTGGAGGAGCGCTTgcgcgc                                                                                                                                                >  1:965068/1‑150 (MQ=255)   aGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGccc                                                                                                                                              <  1:712192/150‑1 (MQ=255)    ggCGGTTCGGGGCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCt                                                                                                                                             <  1:47253/150‑1 (MQ=255)       ggTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGt                                                                                                                                          <  1:285999/150‑1 (MQ=255)       ggTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGt                                                                                                                                          <  1:372113/150‑1 (MQ=255)                         ttACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGa                                                                                                                        <  1:571499/150‑1 (MQ=255)                           aCAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAgc                                                                                                                      <  1:201397/150‑1 (MQ=255)                           aCAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAgc                                                                                                                      <  2:750854/150‑1 (MQ=255)                             aaGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACGGTCGGGGTGACAGCTCGCTCATGTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAgcgc                                                                                                                    >  1:78863/1‑150 (MQ=255)                                       gcgAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGcgacg                                                                                                                                >  2:424642/1‑128 (MQ=255)                                       gcgAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGcgacg                                                                                                                                <  1:424642/128‑1 (MQ=255)                                           gaTACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGAACATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCtc                                                                                                      >  2:357016/1‑149 (MQ=255)                                               cccAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCgg                                                                                                  >  1:832423/1‑150 (MQ=255)                                               cccAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCgg                                                                                                  >  2:251968/1‑150 (MQ=255)                                                ccAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGc                                                                                                 >  2:861766/1‑150 (MQ=255)                                                ccAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGGGAAGGCTCGAACATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGc                                                                                                 >  2:130630/1‑150 (MQ=255)                                                 cAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGt                                                                                                                                          <  1:994304/108‑1 (MQ=255)                                                 cAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGt                                                                                                                                          >  2:994304/1‑108 (MQ=255)                                                 cAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGcga                                                                                                                                  <  2:316114/116‑1 (MQ=255)                                                 cAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGcga                                                                                                                                  >  1:316114/1‑116 (MQ=255)                                                  aaCGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTGGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTc                                                                                               >  1:457742/1‑150 (MQ=255)                                                          gcCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCg                                                                                       <  1:581565/150‑1 (MQ=255)                                                          gcCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCg                                                                                       >  2:885342/1‑150 (MQ=255)                                                                   gTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGGAATTCGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTcggcga                                                                              >  1:674276/1‑150 (MQ=255)                                                                    ttCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCgacgac                                                                             <  2:468798/150‑1 (MQ=255)                                                                    ttCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCgacgac                                                                             >  1:772106/1‑150 (MQ=255)                                                                     tCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGacgacg                                                                            >  2:6310/1‑150 (MQ=255)                                                                       cAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCCCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGc                                                                          >  2:858313/1‑150 (MQ=255)                                                                        aGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGcc                                                                         >  1:182442/1‑150 (MQ=255)                                                                             aCGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGc                                                                          <  2:1098351/144‑1 (MQ=255)                                                                             aCGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGc                                                                          >  1:1098351/1‑144 (MQ=255)                                                                             aCGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCtggtg                                                                    <  1:251968/150‑1 (MQ=255)                                                                              cGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCtggtgg                                                                   <  1:6310/150‑1 (MQ=255)                                                                                  ttGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTAc                                                               <  1:487535/150‑1 (MQ=255)                                                                                  ttGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTAc                                                               >  2:623405/1‑150 (MQ=255)                                                                                    gggCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCg                                                             <  2:654987/150‑1 (MQ=255)                                                                                                aCCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCACTGGTAAACGCAACGAAGATCTAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGGCGGCCTGGTGGCTACCGGAGCCAGTCAgg                                                 <  1:688279/150‑1 (MQ=255)                                                                                                    gccggGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCa                                             <  2:832423/150‑1 (MQ=255)                                                                                                    gccggGTGAAGGCTCGATCATTTCTCCGCAGATGGAGATGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCa                                             >  1:213805/1‑150 (MQ=255)                                                                                                       gggTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAgct                                          <  1:517309/150‑1 (MQ=255)                                                                                                             aGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGAGGCCGGAGTCGTTGCGGCCGGAGCCAGTCAGGGCCAGCtgctgcc                                    >  2:88805/1‑149 (MQ=255)                                                                                                             aGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCa                                    <  1:698903/150‑1 (MQ=255)                                                                                                               gCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCAcc                                  >  2:123433/1‑150 (MQ=255)                                                                                                                 tCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCTCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGc                                >  1:930262/1‑150 (MQ=255)                                                                                                                 tCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCTCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGc                                >  1:930326/1‑150 (MQ=255)                                                                                                                     tCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGc                                               <  2:469080/131‑1 (MQ=255)                                                                                                                     tCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGc                                               >  1:469080/1‑131 (MQ=255)                                                                                                                        tttCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACgg                         >  1:928290/1‑150 (MQ=255)                                                                                                                            tcCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGa                     >  1:158948/1‑150 (MQ=255)                                                                                                                               gCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGc                                               >  2:322203/1‑121 (MQ=255)                                                                                                                               gCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGc                                               <  1:322203/121‑1 (MQ=255)                                                                                                                                  gCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCGGGTTCGACGGCGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAAACTGTGGAGCCTgc               >  1:228718/1‑150 (MQ=255)                                                                                                                                       aggagCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGc                                                         >  1:757267/1‑103 (MQ=255)                                                                                                                                       aggagCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGc                                                         <  2:757267/103‑1 (MQ=255)                                                                                                                                        ggagCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGt         >  2:43524/1‑150 (MQ=255)                                                                                                                                            cgcTTGCTCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCg     <  2:930326/150‑1 (MQ=255)                                                                                                                                            cgcTTGCTCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCg     <  2:930262/150‑1 (MQ=255)                                                                                                                                             gcTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCgg    <  2:78863/150‑1 (MQ=255)                                                                                                                                              cTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGt   <  1:417539/150‑1 (MQ=255)                                                                                                                                              cTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGt   >  1:743142/1‑150 (MQ=255)                                                                                                                                               ttGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCAcc                                  <  1:190099/118‑1 (MQ=255)                                                                                                                                               ttGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCAcc                                  >  2:190099/1‑118 (MQ=255)                                                                                                                                               ttGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGtt  <  2:928290/150‑1 (MQ=255)                                                                                                                                               ttGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGtt  <  2:1014005/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CAAGGCGGTTCGGGTCAACGCGGATTACAAGGAATGGCGCGAGATACCCAACGCGGCGCCGTACACGTTCCAGGTCACGGCTTGGGCCAACGAGAACCGGCCGGGTGAAGGCTCGATCATTTCTCCGCAGCTGGAGGAGCGCTTGCGCGCCCTGGTCGATGCGACGAAGAACGAGCGCAGCCACTACACCTTCCGGCTCGGTGCCCGGTTCGACGACGGCCTGGTGGCTACCGGAGCCAGTCAGGGCCAGCTGCTGCACCGCAAGACGGTGGAGCCTGCGCCCGTTCCGGTT  >  NC_002947/3436369‑3436660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑