Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,049,646	A→T	5.1%	K718M (AAG→ATG) ‡	quiA →	quinate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,049,646	0	A	T	5.1%	238.2 / 3.6	79	K718M (AAG→ATG) ‡	quiA	quinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

CTGTCGGTCGCCGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGG  >  NC_002947/4049498‑4049793                                                                                                                                                     |                                                                                                                                                    cTGTCGGTCGCCGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAg                                                                                                                                                    >  2:1127743/1‑150 (MQ=255)         cgCCGGCATGCCCGGCCAGGCACCGCCGTGTGGCACGCTGACCGCCAGCGTGAGGAGGGCCCGCCAGGTGGCCTGGCAGGGGCCGGTAGTCACCGTGGAAGACACCGTCCCGCTCGGCAGCCGCATGCACCTGCCGAGCTTGATCGGcag                                                                                                                                            <  2:256676/150‑3 (MQ=255)          gCCGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTg                                                                                                                                           >  1:298336/1‑150 (MQ=255)            cGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGcc                                                                                                                                         <  1:812304/150‑1 (MQ=255)            cGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGAGCAAGATCGGCCTGcc                                                                                                                                         >  1:222090/1‑150 (MQ=255)              gCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCaa                                                                                                                                       <  1:1219534/150‑1 (MQ=255)                 ttCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCaa                                                                                                                                                     <  1:483994/133‑1 (MQ=255)                 ttCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCaa                                                                                                                                                     >  2:483994/1‑133 (MQ=255)                     aTGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTgg                                                                                                                                >  1:965188/1‑150 (MQ=255)                     aTGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTgg                                                                                                                                <  2:23204/150‑1 (MQ=255)                            gCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCAcc                                                                                                                         >  1:1885964/1‑150 (MQ=255)                                cgccgTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGCTCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGt                                                                                                                     >  2:1680660/1‑150 (MQ=255)                                cgccgTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGt                                                                                                                     >  2:1780421/1‑150 (MQ=255)                                 gccgTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTc                                                                                                                    <  2:1442585/150‑1 (MQ=255)                                   cgTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGa                                                                                                                  <  2:1885964/150‑1 (MQ=255)                                     tttGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTCTCGAcc                                                                                                                >  2:360912/1‑150 (MQ=255)                                            cccTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAggcg                                                                                                         <  2:1530607/150‑1 (MQ=255)                                             ccTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAggcgg                                                                                                        >  2:489859/1‑150 (MQ=255)                                               tGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGcc                                                                                                      >  1:1406713/1‑150 (MQ=255)                                                gACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCt                                                                                                     >  1:929643/1‑150 (MQ=255)                                                   cGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGgt                                                                                                  >  2:1742198/1‑150 (MQ=255)                                                     ccATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCt                                                                                                     <  2:125779/145‑1 (MQ=255)                                                     ccATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCt                                                                                                     >  1:125779/1‑145 (MQ=255)                                                     ccATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGgtgt                                                                                                >  2:1558572/1‑150 (MQ=255)                                                     ccATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGgtgt                                                                                                <  1:1724985/150‑1 (MQ=255)                                                      cATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGtt                                                                                               >  2:714866/1‑150 (MQ=255)                                                               gAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAgg                                                                                      >  1:293643/1‑150 (MQ=255)                                                               gAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAgg                                                                                      >  1:1837890/1‑150 (MQ=255)                                                                  gACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTTGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTAc                                                                                   >  2:1179352/1‑150 (MQ=255)                                                                   aCCTGCCAGGCCACCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTAcc                                                                                  >  1:1678749/1‑150 (MQ=255)                                                                   aCCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCCGGTGTTCATCGCAGGTAcc                                                                                  <  1:701590/150‑1 (MQ=255)                                                                     ccGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCAc                                                                                                                          <  1:667100/108‑1 (MQ=255)                                                                     ccGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCAc                                                                                                                          >  2:667100/1‑108 (MQ=255)                                                                     ccGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCa                                                                                >  1:1055035/1‑150 (MQ=255)                                                                     ccGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCa                                                                                <  2:929643/150‑1 (MQ=255)                                                                      cGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAg                                                                               >  1:302089/1‑150 (MQ=255)                                                                      cGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAg                                                                               >  1:1914085/1‑150 (MQ=255)                                                                      cGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAg                                                                               <  2:1914085/150‑1 (MQ=255)                                                                         cAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGAc                                                                            <  2:1406713/150‑1 (MQ=255)                                                                             gggCCGGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCGCGATCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCt                                                                        <  2:1684914/149‑1 (MQ=255)                                                                              ggCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTa                                                                       <  1:1779066/150‑1 (MQ=255)                                                                               gCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAg                                                                               <  2:1054286/141‑1 (MQ=255)                                                                               gCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAg                                                                               >  1:1054286/1‑141 (MQ=255)                                                                                         tGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCt                                                            <  2:574116/150‑1 (MQ=255)                                                                                          gCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTa                                                           <  2:1678749/150‑1 (MQ=255)                                                                                             ggTAGGCACCGTTGAAGACGCCGGCCCGCTCGGCCTCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGa                                                        <  2:310796/150‑1 (MQ=255)                                                                                                  gCACCGTTGAAGACACCGGCCCGCTCCGAATCCGCATGCATCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGAcagca                                                   <  1:1866089/150‑1 (MQ=255)                                                                                                   cACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAg                                                                               >  1:1498522/1‑121 (MQ=255)                                                                                                   cACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAg                                                                               <  2:1498522/121‑1 (MQ=255)                                                                                                    aCCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGa                                                                             >  1:1511670/1‑122 (MQ=255)                                                                                                    aCCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGa                                                                             <  2:1511670/122‑1 (MQ=255)                                                                                                            cgACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTaa                                         <  2:298336/149‑1 (MQ=255)                                                                                                            aGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTaa                                         <  1:1127743/150‑1 (MQ=255)                                                                                                             gACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCCACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAAc                                        >  1:1749990/1‑150 (MQ=255)                                                                                                                acCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGaa                                     >  1:1110971/1‑150 (MQ=255)                                                                                                                acCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGaa                                     >  1:1283025/1‑150 (MQ=255)                                                                                                                acCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTAATCGCAGGTACCCAGCAATTCTACTTGCTCGCCTACGACAGCAGCAACGGTAACGaa                                     >  2:449123/1‑150 (MQ=255)                                                                                                                acCGGCCCGCCCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCACAGCCGCCTGGTGTTCATCGCAGGCACCCAGGACTTCTACCTGCGCGCCTCCGACAGCCGCAACGGTAACGaa                                     >  2:1309735/1‑150 (MQ=255)                                                                                                                    gCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCt                                 <  2:513609/150‑1 (MQ=255)                                                                                                                     cccGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTg                                <  2:683914/150‑1 (MQ=255)                                                                                                                     cccGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTg                                <  2:1110971/150‑1 (MQ=255)                                                                                                                       cGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACGCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGa                              <  2:1033565/150‑1 (MQ=255)                                                                                                                       cGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCa                                                                                             >  1:77429/1‑87 (MQ=255)                                                                                                                       cGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCa                                                                                             <  2:77429/87‑1 (MQ=255)                                                                                                                             gCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTc                        <  1:714866/150‑1 (MQ=255)                                                                                                                                 ccGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCt                    >  2:134701/1‑150 (MQ=255)                                                                                                                                   gCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTc                                                                         <  2:924168/95‑1 (MQ=255)                                                                                                                                   gCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTc                                                                         >  1:924168/1‑95 (MQ=255)                                                                                                                                    cATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCt                                 >  2:511153/1‑134 (MQ=255)                                                                                                                                    cATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTCCGCGCCTACGACAGCAGCAACGGTAACGAAATCTCGAAGGCTCACCTccc                 >  2:105562/1‑150 (MQ=255)                                                                                                                                    cATGCACCGGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCt                                 <  1:511153/134‑1 (MQ=255)                                                                                                                                     aTGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTAc                                                                      <  2:773374/96‑1 (MQ=255)                                                                                                                                     aTGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTAc                                                                      >  1:773374/1‑96 (MQ=255)                                                                                                                                         aCCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGCCTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTgg            <  1:1742198/150‑1 (MQ=255)                                                                                                                                         aCCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTgg            >  1:630922/1‑150 (MQ=255)                                                                                                                                           cTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGc          >  2:702066/1‑150 (MQ=255)                                                                                                                                           aggCCGATCTTGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTgcgc                                                                >  1:1002455/3‑96 (MQ=255)                                                                                                                                           aggCCGATCTTGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGc          >  1:898868/3‑150 (MQ=255)                                                                                                                                           aggCCGATCTTGATCGGCCTGCCAACCCTTGGCGGCACCATGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTgcgc                                                                <  2:1002455/94‑1 (MQ=255)                                                                                                                                                gATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCa     >  1:1705150/1‑150 (MQ=255)                                                                                                                                                   cAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAggg  <  2:1705150/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CTGTCGGTCGCCGGCATTCCATGCCAGGCACCGCCGTTTGGCACCCTGACCGCCATCGACATGAAGACCCGCCAGGTGGCCTGGCAGGTGCCGGTAGGCACCGTTGAAGACACCGGCCCGCTCGGCATCCGCATGCACCTGCCGATCAAGATCGGCCTGCCAACCCTTGGCGGCACCCTGTCGACCCAAGGCGGCCTGGTGTTCATCGCAGGTACCCAGGACTTCTACCTGCGCGCCTACGACAGCAGCAACGGTAACGAAATCTGGAAGGCTCGCCTGCCCGTGGGCAGCCAGGG  >  NC_002947/4049498‑4049793

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