Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	81,989	C→T	11.1%	T100M (ACG→ATG)	PP_0071 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	81,989	0	C	T	11.1%	184.0 / 12.2	72	T100M (ACG→ATG)	PP_0071	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base T (4/4);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTACAG  >  NC_002947/81842‑82136                                                                                                                                                    |                                                                                                                                                    cccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGc                                                                                                                                                   >  1:779612/1‑150 (MQ=255)  ccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGcc                                                                                                                                                  >  2:159094/1‑150 (MQ=255)  ccTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGcc                                                                                                                                                  >  2:678786/1‑150 (MQ=255)      taCCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCATCg                                                                                                                                              <  1:159094/150‑1 (MQ=255)         cTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCGTCGACTCGATGCCATCgagt                                                                                                                                           <  2:228810/150‑4 (MQ=255)         cTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCg                                                                                                                                           >  2:107594/1‑150 (MQ=255)             gcAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCAc                                                                                                                                       >  2:1073612/1‑150 (MQ=255)                  gcTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCa                                                                                                                                  <  2:708434/150‑1 (MQ=255)                  gcTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCa                                                                                                                                  <  1:690636/150‑1 (MQ=255)                        gAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTcc                                                                                                                            >  2:460017/1‑150 (MQ=255)                             cGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGc                                                                                                                       >  1:435646/1‑150 (MQ=255)                                   cAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAgcg                                                                                                                 >  2:746647/1‑150 (MQ=255)                                    aaaGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCacg                                                                                                                         <  2:323062/141‑3 (MQ=255)                                    aaaGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                         >  1:323062/1‑141 (MQ=255)                                        cTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCgctgct                                                                                                            <  1:1073612/150‑1 (MQ=255)                                         tttGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTa                                                                                                           >  2:208938/1‑150 (MQ=255)                                           tGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATc                                                                                                         <  1:857421/150‑1 (MQ=255)                                             gACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                                <  1:982520/125‑1 (MQ=255)                                             gACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGc                                                                                                                                >  2:982520/1‑125 (MQ=255)                                               cTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAgg                                                                                                     <  1:107594/150‑1 (MQ=255)                                                tACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGc                                                                                                    <  1:242382/150‑1 (MQ=255)                                                       gTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                             <  1:108382/150‑1 (MQ=255)                                                       gTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                             <  2:779612/150‑1 (MQ=255)                                                       gTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGCCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCa                                                                                             >  1:521799/1‑150 (MQ=255)                                                        tGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAg                                                                                            <  2:457561/150‑1 (MQ=255)                                                        tGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAg                                                                                            >  1:899187/1‑150 (MQ=255)                                                        tGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGCAATAGCGGGCTTGCCCGCGCCTGCACTCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAg                                                                                            >  1:444914/1‑150 (MQ=255)                                                            gaagGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTAcc                                                                                        >  1:219797/1‑150 (MQ=255)                                                               ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGAc                                                                                                                                                     <  2:1132386/86‑1 (MQ=255)                                                               ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGAc                                                                                                                                                     >  1:1132386/1‑86 (MQ=255)                                                               ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCCAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCtgt                                                                                     <  2:345896/150‑1 (MQ=255)                                                               ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGACGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCtgt                                                                                     >  2:424364/1‑150 (MQ=255)                                                               ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGACGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATa                                                                                                                    >  2:904166/1‑119 (MQ=255)                                                               ggCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGACGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATa                                                                                                                    <  1:904166/119‑1 (MQ=255)                                                                      gCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                         <  1:70014/107‑1 (MQ=255)                                                                      gCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                         >  2:70014/1‑107 (MQ=255)                                                                       cATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                             <  2:899187/150‑1 (MQ=255)                                                                       cATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                             <  1:460017/150‑1 (MQ=255)                                                                                gTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCCTCGAGTCGACGCCATAGCCGTCAGGCCCAGCCTTCAAGGCAGAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCATGCCcg                                                                    >  1:537080/1‑150 (MQ=255)                                                                                         gCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGCGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCg                                                           <  1:1150073/150‑1 (MQ=255)                                                                                         gCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCg                                                           <  1:678786/150‑1 (MQ=255)                                                                                           aTGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGc                                                         >  2:549708/1‑150 (MQ=255)                                                                                            tGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCa                                                        >  1:428513/1‑150 (MQ=255)                                                                                             gCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCaa                                                       >  2:47372/1‑150 (MQ=255)                                                                                              ggtttGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCaaa                                                      <  1:369720/147‑1 (MQ=255)                                                                                               gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTGGCCGCCATCGCCGTCAAGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTCACTGTGGCAGCGGGCGTGCCCGCGAAAACCGGCAAAg                                                     >  2:369720/1‑150 (MQ=255)                                                                                               gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                         <  1:1061094/82‑1 (MQ=255)                                                                                               gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAg                                                                                                                         >  2:1061094/1‑82 (MQ=255)                                                                                               gATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAg                                                     >  1:245581/1‑150 (MQ=255)                                                                                                aTTGAGGTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGc                                                    >  2:1130884/1‑150 (MQ=255)                                                                                                    agaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGt                                                <  2:219797/150‑1 (MQ=255)                                                                                                     gaTGGCAGAGGGTGTGATCGAGGCGATGCTCGAGCGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTg                                               <  2:427868/150‑1 (MQ=255)                                                                                                     gaTGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTg                                               >  2:956609/1‑150 (MQ=255)                                                                                                                   tgATCGAGGCGATCCTCCATGGCATCGAGTCGACGCCAGCGCCGTCACGCCCAGCCTCCAAGGTAGATTACTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGt                                 <  1:208938/150‑1 (MQ=37)                                                                                                                    gATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTg                                <  2:444914/150‑1 (MQ=255)                                                                                                                        gAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACCCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGctt                            >  2:1055099/1‑150 (MQ=255)                                                                                                                          ggCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGcttct                          <  2:428513/150‑1 (MQ=255)                                                                                                                             gATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttcc                       >  2:406183/1‑149 (MQ=255)                                                                                                                              aTACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTcgc                      <  1:746647/150‑1 (MQ=255)                                                                                                                               tACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTGCTTcgcg                     >  1:486617/1‑150 (MQ=255)                                                                                                                                    gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                             <  1:630010/89‑1 (MQ=39)                                                                                                                                    gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCgg                                                                             >  2:630010/1‑89 (MQ=39)                                                                                                                                    gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGcca                                            >  1:1004094/1‑122 (MQ=255)                                                                                                                                    gATGGCGTCGACTCGATGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGcca                                            <  2:1004094/122‑1 (MQ=255)                                                                                                                                    gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTTGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttc                        <  2:58548/142‑1 (MQ=255)                                                                                                                                    gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTTGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCttcttc                        >  1:58548/1‑142 (MQ=255)                                                                                                                                    gATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTaa                >  2:431936/1‑150 (MQ=255)                                                                                                                                         cATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCg           >  1:637175/1‑150 (MQ=255)                                                                                                                                             gAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTcc       <  1:431936/150‑1 (MQ=255)                                                                                                                                               gTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTACACCCGCTCCTa     >  1:872540/1‑150 (MQ=255)                                                                                                                                                tCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTAc    >  2:517488/1‑150 (MQ=255)                                                                                                                                                tCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTAc    <  1:1130884/150‑1 (MQ=255)                                                                                                                                                  gACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTACAg  <  2:637175/150‑1 (MQ=255)                                                                                                                                                  gACGCCATAGACGTCACGCCCCGCCTCCAGGGCATAGTTCTGCTATCGATGCGGGCTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTACAg  <  2:435646/150‑1 (MQ=21)                                                                                                                                                    |                                                                                                                                                    CCCTATACCTGCGCAGCGCTTACGAAACCGGCTACAAAGCTTTGGACTACATGCGTGAAGAAGGCCGTGGCATGTTCTGGTCGAACCTGCATGCGATTGAGATGGCAGAGGGTGTGATCGAGGCGATACTCGATGGCATCGAGTCGACGCCATCGCCGTCACGCCCAGCCTCCAAGGCATAGCGCTGCTATCGAGGCCGACTCAGTACCTGTGGGAGCGGGCGTGCCCGCGAACACCGGCAAAGCCGGTGCCATCCATCGCGGTGGCTTCTTCGCGGGTAAACCCGCTCCTACAG  >  NC_002947/81842‑82136

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