Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,106,043	T→G	6.4%	D139A (GAC→GCC)	PP_5355 ←	sodium/proton antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,106,043	0	T	G	6.4%	272.3 / 3.8	88	D139A (GAC→GCC)	PP_5355	sodium/proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (42/40);  new base G (4/2);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACTGTGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAA  >  NC_002947/6105897‑6106191                                                                                                                                                   |                                                                                                                                                         cGACTGTGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCg                                                                                                                                                        >  7:291827/1‑149 (MQ=255)   aCTGTGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGtt                                                                                                                                                      <  5:42674/149‑1 (MQ=255)       tgGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAgc                                                                                                                                                  <  5:169830/149‑1 (MQ=255)        ggAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGATGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAgcg                                                                                                                                                 >  7:160702/1‑149 (MQ=255)         gagaACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAgcgc                                                                                                                                                >  5:393659/1‑149 (MQ=255)            aaCACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCtt                                                                                                                                             >  2:399147/1‑149 (MQ=255)              cacCTCGCCCGAGCTGACCGTGCCGTTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTg                                                                                                                                            >  4:385823/1‑148 (MQ=255)                aatCGCCCGATCTTACCGTGCCGCTTATCTCTATACCCGACAGAAGGGCAATGATACCGATGCCAAGAAAGTCATCTACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCTTTCTTCATCTTCAAGTCGTTGAGAGCCTTGACg                                                                                                                                         <  5:188161/147‑1 (MQ=255)                 cTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACga                                                                                                                                        >  2:67963/1‑149 (MQ=255)                 cTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAAAGATACCGATGCCAAGAATGTCCTCGACGATCAGCCCGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACga                                                                                                                                        >  5:362200/1‑149 (MQ=255)                           cTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑Atggtggt                                                                                                                             >  3:414151/1‑149 (MQ=255)                              aCCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCACGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTgga                                                                                                                          <  3:8239/149‑1 (MQ=255)                              aCCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑Atggtggtg                                                                                                                            >  5:365630/1‑147 (MQ=255)                              aCCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTgga                                                                                                                          <  2:422240/149‑1 (MQ=255)                                      gCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGa                                                                                                                                           <  4:146887/125‑1 (MQ=255)                                      gCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGa                                                                                                                                           >  3:146887/1‑125 (MQ=255)                                       cTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑Gc                                                                                                              <  6:393659/149‑1 (MQ=255)                                       cTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑Gc                                                                                                              <  7:339912/149‑1 (MQ=255)                                       cTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACCATG‑ATGGTGGTGGAGGAAATG‑‑‑Gc                                                                                                              >  7:185863/1‑149 (MQ=255)                                             gcgATCCCCGACAGCATGGCAATGATACCGATGCCAAGAATGTCCTCGACGATAAGCACTCCGAATATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATTGTGGTGGAGGAAATG‑‑‑GCCAGGAt                                                                                                        <  4:50313/149‑1 (MQ=255)                                               gATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATgg                                                                                                      >  2:110850/1‑149 (MQ=255)                                                 tACCCGCCAGCATTGCAATGATACCTATCCCAAGCCAGTCCTCGACGATCAGCACGCCGAAGATCATCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGgc                                                                                                     <  6:136466/148‑1 (MQ=255)                                                 tACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGgcg                                                                                                    <  4:414151/149‑1 (MQ=255)                                                 tACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGgcg                                                                                                    <  8:291827/149‑1 (MQ=255)                                                   cccGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGgcgc                                                                                                   <  3:427180/148‑1 (MQ=255)                                                    ccGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGcc                                                                                                  >  3:369993/1‑148 (MQ=255)                                                         agcagGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGaa                                                                                            >  6:324409/1‑149 (MQ=255)                                                         agcagGGCAATGATACCCATGCCAAGAATGTCCTCGACGAACAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTACTTCATCATCAAGTCGATCAGCGCCTACACGATG‑ATCGGCGGCGAGGAAATC‑‑‑GACAGGATCGCGCCCAgaa                                                                                             >  5:248237/1‑146 (MQ=255)                                                             ttgCAATGATACCGATGCCAAGAATGTCCTCGACTATCAGCACGCCGAAGATCAGCTGGGGAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGTCGATG‑ATGTTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGt                                                                                        <  4:130379/147‑1 (MQ=255)                                                             gggCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGt                                                                                        <  5:444922/149‑1 (MQ=255)                                                             gggCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGt                                                                                        <  5:54104/149‑1 (MQ=255)                                                              ggCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGt                                                                                        >  8:61621/1‑148 (MQ=255)                                                                    gATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTgga                                                                                                                          <  6:104206/111‑1 (MQ=255)                                                                    gATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTgga                                                                                                                          >  5:104206/1‑111 (MQ=255)                                                                      tACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCTAt                                                                                <  2:37900/148‑1 (MQ=255)                                                                      tACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCAt                                                                                <  4:218717/148‑1 (MQ=255)                                                                      tACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATg                                                                               >  2:140838/1‑149 (MQ=255)                                                                      tACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATg                                                                               <  6:365630/149‑1 (MQ=255)                                                                      tACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATg                                                                               <  1:399147/149‑1 (MQ=255)                                                                             gCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATTTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTcc                                                                        >  4:410122/1‑149 (MQ=255)                                                                             gCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTc                                                                         >  4:424936/1‑148 (MQ=255)                                                                              ccAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCa                                                                       <  3:427829/149‑1 (MQ=255)                                                                               cAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAg                                                                      >  6:126621/1‑149 (MQ=255)                                                                                    aTGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCaaa                                                                 >  3:160985/1‑149 (MQ=255)                                                                                          tCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGc                                                                     >  6:131192/1‑139 (MQ=255)                                                                                          tCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGc                                                                     <  5:131192/139‑1 (MQ=255)                                                                                           cgacgaTCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCgg                                                          >  5:166649/1‑149 (MQ=255)                                                                                                aTCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAt                                                     <  5:161520/149‑1 (MQ=255)                                                                                                   aGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtct                                                   <  4:160985/148‑1 (MQ=255)                                                                                                    gCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtctc                                                  >  8:61824/1‑148 (MQ=255)                                                                                                    gCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAAtctc                                                  >  4:355108/1‑148 (MQ=255)                                                                                                     cACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCaa                                                <  1:81194/149‑1 (MQ=255)                                                                                                     cACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCa                                                 >  3:167366/1‑148 (MQ=255)                                                                                                       cGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCaaaa                                              <  3:355108/149‑1 (MQ=255)                                                                                                        gCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAAc                                             <  3:211784/149‑1 (MQ=255)                                                                                                           gAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGAAGAATGGTGTTGGAGGAAATG‑‑‑GCCAGGATGGCGCCAAGGAAAAGTGAATCAATGGTGGTCCAGCCAAACCAGCGGCCAAACTCAAAAcc                                            >  1:125122/1‑148 (MQ=255)                                                                                                                tCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAg                                     <  1:67963/149‑1 (MQ=255)                                                                                                                  aGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATc                                        >  4:82110/1‑144 (MQ=255)                                                                                                                  aGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATc                                        >  6:147123/1‑144 (MQ=255)                                                                                                                    cTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATGGGAGCCAGAATGGCGCCCAGGAAAAGTCAAACCATGTTCATCCAGCCAAACCAGCGGCCAATAACAAAACCAATCCAg                                     >  1:136644/1‑148 (MQ=255)                                                                                                                    cTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGatca                                 >  4:379486/1‑149 (MQ=255)                                                                                                                       ggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGATCGCCTTGACGATG‑ATGGTTGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGatcatc                               <  7:172344/148‑1 (MQ=255)                                                                                                                       ggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGatcatc                               <  5:126621/148‑1 (MQ=255)                                                                                                                       ggCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGAtcatca                              <  1:110850/149‑1 (MQ=255)                                                                                                                       ggCAAAGCGCTCGTTCTTCATCATCCAGGCGTGGAGCGCCTTGACGAGC‑ATGGTGGGGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCCGTtcatca                              <  6:362200/149‑1 (MQ=255)                                                                                                                         cAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGatc                                  <  4:182410/143‑1 (MQ=255)                                                                                                                         cAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGatc                                  >  3:182410/1‑143 (MQ=255)                                                                                                                          aaaGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGa                           >  3:223603/1‑149 (MQ=255)                                                                                                                               gcTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCGAAACCGATCCAGATCATCAGGACGAtt                      >  5:21474/1‑149 (MQ=255)                                                                                                                                cTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCa                                      >  3:281619/1‑132 (MQ=255)                                                                                                                                cTCGTTCTTCATCTTCAAGTCGTTCAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGAttt                     >  7:269873/1‑149 (MQ=255)                                                                                                                                cTCGGTCTTCATCTTCAAGTCGTTGCGCGCCTTGACGATG‑ATTGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCa                                      <  4:281619/132‑1 (MQ=255)                                                                                                                                 tCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTc                    >  3:331156/1‑149 (MQ=255)                                                                                                                                   gTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTcc                   <  5:324409/148‑1 (MQ=255)                                                                                                                                   gTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCa                  <  4:331156/149‑1 (MQ=255)                                                                                                                                   gTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCa                  <  1:140838/149‑1 (MQ=255)                                                                                                                                    ttcttcATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAg                 <  7:61824/149‑1 (MQ=255)                                                                                                                                      catcatcgtcaAGGCGCTCAACGACTTGAAGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGa                           <  4:262066/135‑1 (MQ=38)                                                                                                                                      catcatcgtcaAGGCGCTCAACGACTTGAAGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGa                           >  3:262066/3‑137 (MQ=38)                                                                                                                                        tcatcgtcaAGGCGCTCAACGACTTGAAGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAAc             >  5:351971/1‑149 (MQ=255)                                                                                                                                        tcatcgtcaAGGCGCTCAACGACTTGAAGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCAAGGAAc             >  1:322605/1‑149 (MQ=255)                                                                                                                                           tcgtcaAGGCGCTCAACGACTTGAAGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGCCGATTTCCAGGAACGCt          >  2:63075/4‑149 (MQ=255)                                                                                                                                             ttcaAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGc        >  8:214291/1‑149 (MQ=255)                                                                                                                                               caAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCt          <  6:166649/145‑1 (MQ=255)                                                                                                                                               caAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGa      >  7:247885/1‑149 (MQ=255)                                                                                                                                                aaGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGAt     <  8:45583/149‑1 (MQ=255)                                                                                                                                                 aGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGAt     >  3:117597/1‑148 (MQ=255)                                                                                                                                                  gTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATg    >  1:212790/1‑148 (MQ=255)                                                                                                                                                   tCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGaa  >  1:31401/1‑149 (MQ=255)                                                                                                                                                   tCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGaa  <  8:269873/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                         CGACTGTGGAGAACACCTCGCCCGAGCTGACCGTGCCGCTGACCGCGATACCCGACAGCAGGGCAATGATACCGATGCCAAGAATGTCCTCGACGATCAGCACGCCGAAGATCAGCTGGGCAAAGCGCTCGTTCTTCATCTTCAAGTCGTTGAGCGCCTTGACGATG‑ATGGTGGTGGAGGAAATG‑‑‑GCCAGGATGGCGCCCAGGAACAGTGAATCCATGGTGCTCCAGCCAAACCAGCGGCCAATCTCAAAACCGATCCAGATCATCAGGACGATTTCCAGGAACGCTGCGATGAA  >  NC_002947/6105897‑6106191

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