Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,611,149 A→T 5.6% V324E (GTG→GAG)  PP_1412 ← 3‑hydroxyisobutyryl‑CoA hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,611,1490AT 5.6% 315.0 / 6.2 103V324E (GTG→GAG) PP_14123‑hydroxyisobutyryl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base A (48/48);  new base T (3/3);  total (52/51)
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.99e-01

Only 100 of 104 total aligned reads displayed.

CCGGACCAGCAGGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGA  >  NC_002947/1611000‑1611292
                                                                                                                                                     |                                                                                                                                               
ccGGACCAGCAGGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGAATTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCa                                                                                                                                                 >  1:1996716/1‑150 (MQ=255)
        gcagGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCa                                                                                                                                         <  2:261272/150‑1 (MQ=255)
                 cccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCa                                                                                                                                >  2:2282707/1‑150 (MQ=255)
                  ccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGTGCGCGCTCGCCtgct                                                                                                                                         >  1:1179963/1‑137 (MQ=255)
                  ccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGTGCGCGCTCGCCtgct                                                                                                                                         <  2:1179963/140‑4 (MQ=255)
                  ccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCa                                                                                                                                <  2:811740/149‑1 (MQ=255)
                  ccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCa                                                                                                                                >  1:811740/1‑149 (MQ=255)
                  ccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAg                                                                                                                               <  1:2302432/150‑1 (MQ=255)
                  ccACATCCGCCAACGGGTGCCGCCCTACCGACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAg                                                                                                                               >  1:1025028/1‑150 (MQ=255)
                  ccACATCAGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTGGACTACCGCCGTCGGCACTTGCGCCCCATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAg                                                                                                                               >  2:2365350/1‑150 (MQ=255)
                      aTCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTg                                                                                                                           >  1:2273991/1‑150 (MQ=255)
                      aTCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTg                                                                                                                           >  1:1867845/1‑150 (MQ=255)
                          gCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTggcgg                                                                                                                       >  1:2171134/1‑150 (MQ=255)
                          gCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTggcgg                                                                                                                       >  1:72010/1‑150 (MQ=255)
                          gCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTggcgg                                                                                                                       >  1:2437095/1‑150 (MQ=255)
                          gCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTggcgg                                                                                                                       >  2:1903476/1‑150 (MQ=255)
                           ccAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTggcggc                                                                                                                      >  2:235190/1‑150 (MQ=255)
                                    gccgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCa                                                                                                             <  1:1965046/150‑1 (MQ=255)
                                    gccgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACCTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCa                                                                                                             <  1:980008/150‑1 (MQ=255)
                                      cgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGTAGGTGAGCGCGCACGCTTTCGCTTGATTCTTTGTGGCGGCAACAGTTCAgg                                                                                                           >  2:1946112/1‑150 (MQ=255)
                                      cgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGTGCGCGCTCGCCtgctg                                                                                                                                        >  1:1516925/1‑117 (MQ=255)
                                      cgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGTGCGCGCTCGCCtgctg                                                                                                                                        <  2:1516925/121‑5 (MQ=255)
                                      cgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAgg                                                                                                           <  2:2273991/150‑1 (MQ=255)
                                       gccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGc                                                                                                          <  2:122698/150‑1 (MQ=255)
                                          cTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCa                                                                                                       <  2:1996716/150‑1 (MQ=255)
                                           ttCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCAt                                                                                                      >  1:435193/1‑150 (MQ=255)
                                             cccACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGc                                                                                                    >  2:1249117/1‑150 (MQ=255)
                                              ccACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAg                                                                                                                               >  1:791349/1‑122 (MQ=255)
                                              ccACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAg                                                                                                                               <  2:791349/122‑1 (MQ=255)
                                              ccACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCt                                                                                                   >  2:494327/1‑150 (MQ=255)
                                              ccACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCt                                                                                                   <  2:1867845/150‑1 (MQ=255)
                                              ccACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTAGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCt                                                                                                   >  1:69146/1‑150 (MQ=255)
                                               cACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTg                                                                                                  >  1:1512410/1‑150 (MQ=255)
                                                 cGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTa                                                                                                <  2:1630329/150‑1 (MQ=255)
                                                 cGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTa                                                                                                <  2:1630241/150‑1 (MQ=255)
                                                  gCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTAt                                                                                               <  1:1491274/150‑1 (MQ=255)
                                                    cTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTc                                                                                             >  2:1277200/1‑150 (MQ=255)
                                                         cAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCt                                                                                        <  2:1558918/150‑1 (MQ=255)
                                                         cAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCt                                                                                        <  2:2297991/150‑1 (MQ=255)
                                                         cAAAACGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCt                                                                                        <  1:245194/150‑1 (MQ=255)
                                                          aaaaTGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTg                                                                                       >  1:1222492/1‑150 (MQ=255)
                                                                   cTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTg                                                                              >  2:567280/1‑150 (MQ=255)
                                                                   cTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGCCATTGTCCAGCCGATAACCGCGCGCGTGTGCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTg                                                                              >  1:1484596/1‑150 (MQ=255)
                                                                    tCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgc                                                                             <  1:2333040/150‑1 (MQ=255)
                                                                    tCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgc                                                                             >  1:2392698/1‑150 (MQ=255)
                                                                    tCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgc                                                                             <  1:539844/150‑1 (MQ=255)
                                                                    tCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTgc                                                                             <  1:498540/150‑1 (MQ=255)
                                                                         tACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCaa                                                                                                                    >  2:754768/1‑106 (MQ=255)
                                                                         tACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAgg                                                                                                           >  1:135610/1‑115 (MQ=255)
                                                                         tACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAgg                                                                                                           <  2:135610/115‑1 (MQ=255)
                                                                         tACCGCCGTCGGCACTTGCGCCACACCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCaa                                                                                                                    <  1:754768/106‑1 (MQ=255)
                                                                          aCCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAAc                                                                       <  1:567280/150‑1 (MQ=255)
                                                                             gccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGAc                                                                    <  1:1249117/150‑1 (MQ=255)
                                                                                    gCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTAcc                                                             >  1:1093532/1‑150 (MQ=255)
                                                                                      aCTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACcgc                                                           <  1:2433551/150‑1 (MQ=255)
                                                                                        ttGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCgcgc                                                         >  1:1727479/1‑150 (MQ=255)
                                                                                            gcCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCt                                                                                        >  2:1672662/1‑115 (MQ=255)
                                                                                            gcCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCt                                                                                        <  1:1672662/115‑1 (MQ=255)
                                                                                            gcCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcgg                                                     >  1:2221414/1‑150 (MQ=255)
                                                                                             ccACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcggc                                                    >  2:1213671/1‑150 (MQ=255)
                                                                                              cacaTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACa                                                                                                                  >  1:416394/1‑87 (MQ=255)
                                                                                              cacaTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACa                                                                                                                  <  2:416394/87‑1 (MQ=255)
                                                                                                 aTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGAt                                                >  1:1625721/1‑150 (MQ=255)
                                                                                                 aTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGAt                                                >  1:2284370/1‑150 (MQ=255)
                                                                                                  tCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACa                                                                                                                  >  2:2142992/1‑83 (MQ=255)
                                                                                                  tCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACa                                                                                                                  <  1:2142992/83‑1 (MQ=255)
                                                                                                       cgactGCCAGCGCGGCTGGGTGGCGTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTc                                          <  2:1093532/146‑1 (MQ=255)
                                                                                                       ccAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTgg                                                                                                                          <  1:567069/70‑1 (MQ=255)
                                                                                                       ccAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTgg                                                                                                                          >  2:567069/1‑70 (MQ=255)
                                                                                                        cAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTcc                                         >  1:1712438/1‑150 (MQ=255)
                                                                                                        cAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTcc                                         <  2:1625721/150‑1 (MQ=255)
                                                                                                             ccAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCTCGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagac                                    <  2:1484596/150‑1 (MQ=255)
                                                                                                             ccAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagac                                    <  2:969215/150‑1 (MQ=255)
                                                                                                               aGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagacca                                  <  2:2284370/150‑1 (MQ=255)
                                                                                                                gcgcGGCTGGTTGTCTTTGTACAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagaccag                                 >  1:1253502/1‑150 (MQ=255)
                                                                                                                  gcgGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGAt                               <  2:1222492/150‑1 (MQ=255)
                                                                                                                  gcgGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGAt                               <  2:1856114/150‑1 (MQ=255)
                                                                                                                    ggCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGNTgc                             <  2:1512410/150‑1 (MQ=255)
                                                                                                                     gCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAg                                                                                                            <  1:2047268/70‑1 (MQ=255)
                                                                                                                     gCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAg                                                                                                            >  2:2047268/1‑70 (MQ=255)
                                                                                                                       tGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGcc                          <  2:2221414/150‑1 (MQ=255)
                                                                                                                          ttGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTc                                          >  1:1372714/1‑131 (MQ=255)
                                                                                                                          ttGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTc                                          <  2:1372714/131‑1 (MQ=255)
                                                                                                                          ttGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTc                       >  1:225324/1‑150 (MQ=255)
                                                                                                                              cTTTGTCCAGCAGGCGAGCGCGCACGCCTTGACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCg                   >  2:1407488/1‑150 (MQ=255)
                                                                                                                              cTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCg                   <  1:1213671/150‑1 (MQ=255)
                                                                                                                               tttGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATgcgc                           <  2:2140891/141‑1 (MQ=255)
                                                                                                                               tttGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATgcgc                           >  1:2140891/1‑141 (MQ=255)
                                                                                                                                  gTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagacca                                  <  1:1764287/131‑1 (MQ=255)
                                                                                                                                  gTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagacca                                  >  2:1764287/1‑131 (MQ=255)
                                                                                                                                   tCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGcc                          <  1:2371331/138‑1 (MQ=255)
                                                                                                                                   tCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGcc                          >  2:2371331/1‑138 (MQ=255)
                                                                                                                                     cagcagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAAcc            >  2:1704549/1‑150 (MQ=255)
                                                                                                                                       gcagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCaa                                                                                                                    >  2:1259424/1‑44 (MQ=255)
                                                                                                                                       gcagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCaa                                                                                                                    <  1:1259424/44‑1 (MQ=255)
                                                                                                                                        cagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTc         <  1:1903476/150‑1 (MQ=255)
                                                                                                                                        cagGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTc         <  2:72010/150‑1 (MQ=255)
                                                                                                                                        aagGCGTGCGCGCTCGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACTTGTGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTc         >  1:1972901/2‑150 (MQ=255)
                                                                                                                                               gcgcgcACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGa  <  2:2171134/150‑1 (MQ=255)
                                                                                                                                                   gcACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGAc                                                                    <  2:90970/80‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CCGGACCAGCAGGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCAGA  >  NC_002947/1611000‑1611292

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: