Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,290,063	A→T	5.1%	Q768L (CAG→CTG)	pepN →	aminopeptidase N

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,290,063	0	A	T	5.1%	307.3 / 5.9	99	Q768L (CAG→CTG)	pepN	aminopeptidase N
Reads supporting (aligned to +/- strand):  ref base A (54/39);  new base T (2/3);  total (57/42)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

ATGACCGAGCGCCTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCATGCTGCCGATGGCT  >  NC_002947/2289930‑2290212                                                                                                                                      |                                                                                                                                                      aTGAACGAGCGCTTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTg                                                                                                                                       >  2:2368573/1‑150 (MQ=255)   gACCGAGCGCCTGACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCGCTTCAAGGACAAACCGCTGGTAATGGACCCATGGTGCAGAGTGCTGGACGCCTGCACGCTgaa                                                                                                                                     <  2:1331663/150‑3 (MQ=255)     ccGAGCGCCTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCcgg                                                                                                                                   <  2:20328/150‑1 (MQ=255)       gAGCGCCTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTAGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCcgggc                                                                                                                                 >  1:66465/1‑150 (MQ=255)         gcgcCTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCcgggcgg                                                                                                                               <  1:1747779/150‑1 (MQ=255)          cgcCTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCt                                                                                                                                        >  2:2393311/1‑140 (MQ=255)          cgcCTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCt                                                                                                                                        <  1:2393311/140‑1 (MQ=255)                        ggCGGTACTGGACCTCTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCaa                                                                                                                <  1:2332380/150‑1 (MQ=255)                           ggTACTGGTGAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGc                                                                                                             <  1:953398/150‑1 (MQ=255)                           ggTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGc                                                                                                             >  1:1895960/1‑150 (MQ=255)                            gTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCg                                                                                                            >  2:974525/1‑150 (MQ=255)                             tACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCGGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGGTGCCGGGCGGGCCGGCGCGGGTGAAGGCgg                                                                                                           >  2:680018/1‑149 (MQ=255)                                tGGACAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACGCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGa                                                                                                        >  2:2299999/1‑150 (MQ=255)                                 ggTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGAt                                                                                                       <  2:2245889/150‑1 (MQ=255)                                  gTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCa                                                                                                                                               <  2:905561/109‑1 (MQ=255)                                  gTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCa                                                                                                                                               >  1:905561/1‑109 (MQ=255)                                  gTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAg                                                                                                               <  2:1599895/141‑1 (MQ=255)                                  gTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAg                                                                                                               >  1:1599895/1‑141 (MQ=255)                                  gTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATg                                                                                                      >  1:1525648/1‑150 (MQ=255)                                  gTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATg                                                                                                      >  2:1278394/1‑150 (MQ=255)                                   tCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCTCTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATgc                                                                                                     >  2:948561/1‑150 (MQ=255)                                   tCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATgc                                                                                                     >  1:1583060/1‑150 (MQ=255)                                    cAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATgca                                                                                                    <  2:201238/150‑1 (MQ=255)                                        tCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGTGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCAc                                                                                                >  2:959624/1‑150 (MQ=255)                                          gCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCAccc                                                                                              >  1:1608522/1‑150 (MQ=255)                                                 gAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTAATGCAGCACCCGGCCTTc                                                                                       <  1:596770/150‑1 (MQ=255)                                                  aGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGAACGTACTACTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCa                                                                                      >  1:1240929/1‑150 (MQ=255)                                                   ggCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCAc                                                                                     <  1:1295002/150‑1 (MQ=255)                                                   ggCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGGGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCAc                                                                                     >  1:1379387/1‑150 (MQ=255)                                                    gCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCAcc                                                                                    >  1:1826748/1‑150 (MQ=255)                                                     ccGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCAccc                                                                                   >  1:1827020/1‑150 (MQ=255)                                                       gAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAg                                                                                                               >  2:1363780/1‑120 (MQ=255)                                                       gAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAg                                                                                                               <  1:1363780/120‑1 (MQ=255)                                                       gAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTg                                                                                 >  1:2366926/1‑150 (MQ=255)                                                         gcgcgcCCAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCACTCAATTAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTgaa                                                                               <  1:746625/150‑1 (MQ=255)                                                              ccAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGa                                                                                                        >  2:1734904/1‑120 (MQ=255)                                                              ccAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGa                                                                                                        <  1:1734904/120‑1 (MQ=255)                                                              ccAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAAcc                                                                          <  2:1895960/150‑1 (MQ=255)                                                               cAAAGCGCTGGAGGCCTTTGCCGAGCACTTCATGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAAccc                                                                         >  2:1345678/1‑150 (MQ=255)                                                                      cTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAg                                                                  <  1:1692880/150‑1 (MQ=255)                                                                       tGGAGGCCTTCGCCGAACCTTCCAAATATTACCCTCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGAACAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAgg                                                                 <  1:2368573/150‑1 (MQ=255)                                                                          aGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGgtgc                                                              <  2:1827020/150‑1 (MQ=255)                                                                          aGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGgtgc                                                              <  1:1692100/150‑1 (MQ=255)                                                                              cTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGgtgcgtgc                                                          >  2:1925483/1‑150 (MQ=255)                                                                               tttGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATgca                                                                                                    >  2:2177741/1‑107 (MQ=255)                                                                               tttGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATgca                                                                                                    <  1:2177741/107‑1 (MQ=255)                                                                               tttGCCGAGCACTTCAAGGACAACCAGCTGCTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGtgcgtgcg                                                         >  1:1880543/1‑150 (MQ=255)                                                                                ttGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTgcgc                                                        <  2:1583060/150‑1 (MQ=255)                                                                                ttGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTgcgc                                                        <  1:1928955/150‑1 (MQ=255)                                                                                   ccGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGa                                                     >  2:1432913/1‑150 (MQ=255)                                                                                     gAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCAAGCTGCAGGGCGGGCTGGCGCGGGTAAAGGAGTTGATGCAGCACCCGGTCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATc                                                   <  2:2102933/150‑1 (MQ=255)                                                                                        cACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGgc                                                <  1:1407931/150‑1 (MQ=255)                                                                                          cTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCggg                                                                                                                    >  2:1894131/1‑80 (MQ=255)                                                                                          cTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCggg                                                                                                                    <  1:1894131/80‑1 (MQ=255)                                                                                          cTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGgcgc                                              >  2:1335712/1‑150 (MQ=255)                                                                                           ttCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAg                                                                  <  1:1399323/129‑1 (MQ=255)                                                                                           ttCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAg                                                                  >  2:1399323/1‑129 (MQ=255)                                                                                           ttCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGgtgc                                                              >  2:1765299/1‑133 (MQ=255)                                                                                           ttCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGgtgc                                                              <  1:1765299/133‑1 (MQ=255)                                                                                           ttCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGcc                                             >  2:921679/1‑150 (MQ=255)                                                                                           ttCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGcc                                             >  2:88404/1‑150 (MQ=255)                                                                                            tCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCt                                            >  1:815354/1‑150 (MQ=255)                                                                                               aGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTg                                         >  1:2280204/1‑150 (MQ=255)                                                                                                 gACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGcc                                       >  2:995944/1‑150 (MQ=255)                                                                                                 gACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGcc                                       >  2:70034/1‑150 (MQ=255)                                                                                                 gACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGcc                                       <  1:88404/150‑1 (MQ=255)                                                                                                   cAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTg                                                      >  2:912218/1‑133 (MQ=255)                                                                                                   cAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTg                                                      <  1:912218/133‑1 (MQ=255)                                                                                                    aaCCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGa                                                                                                        <  2:2262214/82‑1 (MQ=255)                                                                                                    aaCCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGa                                                                                                        >  1:2262214/1‑82 (MQ=255)                                                                                                       ccGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCAACCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAg                                 >  1:880263/1‑150 (MQ=255)                                                                                                          cTGGTCATGGACCAATGGTTCAGGGAGCGGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAAc                              >  1:917646/1‑150 (MQ=255)                                                                                                          cTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGTTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTTCGTGCGCTGATCGGCGCCTTTGCCTGGCAGAAc                              >  2:1331363/1‑150 (MQ=255)                                                                                                             gTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCTAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTg                           <  2:778404/150‑1 (MQ=255)                                                                                                              tCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTgg                          >  2:2427247/1‑150 (MQ=255)                                                                                                              tCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTgg                          <  1:2103765/150‑1 (MQ=255)                                                                                                                aTGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTc                        >  1:1703076/1‑150 (MQ=255)                                                                                                                aTGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTc                        <  1:2299999/150‑1 (MQ=255)                                                                                                                  ggACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCgg                                                                                                                     >  2:1805263/1‑55 (MQ=255)                                                                                                                  ggACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCgg                                                                                                                     <  1:1805263/55‑1 (MQ=255)                                                                                                                    aCCAATGGTTCAGCGTGCAGGCGGCCAGCACGGTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACt                    >  2:1344224/1‑150 (MQ=255)                                                                                                                    aCCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACt                    >  1:1978521/1‑150 (MQ=255)                                                                                                                      cAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACtgc                  >  1:1409080/1‑148 (MQ=255)                                                                                                                       aaTGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTcc                 <  1:959624/150‑1 (MQ=255)                                                                                                                        aTGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCa                >  1:351158/1‑150 (MQ=255)                                                                                                                        aTGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCa                >  2:1976176/1‑150 (MQ=255)                                                                                                                           gTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGAt                                                                                                       >  1:539830/1‑60 (MQ=255)                                                                                                                           gTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGAt                                                                                                       <  2:539830/60‑1 (MQ=255)                                                                                                                           gTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAActgg                                                                  >  1:1641996/1‑94 (MQ=255)                                                                                                                           gTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAActgg                                                                  <  2:1641996/97‑4 (MQ=255)                                                                                                                            ttCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATGGTCGCCTTTGCCGTGCAGAACCTGGTCAACTTCCAtgct            >  1:2342978/1‑150 (MQ=255)                                                                                                                               aGCGTGCTGGCCGCCTGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGgcgc                                              <  2:1582762/113‑1 (MQ=255)                                                                                                                               aGCGTGCTGGCCGCCTGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGgcgc                                              >  1:1582762/1‑113 (MQ=255)                                                                                                                                gCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCATGCTGCCg        >  1:1569491/1‑150 (MQ=255)                                                                                                                                 cGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCATGCTGCCGa       <  1:941833/150‑1 (MQ=255)                                                                                                                                 cGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCATGCTGCCGa       <  1:941823/150‑1 (MQ=255)                                                                                                                                  gTGCTGGCCGCCTGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCAcc                                                                                    >  1:2047390/1‑72 (MQ=38)                                                                                                                                  gTGCTGGCCGCCTGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCAcc                                                                                    <  2:2047390/72‑1 (MQ=38)                                                                                                                                   tGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCATGCTGCCGATg     >  1:1728011/1‑150 (MQ=255)                                                                                                                                      aGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCATGCTGCCGATGGCt  >  2:1917771/1‑150 (MQ=255)                                                                                                                                      |                                                                                                                                                      ATGACCGAGCGCCTTACCGCGCTGGCGGTACTGGTCAACTCGCCGTTCGAGGCCGAGCGCGCCAAAGCGCTGGAGGCCTTTGCCGAGCACTTCAAGGACAACCCGCTGGTCATGGACCAATGGTTCAGCGTGCAGGCGGCCAGCACGCTGCCGGGCGGGCTGGCGCGGGTCAAGGCGTTGATGCAGCACCCGGCCTTCACCCTGAAGAACCCGAACAAGGTGCGTGCGCTGATCGGCGCCTTTGCCGGGCAGAACCTGGTCAACTTCCATGCTGCCGATGGCT  >  NC_002947/2289930‑2290212

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