Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,823,170 C→T 5.2% P17L (CCG→CTG)  yggW → coproporphyrinogen/heterocyclic compound oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,823,1700CT 5.2% 226.0 / 2.1 77P17L (CCG→CTG) yggWcoproporphyrinogen/heterocyclic compound oxidase
Reads supporting (aligned to +/- strand):  ref base C (32/41);  new base T (2/2);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.95e-01

AGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAGCTGGCTGC  >  NC_002947/5823021‑5823314
                                                                                                                                                     |                                                                                                                                                
aGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAANTGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTgcc                                                                                                                                                  <  1:2433649/150‑1 (MQ=255)
   gcAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTgccgcc                                                                                                                                               >  1:211123/1‑150 (MQ=255)
     aaCTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCcgccgc                                                                                                                                             >  1:2382408/1‑150 (MQ=255)
           tCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGgcgc                                                                                                                                       >  1:465458/1‑150 (MQ=255)
            ccAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCt                                                                                                                                      <  2:1143572/150‑1 (MQ=255)
             cAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCTGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCtt                                                                                                                                     >  1:1705980/1‑149 (MQ=255)
             cAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTg                                                                                                                                     >  1:750686/1‑150 (MQ=255)
              aGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGt                                                                                                                                    <  2:1391574/150‑1 (MQ=255)
              aGCGCCGATCTAGCCCCAGCAGACAAGAACAAGCACATCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGt                                                                                                                                    <  2:2017652/150‑1 (MQ=255)
                cgcCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTAc                                                                                                                                  <  2:1409571/150‑1 (MQ=255)
                 gcCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACa                                                                                                                                 >  2:753450/1‑150 (MQ=255)
                         aGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACAtt                                                                                                                         <  1:930869/150‑1 (MQ=255)
                             ccAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgt                                                                                                                     <  1:269251/150‑1 (MQ=255)
                              cagcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtg                                                                                                                    >  1:465611/1‑150 (MQ=255)
                              cagcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtg                                                                                                                    >  2:1834244/1‑150 (MQ=255)
                              cagcagACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTAGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCGGGCGCTGTACATCCACATTCCgtg                                                                                                                    <  2:27151/150‑1 (MQ=255)
                               agcacACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGCAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGTCGCTGTACATCCACATTCCgtgg                                                                                                                   <  2:1184674/150‑1 (MQ=255)
                                  agACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGc                                                                                                                <  2:2303089/150‑1 (MQ=255)
                                  agACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGc                                                                                                                <  2:211123/150‑1 (MQ=255)
                                     cAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGAAGGAAGGACCCGGCGCGGCAGGCTTTAGCAGCGTGCCCCCCCCGGCCGGGCACAACCACATGCCGTTGTGCGTa                                                                                                             >  1:1615111/1‑150 (MQ=255)
                                         aaCCAGCTCAGTCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGTCGTGCATCCACATTCCGTGGTGCGTACGCa                                                                                                         >  1:1818899/1‑150 (MQ=255)
                                          aCCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCaa                                                                                                        >  2:57452/1‑150 (MQ=255)
                                            cAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtggt                                                                                                                  >  2:1916813/1‑138 (MQ=255)
                                            cAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCgtggt                                                                                                                  <  1:1916813/138‑1 (MQ=255)
                                                   gCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTa                                                                                               >  2:221812/1‑150 (MQ=255)
                                                    ccACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGc                                                                                                                >  2:829702/1‑132 (MQ=255)
                                                    ccACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGc                                                                                                                <  1:829702/132‑1 (MQ=255)
                                                         gcgcCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGa                                                                                         <  2:465458/150‑1 (MQ=255)
                                                          cgcCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGATGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGAc                                                                                        <  1:1435406/150‑1 (MQ=255)
                                                                tGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAAc                                                                                  >  1:1435819/1‑150 (MQ=255)
                                                                tGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAAc                                                                                  >  1:1435826/1‑150 (MQ=255)
                                                                     tGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCa                                                                             <  1:1834244/150‑1 (MQ=255)
                                                                         ccTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCt                                                                                                 <  2:726751/126‑1 (MQ=255)
                                                                         ccTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCt                                                                                                 >  1:726751/1‑126 (MQ=255)
                                                                            gctgcGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGcc                                                                      <  1:932335/150‑1 (MQ=255)
                                                                                 gcAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACt                                                                                 >  2:1828477/1‑134 (MQ=255)
                                                                                 gcAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACt                                                                                 <  1:1828477/134‑1 (MQ=255)
                                                                                 gcAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGcc                                                                 <  1:57452/150‑1 (MQ=255)
                                                                                     ccGTCGGGGCCTGGCATGATCGAAATGCTGCCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGGACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGaa                                                             <  2:465611/149‑1 (MQ=255)
                                                                                      cGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAAc                                                            <  2:665617/150‑1 (MQ=255)
                                                                                          tcggCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGAGACTTCAACTCCCACGCAGCCGGCCCTGAACTAcc                                                        <  2:1277174/148‑1 (MQ=255)
                                                                                          tGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAAc                                                            <  1:2297195/146‑1 (MQ=255)
                                                                                          tGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAAc                                                            >  2:2297195/1‑146 (MQ=255)
                                                                                            ggCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCgg                                                      >  1:2287894/1‑150 (MQ=255)
                                                                                                tGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGa                                                  <  2:2287894/150‑1 (MQ=255)
                                                                                                 ggCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCgg                                                                    >  1:410219/1‑131 (MQ=255)
                                                                                                 ggCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCgg                                                                    <  2:410219/131‑1 (MQ=255)
                                                                                                    aTGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTa                                                                                                             <  2:1318234/87‑1 (MQ=255)
                                                                                                    aTGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTa                                                                                                             >  1:1318234/1‑87 (MQ=255)
                                                                                                     tGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCt                                             >  1:2139579/1‑150 (MQ=255)
                                                                                                       aTCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTAc                                           <  2:2139579/150‑1 (MQ=255)
                                                                                                       aTCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTAc                                           <  1:221812/150‑1 (MQ=255)
                                                                                                       aTCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTAc                                           <  1:753450/150‑1 (MQ=255)
                                                                                                            aaTGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCATGACGACGCTGGCGCTGTACATCCCCATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTcac                                      <  2:1705980/150‑1 (MQ=255)
                                                                                                                cTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCACCTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGcc                                  <  1:582701/150‑1 (MQ=255)
                                                                                                                cTGTCGCCCCCCGGCGCGGCAGGATTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGCGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGcc                                  <  2:100057/150‑1 (MQ=255)
                                                                                                                 tGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGccc                                 >  1:1234882/1‑150 (MQ=255)
                                                                                                                     gCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGcc                                                                                                   >  2:723815/1‑80 (MQ=255)
                                                                                                                     gCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGcc                                                                                                   <  1:723815/80‑1 (MQ=255)
                                                                                                                     gCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCa                                                                         <  1:31524/106‑1 (MQ=255)
                                                                                                                     gCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCa                                                                         >  2:31524/1‑106 (MQ=255)
                                                                                                                     gCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCt                             >  2:2163471/1‑150 (MQ=255)
                                                                                                                      ccccccGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTc                            >  2:1342342/1‑150 (MQ=255)
                                                                                                                                  gCAGGCTTTAGCAGCCTGCCGCCGCTGGCGATGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAACATGCCTACGTCACAGCCCTGCTCGACGACCTCGAc                <  2:1277228/150‑1 (MQ=255)
                                                                                                                                      gCTNTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAgg            >  2:1514130/1‑150 (MQ=255)
                                                                                                                                       cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGa                                                                                         <  1:114752/72‑1 (MQ=255)
                                                                                                                                       cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGa                                                                                         >  2:114752/1‑72 (MQ=255)
                                                                                                                                       cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGa           <  2:2060402/150‑1 (MQ=255)
                                                                                                                                       cTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGa           >  1:2060402/1‑150 (MQ=255)
                                                                                                                                        tttAGCAGCCTGCTGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAg                                    <  1:803479/124‑1 (MQ=255)
                                                                                                                                        tttAGCAGCCTGCTGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAg                                    <  1:803811/124‑1 (MQ=255)
                                                                                                                                        tttAGCAGCCTGCTGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAg                                    >  2:803479/1‑124 (MQ=255)
                                                                                                                                        tttAGCAGCCTGCTGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAg                                    >  2:803811/1‑124 (MQ=255)
                                                                                                                                          tAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgct        <  2:1435819/150‑1 (MQ=255)
                                                                                                                                          tAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgct        <  2:1435826/150‑1 (MQ=255)
                                                                                                                                           agcagcCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAg                                                   >  1:1959189/1‑106 (MQ=255)
                                                                                                                                           agcagcCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAg                                                   <  2:1959189/106‑1 (MQ=255)
                                                                                                                                            gcagcCTGTCGCCGCTGGCTCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTAAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgctgg      <  2:684105/150‑1 (MQ=255)
                                                                                                                                            cgatccTGTCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAgctgg      >  1:684105/5‑150 (MQ=255)
                                                                                                                                                ccTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGCCCACGACCAGGAGCCGGCTGc  >  2:1831433/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
AGCGCAACTGCTCCAGCGCCGATCTAGCCCCAGCAGACAAGAACCAGCTCAGCCACCGCGCCCGTGCCATGGCCCTGCTGCGCAAACGTCTGGGCCTGGCATGATCGAAATGCTGTCGCCCCCCGGCGCGGCAGGCTTTAGCAGCCTGCCGCCGCTGGCGCTGTACATCCACATTCCGTGGTGCGTACGCAAATGCCCTTATTGCGACTTCAACTCCCACGCAGCCGGGCCTGAACTACCGGAAGATGCCTACGTCACAGCCCTGCTCGACGACCTCGACCAGGAGCTGGCTGC  >  NC_002947/5823021‑5823314

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