Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,611,142 A→T 6.2% A326A (GCT→GCA 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,1420AT 6.2% 190.4 / 4.6 65A326A (GCT→GCAPP_14123‑hydroxyisobutyryl‑CoA hydrolase
Reads supporting (aligned to +/- strand):  ref base A (29/32);  new base T (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CCCGGACCAGCAGGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCA  >  NC_002947/1610999‑1611290
                                                                                                                                               |                                                                                                                                                    
cccGGACCAGCAGGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAgcgcgc                                                                                                                                                >  1:912046/1‑150 (MQ=255)
     aCCAGCAGGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGcc                                                                                                                                           <  2:1204363/150‑1 (MQ=255)
           agGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACt                                                                                                                                     <  2:608534/150‑1 (MQ=255)
                  cccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCa                                                                                                                              <  2:912046/150‑1 (MQ=255)
                   ccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGTGTGCGCGCTCGCCtgctg                                                                                                                                      >  1:479578/1‑137 (MQ=255)
                   ccACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGTGTGCGCGCTCGCCtgctg                                                                                                                                      <  2:479578/141‑5 (MQ=255)
                           gCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTggcgg                                                                                                                     >  1:200653/1‑150 (MQ=255)
                            ccAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGa                                                                                                                                                      <  2:598239/116‑1 (MQ=255)
                            ccAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGa                                                                                                                                                      >  1:598239/1‑116 (MQ=255)
                                cGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACa                                                                                                                <  2:109650/150‑1 (MQ=255)
                                cGGGTGCCGCCCTTCCCACGACTTGGCAAAATGCGCCTCGACTACCGCCGTCGGAACTTGCGCGCCCGCCTGGCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACa                                                                                                                <  1:526929/150‑1 (MQ=255)
                                     gccgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCa                                                                                                           >  1:1169634/1‑150 (MQ=255)
                                     gccgccCTTCCCACGCCTTCGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAACCCTTTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCa                                                                                                           <  1:612363/150‑1 (MQ=255)
                                       cgccCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAgg                                                                                                         >  1:316697/1‑150 (MQ=255)
                                           cTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCa                                                                                                     >  1:605805/1‑150 (MQ=255)
                                                 aCGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCCGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGt                                                                                               >  2:1233953/1‑150 (MQ=255)
                                                  cGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTa                                                                                              <  1:257052/150‑1 (MQ=255)
                                                  cGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTa                                                                                              <  2:994100/150‑1 (MQ=255)
                                                       tGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCa                                                                                         <  2:1720809/150‑1 (MQ=255)
                                                          cAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCt                                                                                      <  1:1064973/150‑1 (MQ=255)
                                                          cAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCt                                                                                      <  1:1446044/150‑1 (MQ=255)
                                                            aaaTGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTgg                                                                                    >  2:171243/1‑150 (MQ=255)
                                                                          tACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACt                                                                                                                                     >  2:1521204/1‑87 (MQ=255)
                                                                          tACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACt                                                                                                                                     <  1:1521204/87‑1 (MQ=255)
                                                                          tACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCaa                                                                      >  1:1259120/1‑150 (MQ=255)
                                                                           aCCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAAc                                                                     <  1:984741/150‑1 (MQ=255)
                                                                            ccgccgTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCAAACg                                                                    >  1:972434/1‑150 (MQ=255)
                                                                                 gTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTAGTCTCTGTCCAGCAGGCGAGCGCGCACGCCTCGGCTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAgg                                                               <  1:435884/150‑1 (MQ=255)
                                                                                    ggCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTg                                                                                                <  1:1017641/114‑1 (MQ=255)
                                                                                    ggCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTg                                                                                                >  2:1017641/1‑114 (MQ=255)
                                                                                    ggCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTAc                                                            >  1:451930/1‑150 (MQ=255)
                                                                                    ggCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTAc                                                            >  1:1940360/1‑150 (MQ=255)
                                                                                     gCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTAcc                                                           >  2:1266841/1‑150 (MQ=255)
                                                                                        cTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACcgcg                                                        <  2:1259120/150‑1 (MQ=255)
                                                                                         ttGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCgcgc                                                       >  1:1238830/1‑150 (MQ=255)
                                                                                             gcCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCCACGCCAGGTACCGCGCCcgc                                                   >  2:798058/1‑149 (MQ=255)
                                                                                              ccACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcggc                                                  >  1:6294/1‑150 (MQ=255)
                                                                                               cacaTCAGGCCAGTACCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcggcg                                                 >  1:984022/1‑150 (MQ=255)
                                                                                                 caTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcgg                                                   >  2:1361550/1‑146 (MQ=255)
                                                                                                 caTCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcgg                                                   <  1:1361550/146‑1 (MQ=255)
                                                                                                   tCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATc                                             <  1:773539/150‑1 (MQ=255)
                                                                                                         cAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTcc                                       <  1:798058/150‑1 (MQ=255)
                                                                                                          aGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTccc                                      <  2:451930/150‑1 (MQ=255)
                                                                                                          aGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTccc                                      <  1:171243/150‑1 (MQ=255)
                                                                                                          aGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTccc                                      <  2:1594711/150‑1 (MQ=255)
                                                                                                           gTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCcca                                     <  2:200653/150‑1 (MQ=255)
                                                                                                                  cgcgGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCccagaccaga                              >  1:768557/1‑150 (MQ=255)
                                                                                                                        tGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGc                                          <  2:1297511/132‑1 (MQ=255)
                                                                                                                        tGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGc                                          >  1:1297511/1‑132 (MQ=255)
                                                                                                                        tGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGcc                        <  2:984022/150‑1 (MQ=255)
                                                                                                                           ttGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCACCAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGc                                          <  1:1106420/129‑1 (MQ=255)
                                                                                                                           ttGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGc                                          >  2:1106420/1‑129 (MQ=255)
                                                                                                                            tGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCa                    <  1:1906267/150‑1 (MQ=255)
                                                                                                                              tCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGCCCAGATGCGCCGTCAGc                  >  1:274636/1‑150 (MQ=255)
                                                                                                                               cTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGaa                                                                                  <  2:562255/85‑1 (MQ=255)
                                                                                                                               cTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGaa                                                                                  >  1:562255/1‑85 (MQ=255)
                                                                                                                               cTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCg                 >  1:1515245/1‑150 (MQ=255)
                                                                                                                               cTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCg                 <  2:768557/150‑1 (MQ=255)
                                                                                                                                      cagcagGCGAGCGCGCACGCCTTAGCGGAATGCTGGGGTGCTGCTACAGTTGAGGCTCATGCTGTATTCCATCTGGAACACATGCGCCAACGAAAATCAGAACACCGTCCTGATGGGCTCCCAGACCAGATGCGCCGTCAGCGGGCAAcc          >  2:1876604/1‑150 (MQ=255)
                                                                                                                                       tgaagGCGTGCGCGCTCGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcg                                                    <  1:457932/104‑1 (MQ=255)
                                                                                                                                       tgaagGCGTGCGCGCTCGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcg                                                    >  2:457932/4‑107 (MQ=255)
                                                                                                                                          agGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcggc                                                  <  2:807541/106‑1 (MQ=255)
                                                                                                                                          agGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCcggc                                                  >  1:807541/1‑106 (MQ=255)
                                                                                                                                          agGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCa      >  1:1210592/1‑150 (MQ=255)
                                                                                                                                          agGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCa      <  2:316697/150‑1 (MQ=255)
                                                                                                                                           ggCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCAt     <  2:6294/150‑1 (MQ=255)
                                                                                                                                              gAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCa  >  1:476645/1‑150 (MQ=255)
                                                                                                                                               aGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGc              >  2:370909/1‑137 (MQ=255)
                                                                                                                                               aGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGc              <  1:370909/137‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CCCGGACCAGCAGGCTAACCCACATCCGCCAACGGGTGCCGCCCTTCCCACGCCTTGGCAAAATGCGCCTCGACTACCGCCGTCGGCACTTGCGCCACATCAGGCCAGTGCCAGCGCGGCTGGTTGTCTTTGTCCAGCAGGCGAGCGCGCACGCCTTCACTGAATTCAGGGTGGCGGCAACAGTTCAGGCTCATGCTGTATTCCATCTGGAACACCTGCGCCAACGACAGGTACCGCGCCCGGCGGATCTGCTCCCAGACCAGATGCGCCGTCAGCGGGCAACCCTCATGCA  >  NC_002947/1610999‑1611290

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: