Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,863,240	T→A	7.6%	L309Q (CTG→CAG)	PP_1666 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,863,240	0	T	A	7.6%	176.4 / 3.3	62	L309Q (CTG→CAG)	PP_1666	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (34/23);  new base A (3/2);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGA  >  NC_002947/1863092‑1863385                                                                                                                                                      |                                                                                                                                                   gTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcct                                                                                                                                                    >  6:776/1‑149 (MQ=255)    cAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTACAGGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctgg                                                                                                                                                  <  1:82997/149‑1 (MQ=255)        tGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAg                                                                                                                                             <  2:81893/149‑1 (MQ=255)             ggCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTg                                                                                                                                          >  2:372320/1‑147 (MQ=255)             ggCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTg                                                                                                                                          >  2:350829/1‑147 (MQ=255)              gCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAg                                                                                                                                       >  1:184417/1‑149 (MQ=255)              gCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAg                                                                                                                                       <  1:236267/149‑1 (MQ=255)              gCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAg                                                                                                                                       <  5:312116/149‑1 (MQ=255)               cATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAg                                                                                                                                       >  1:155720/1‑148 (MQ=255)                aTGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAgg                                                                                                                                      >  1:254857/1‑148 (MQ=255)                 tGAACCTGCAGCGCTATGCCGACCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGaa                                                                                                                                    <  5:333453/149‑1 (MQ=255)                      cTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGcc                                                                                                                               >  4:292455/1‑149 (MQ=255)                         cAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCgg                                                                                                                             <  8:299892/148‑1 (MQ=255)                         cAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCg                                                                                                                              <  3:329839/147‑1 (MQ=255)                                 tGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCa                                                                                                                      <  1:233247/147‑1 (MQ=255)                                    cGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCTCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCAc                                                                                                                 >  3:299573/1‑149 (MQ=255)                                        cTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGCCCAAGCTGCAGGAAGTGcc                                                                                                                               >  2:55319/1‑131 (MQ=255)                                        cTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGCCCAAGCTGCAGGAAGTGcc                                                                                                                               <  1:55319/131‑1 (MQ=255)                                             ccGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGc                                                                                                        >  8:172401/1‑149 (MQ=255)                                             ccGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGATACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGc                                                                                                        >  3:221100/1‑149 (MQ=255)                                                     ttGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGc                                                                                                    <  4:66376/145‑1 (MQ=255)                                                     ttGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGc                                                                                                    >  3:66376/1‑145 (MQ=255)                                                                   ccACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctggc                                                                                                                                                 <  6:133217/86‑1 (MQ=255)                                                                   ccACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTtggcctggc                                                                                                                                                 >  5:133217/1‑86 (MQ=255)                                                                   ccACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGt                                                                                 <  7:172401/149‑1 (MQ=255)                                                                   ccACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGt                                                                                 >  7:172442/1‑149 (MQ=255)                                                                    cACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCAGGCCAAGCAGCAGGAAGTGCCGGTCGACCAGGCCCCGGCGGTGCCGGCCATGCCAGAAC‑CGGCaaac                                                                                >  7:234570/1‑146 (MQ=255)                                                                       gCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCTCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGAc                                                                             <  4:299573/149‑1 (MQ=255)                                                                         ttGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACTGCCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGAcgc                                                                           >  6:364304/1‑149 (MQ=255)                                                                         ttGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGAcgc                                                                           >  5:223262/1‑149 (MQ=255)                                                                            cAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGAcgccgc                                                                        <  5:776/149‑1 (MQ=255)                                                                            cAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGAcgccg                                                                         <  7:329778/148‑1 (MQ=255)                                                                                   cTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCAGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCcgcg                                                                 >  6:376438/1‑149 (MQ=255)                                                                                    tGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCAGGCCAAGCTgg                                                                                                                                         >  7:290668/1‑76 (MQ=255)                                                                                    tGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCAGGCCAAGCTgg                                                                                                                                         <  8:290668/77‑2 (MQ=255)                                                                                      aaGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAg                                                              >  6:171582/1‑149 (MQ=255)                                                                                      aaGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAg                                                              >  6:165908/1‑149 (MQ=255)                                                                                        ggCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGTTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGcc                                                            <  1:372320/149‑1 (MQ=255)                                                                                           aCTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGa                                                                              <  4:301250/128‑1 (MQ=255)                                                                                           aCTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGa                                                                              >  3:301250/1‑128 (MQ=255)                                                                                                ccTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTg                                                    >  6:121624/1‑149 (MQ=255)                                                                                                 cTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTgc                                                   >  8:72547/1‑149 (MQ=255)                                                                                                    aCCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGca                                                <  5:12464/149‑1 (MQ=255)                                                                                                      cTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGcagc                                              >  6:105123/1‑149 (MQ=255)                                                                                                                   gggCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTAGGTGGGGTCCCAGTGCGGGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAg                                   >  1:312249/1‑149 (MQ=255)                                                                                                                    ggCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGcc                                 <  8:1478/149‑1 (MQ=255)                                                                                                                     gCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTCCAGCAACTGCCGGTCGACCAGGCACCCGCGGTGCCGGCCATGCCAGTAC‑CGGCAGAGGACGCCGCAGCCCCGCAGCCCCCGCCTGCCCAGCCAGCACCCAAcccc                                >  3:121174/1‑148 (MQ=255)                                                                                                                        aTCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTc                             <  1:350829/149‑1 (MQ=255)                                                                                                                         tCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTc                             >  1:316646/1‑148 (MQ=255)                                                                                                                         tCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTc                             >  4:289626/1‑148 (MQ=255)                                                                                                                           gTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGAc                          >  4:252645/1‑149 (MQ=255)                                                                                                                            tGAACAGCTACGTGCCCAGCTTGGCCTGGCCCAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACg                         >  6:353555/1‑149 (MQ=255)                                                                                                                                ccgCTCCGGGCCCGGGTTGGCCTGGCCCAGCTGCAGGTAGTGCCGGGCGCGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCtt                     <  8:234570/147‑1 (MQ=255)                                                                                                                                 aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGt                               <  7:152109/138‑1 (MQ=255)                                                                                                                                 aGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGAAGGCACCGGCGGTGCCGGCCATGACAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCATCACCCAAcccct                               >  8:152109/1‑136 (MQ=255)                                                                                                                                   cTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCg                  >  8:142211/1‑149 (MQ=255)                                                                                                                                     aCGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGtt                >  6:250662/1‑149 (MQ=255)                                                                                                                                      cGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGttt               <  7:142211/149‑1 (MQ=255)                                                                                                                                          tgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTgc                                                   >  6:294722/3‑109 (MQ=255)                                                                                                                                          tgcAGCTTGGCCAGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTgc                                                   <  5:294722/107‑1 (MQ=255)                                                                                                                                          cccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGc           >  2:291368/1‑149 (MQ=255)                                                                                                                                           ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGttt               <  3:147643/144‑1 (MQ=255)                                                                                                                                           ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGttt               >  4:147643/1‑144 (MQ=255)                                                                                                                                           ccAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCt          <  5:105123/149‑1 (MQ=255)                                                                                                                                            cAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTg         >  2:196006/1‑149 (MQ=255)                                                                                                                                                ttGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGt       >  1:12447/1‑147 (MQ=255)                                                                                                                                                  ggcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCTCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAg    <  4:221100/148‑1 (MQ=255)                                                                                                                                                   gcctggccAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGa  <  3:289626/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GTGCAAGTGCAGGGCATGAACCTGCAGCGCTATGCCGAGCTGGCCCGCGTGCTTGAGCCCTATGGCCCACGCTTGCAGGCGGCTGAAGGCACTACCCTGACCTAC‑GGGGTTACGGGCAATCGTGAACAGCTACGTGCCCAGCTTGGCCTGGCCAAGCTGCAGGAAGTGCCGGTCGAGCAGGCACCGGCGGTGCCGGCCATGCCAGTAC‑CGGCAGTGGACGCCGCAGCCGCGCAGCCCGCGCCTGCGCAGCCAGCACCCAAGCCGTTCGACGGCTTGCGTTTTCGCTGGTAAGGA  >  NC_002947/1863092‑1863385

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