Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,015,317	T→A	6.7%	E57D (GAA→GAT)	PP_1794 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,015,317	0	T	A	6.7%	222.7 / 3.7	74	E57D (GAA→GAT)	PP_1794	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/34);  new base A (3/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

AGTTACCTTCCAGGTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCCAGCAGATCGCTG  >  NC_002947/2015169‑2015461                                                                                                                                                      |                                                                                                                                                 aGTTACCTTCCAGGTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGAt                                                                                                                                                  >  2:298529/1‑149 (MQ=255)     aCCTTCCAGGTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTcc                                                                                                                                               <  8:17131/148‑1 (MQ=255)        ttCCAGGTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCgg                                                                                                                                             >  3:196938/1‑147 (MQ=255)        ttCCAGGTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCgg                                                                                                                                             >  8:276807/1‑147 (MQ=255)           cAGGTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCGGGATTCCGGCAGCt                                                                                                                                        <  3:217291/149‑1 (MQ=255)             ggTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGAATCCGGCactgcc                                                                                                                                      <  6:35196/149‑6 (MQ=255)               ttGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAAc                                                                                                                                    >  5:418703/1‑149 (MQ=255)                   cTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGgtg                                                                                                                                <  5:35352/149‑1 (MQ=255)                      ccAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTg                                                                                                                             >  8:247220/1‑149 (MQ=255)                      ccAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTg                                                                                                                             >  1:71811/1‑149 (MQ=255)                       cAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGc                                                                                                                            >  4:169294/1‑149 (MQ=255)                              gCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGa                                                                                                                      <  6:418703/148‑1 (MQ=255)                                aTCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGtt                                                                                                                              <  2:164471/138‑1 (MQ=255)                                aTCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGtt                                                                                                                              >  1:164471/1‑138 (MQ=255)                                aTCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGc                                                                                                                   >  5:413965/1‑149 (MQ=255)                                   cGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAg                                                                                                                <  7:247220/149‑1 (MQ=255)                                   cGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAg                                                                                                                <  1:298529/149‑1 (MQ=255)                                   cGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAg                                                                                                                <  2:71811/149‑1 (MQ=255)                                    gACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGc                                                                                                               <  7:276807/149‑1 (MQ=255)                                     aCTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAg                                                                                                                                          <  4:7965/121‑1 (MQ=255)                                     aCTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAg                                                                                                                                          >  3:7965/1‑121 (MQ=255)                                         cagcaCGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGAtt                                                                                                           <  7:111343/148‑1 (MQ=255)                                         cagcaCGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGAt                                                                                                            <  5:140963/147‑1 (MQ=255)                                            caCGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCa                                                                                                       <  3:169294/149‑1 (MQ=255)                                               gTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACgg                                                                                                    >  5:30544/1‑149 (MQ=255)                                                  ggCATTCACGCCAGCGGCTTCAGAAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGTCAGTGCCGGATTCAGGAAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGtc                                                                                                  <  4:327308/148‑1 (MQ=255)                                                         aCGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCg                                                                                                              <  8:187145/129‑1 (MQ=255)                                                         aCGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCg                                                                                                              >  7:187145/1‑129 (MQ=255)                                                           gCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTAt                                                                                        >  6:15066/1‑149 (MQ=255)                                                                   cTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGtt                                                                                 >  8:212042/1‑148 (MQ=255)                                                                         cagcagcCAGGTCGACGATGACCTGGGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTcagcag                                                                           >  1:3552/1‑149 (MQ=255)                                                                         cagcagcCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCa                                                                                     >  3:215536/1‑138 (MQ=255)                                                                         cagcagcCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCa                                                                                     <  4:215536/138‑1 (MQ=255)                                                                         cagcagcCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAgg                                                                          >  1:12908/1‑149 (MQ=255)                                                                              gcCAGGTCCACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGtt                                                                     <  2:396934/149‑1 (MQ=255)                                                                               ccAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTctct                                                                                               <  8:116149/122‑1 (MQ=255)                                                                               ccAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTctct                                                                                               >  7:116149/1‑122 (MQ=255)                                                                                    tCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGc                                                               >  5:398611/1‑149 (MQ=255)                                                                                      gacgaTGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCt                                                             <  6:413965/149‑1 (MQ=255)                                                                                         gaTGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGa                                                          <  6:375011/149‑1 (MQ=255)                                                                                         gaTGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGa                                                          <  7:212042/149‑1 (MQ=255)                                                                                             aCCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGt                                                      >  1:403393/1‑149 (MQ=255)                                                                                             aCCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGt                                                      <  2:403393/149‑1 (MQ=255)                                                                                               cT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTgg                                                    >  3:8863/1‑149 (MQ=255)                                                                                                   gTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGa                                                   <  2:12908/147‑1 (MQ=255)                                                                                                    tCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTc                                                >  4:350894/1‑149 (MQ=255)                                                                                                    tCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCACCTGAACGGTGTTGCCAATGATGCCAGCGATTTCCACGGTCGCTACGTTATCCAGGTACAGAAGGCTGTTGATAACCTCGATAGTGAAGTc                                                >  7:252097/1‑149 (MQ=255)                                                                                                      gCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTg                                                     >  7:52625/1‑142 (MQ=255)                                                                                                      gCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTg                                                     <  8:52625/142‑1 (MQ=255)                                                                                                      gCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCg                                               >  8:21069/1‑148 (MQ=255)                                                                                                        cTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGgag                                            >  4:133122/1‑149 (MQ=255)                                                                                                         tttAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGga                                             >  2:296059/1‑147 (MQ=255)                                                                                                           tAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGgaga                                           <  7:21069/147‑1 (MQ=255)                                                                                                           tAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGgaga                                           <  4:196938/147‑1 (MQ=255)                                                                                                              cgcCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTg                                      <  6:398611/149‑1 (MQ=255)                                                                                                                   tCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCt                                 >  5:146066/1‑149 (MQ=255)                                                                                                                       cAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTcg                             <  2:3552/149‑1 (MQ=255)                                                                                                                         gTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGAATTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGAGAGTGGAGTCGGAGAGGTTGGCGCTTTcggc                           >  3:272699/1‑149 (MQ=255)                                                                                                                             cTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATa                                                                 >  2:53965/1‑107 (MQ=255)                                                                                                                             cTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGt                        >  6:366861/1‑148 (MQ=255)                                                                                                                             cTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGACTCTACGTTATCCAGGTTCAGCAGGCTGTTGATa                                                                 <  1:53965/107‑1 (MQ=255)                                                                                                                                gAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCa                      <  5:366861/147‑1 (MQ=255)                                                                                                                                 aaCAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAgg                    <  4:8863/148‑1 (MQ=255)                                                                                                                                   cAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGtt                 >  6:233942/1‑149 (MQ=255)                                                                                                                                    aaTGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTg                >  7:333188/1‑149 (MQ=255)                                                                                                                                        gttcagctGCCGGAATCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCcag            >  7:92112/8‑149 (MQ=255)                                                                                                                                          ggCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGGCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGcc              >  1:333100/1‑145 (MQ=255)                                                                                                                                          ggCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGcc              <  2:333100/145‑1 (MQ=255)                                                                                                                                           gCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAggg                   >  2:54901/1‑139 (MQ=255)                                                                                                                                           gCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAggg                   <  1:54901/139‑1 (MQ=255)                                                                                                                                             gctGCCGGAATCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCcag            >  5:172117/3‑144 (MQ=255)                                                                                                                                             gctGCCGGAATCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCCAGCAGa        >  2:286803/3‑148 (MQ=255)                                                                                                                                                 ccGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCCAGCAGATCgct   >  7:136723/1‑149 (MQ=255)                                                                                                                                                 ccGGAATCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGACTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCCAGCAGATCg     <  1:286803/147‑1 (MQ=255)                                                                                                                                                  cGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCCAGCAGATCgctg  <  5:15066/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                 AGTTACCTTCCAGGTTGACTACCAGGTTGGCATCCGACTGCAGCACGTAGGCATTAACGCCAGCGGCTTCAGCAGCAGCCAGGTCGACGATGACCT‑GGTCGCCTTTAGCGCCATCGACAGTACCTTGAACAATGCTGGCAGTGCCGGATTCCGGCAGCTGAACGGTGTTGCCATTGATGCCAGCGATTTCCACGGTCTCTACGTTATCCAGGTTCAGCAGGCTGTTGATAGCCTCGATAGTGGAGTCGGAGAGGTTGGCGCTTTCGGCGGTCAGGGTTGCCAGCAGATCGCTG  >  NC_002947/2015169‑2015461

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