Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,874,876 A→C 6.3% S451A (TCG→GCG) ‡ ilvA‑II ← threonine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,874,8760AC 6.3% 318.8 / 11.6 111S451A (TCG→GCG) ‡ilvA‑IIthreonine deaminase
Reads supporting (aligned to +/- strand):  ref base A (50/54);  new base C (3/4);  total (53/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.52e-01

Only 100 of 117 total aligned reads displayed.

TTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTG  >  NC_002947/5874727‑5875024
                                                                                                                                                     |                                                                                                                                                    
tttGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCga                                                                                                                                                      <  1:2393102/150‑1 (MQ=255)
       gCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                               >  2:1225735/1‑150 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCCCGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGAt                                                                                                                                                   >  2:148657/1‑144 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGTGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCa                                                                                                                                             >  1:473922/1‑150 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTCCa                                                                                                                                             <  2:1178093/150‑1 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTCCa                                                                                                                                             <  2:2111944/150‑1 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGAt                                                                                                                                                   <  1:148657/144‑1 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCa                                                                                                                                             <  1:1219280/150‑1 (MQ=255)
         cGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCa                                                                                                                                             <  1:1225735/150‑1 (MQ=255)
           ggTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAgc                                                                                                                                           >  2:1033214/1‑150 (MQ=255)
             ttttCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAgcgc                                                                                                                                         >  2:669727/1‑150 (MQ=255)
              tttCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGcgcc                                                                                                                                        <  1:669727/150‑1 (MQ=255)
                        gTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCACCCGCCCAGCtt                                                                                                                              <  2:2172746/150‑1 (MQ=255)
                         tCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttg                                                                                                                             >  2:1863588/1‑150 (MQ=255)
                          cccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgt                                                                                                                            >  1:183134/1‑150 (MQ=255)
                          cccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgt                                                                                                                            >  2:456661/1‑150 (MQ=255)
                           ccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgtt                                                                                                                           <  2:1980209/150‑1 (MQ=255)
                           ccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttgtt                                                                                                                           >  1:732693/1‑150 (MQ=255)
                               tACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                               >  1:1185631/1‑126 (MQ=255)
                                 cGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGaa                                                                                                                     <  2:735139/150‑1 (MQ=255)
                                   ggTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGt                                                                                                                   >  2:1773889/1‑150 (MQ=255)
                                   ggTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGt                                                                                                                   <  1:486599/150‑1 (MQ=255)
                                    gTAACCGATCTTGGCCAGCGCTGCCGGCACAAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGtt                                                                                                                  >  1:961792/1‑150 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGAt                                                                                                                                                   <  2:1275984/112‑1 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGAt                                                                                                                                                   >  1:1275984/1‑112 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttg                                                                                                                             <  2:1458494/134‑1 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCttg                                                                                                                             >  1:1458494/1‑134 (MQ=255)
                                         cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACa                                                                                                             <  2:732693/150‑1 (MQ=255)
                                            tCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTg                                                                                                          >  1:704490/1‑150 (MQ=255)
                                                ggCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGACCGCCATGCTTGTTGAGGAAGTTGAACAGTGcccc                                                                                                      <  1:1863588/150‑1 (MQ=255)
                                                  ccAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTg                                                                                                    >  1:1428528/1‑150 (MQ=255)
                                                  ccAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTg                                                                                                    >  1:1012407/1‑150 (MQ=255)
                                                    aGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGTGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGc                                                                                                  <  2:473922/150‑1 (MQ=255)
                                                       gcTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGc                                                                                               >  2:198086/1‑150 (MQ=255)
                                                         tGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGc                                                                                                                                >  1:2073/1‑115 (MQ=255)
                                                         tGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTc                                                                                             >  2:124620/1‑150 (MQ=255)
                                                         tGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTc                                                                                             >  2:2060353/1‑150 (MQ=255)
                                                           ccGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCTGGTTCCCTCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCgg                                                                                           >  2:64505/1‑150 (MQ=255)
                                                           ccGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCgg                                                                                           <  1:198086/150‑1 (MQ=255)
                                                           ccGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCgg                                                                                           <  1:1773889/150‑1 (MQ=255)
                                                            cGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCggg                                                                                          <  1:1629102/150‑1 (MQ=255)
                                                               cACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGaa                                                                                       <  1:2386089/150‑1 (MQ=255)
                                                                aCCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAt                                                                                      <  1:2319405/150‑1 (MQ=255)
                                                                aCCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAt                                                                                      <  1:1033214/150‑1 (MQ=255)
                                                                aCCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAt                                                                                      <  2:436253/150‑1 (MQ=255)
                                                                 ccAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAtt                                                                                     >  2:2186566/1‑150 (MQ=255)
                                                                  cAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTCCa                                                                                                                                             <  2:1411574/93‑1 (MQ=255)
                                                                     gTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                               <  1:400490/88‑1 (MQ=255)
                                                                     gTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                               >  2:400490/1‑88 (MQ=255)
                                                                          gTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAgcgc                                                                             <  1:136281/149‑1 (MQ=255)
                                                                              gTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGaa                                                                                                                     >  2:488225/1‑105 (MQ=255)
                                                                              gTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGaa                                                                                                                     <  1:488225/105‑1 (MQ=255)
                                                                               tcctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCacta                                                                                                                                          >  2:2256823/1‑80 (MQ=255)
                                                                               tcctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCacta                                                                                                                                          <  1:2256823/83‑4 (MQ=255)
                                                                                 ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAt                                                                     >  2:462926/1‑150 (MQ=255)
                                                                                  tccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAtt                                                                    >  2:1105823/1‑150 (MQ=255)
                                                                                   ccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCCGCACCAttt                                                                   <  2:2168961/150‑1 (MQ=255)
                                                                                     ggCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGa                                                                                                                         >  1:967941/1‑94 (MQ=255)
                                                                                      gCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGaa                                                                                                                     >  2:568802/1‑97 (MQ=255)
                                                                                      gCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGaa                                                                                                                     <  1:568802/97‑1 (MQ=255)
                                                                                        aCCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcgtcg                                                              >  2:1462164/1‑150 (MQ=255)
                                                                                        aCCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTAAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcgtcg                                                              >  1:148502/1‑150 (MQ=255)
                                                                                         ccTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATTGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGc                                                             <  2:191753/150‑1 (MQ=255)
                                                                                          cTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCt                                                            <  2:2075140/150‑1 (MQ=255)
                                                                                          cTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCt                                                            >  1:375727/1‑150 (MQ=255)
                                                                                               gCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCAc                                                       <  2:1428528/150‑1 (MQ=255)
                                                                                                  ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCgg                                                    >  1:1106144/1‑150 (MQ=255)
                                                                                                      caccacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                                <  2:1595636/150‑1 (MQ=255)
                                                                                                      caccacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGGTGAGGAAGTTGATCAGTGCCCCGGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                                <  2:331045/150‑1 (MQ=255)
                                                                                                         cacGCGCCCGTCAGCCGCGGCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCggacctgat                                             <  2:1854887/150‑8 (MQ=255)
                                                                                                         cacGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGCGTg                                             >  2:745504/1‑150 (MQ=255)
                                                                                                          acGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTgg                                            >  2:2155586/1‑150 (MQ=255)
                                                                                                                ccGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGccgccg                                      >  2:1113924/1‑150 (MQ=255)
                                                                                                                ccGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTCCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                           >  1:2007856/1‑129 (MQ=255)
                                                                                                                ccGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTCCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                           <  2:2007856/129‑1 (MQ=255)
                                                                                                                        cgcgCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATg                              >  2:1267459/1‑150 (MQ=255)
                                                                                                                          cgcCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGAc                            <  2:375727/150‑1 (MQ=255)
                                                                                                                            ccATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGa                          <  2:183134/150‑1 (MQ=255)
                                                                                                                             cATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAGTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGaa                         <  2:1012407/150‑1 (MQ=255)
                                                                                                                              aTGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCACCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAat                        >  1:1774296/1‑150 (MQ=255)
                                                                                                                                ggTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGTCGCCGACCATATGACGAatat                      <  1:124620/150‑1 (MQ=255)
                                                                                                                                   tGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCa                   <  1:1113924/150‑1 (MQ=255)
                                                                                                                                    gCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAg                  >  1:1895799/1‑150 (MQ=255)
                                                                                                                                        tAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTg              <  1:2186566/150‑1 (MQ=255)
                                                                                                                                        tAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTg              <  1:2155586/150‑1 (MQ=255)
                                                                                                                                        tAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTg              <  2:1223953/150‑1 (MQ=255)
                                                                                                                                        tAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTg              <  1:1909180/150‑1 (MQ=255)
                                                                                                                                           tGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAGCAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAata                       >  1:2291055/1‑138 (MQ=255)
                                                                                                                                           tGGAACATCTCGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAGCAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAata                       <  2:2291055/138‑1 (MQ=255)
                                                                                                                                           tGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGaa                                                                                                               <  2:1693722/50‑1 (MQ=255)
                                                                                                                                           tGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGcc           >  1:727974/1‑150 (MQ=255)
                                                                                                                                            ggAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGaa                                                                                       <  1:141979/73‑1 (MQ=255)
                                                                                                                                            ggAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGaa                                                                                       >  2:141979/1‑73 (MQ=255)
                                                                                                                                              aaCATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAAt        <  1:2216456/150‑1 (MQ=255)
                                                                                                                                                  tCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCata                                <  2:1952944/122‑1 (MQ=255)
                                                                                                                                                    gagaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAg                                                                                <  2:1263464/72‑1 (MQ=255)
                                                                                                                                                    gagaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAg                                                                                >  1:1263464/1‑72 (MQ=255)
                                                                                                                                                    gagaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTg  <  1:1267459/150‑1 (MQ=255)
                                                                                                                                                     agaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGt                                                                                                                   <  2:25111/36‑1 (MQ=255)
                                                                                                                                                     agaTGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGt                                                                                                                   >  1:25111/1‑36 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTTGTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCAGCGCCCGCCCAGCTTGTTGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTG  >  NC_002947/5874727‑5875024

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: