Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,715,401 A→T 9.1% D166E (GAT→GAA)  paaJ ← 3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,715,4010AT 9.1% 164.4 / 7.5 55D166E (GAT→GAA) paaJ3‑oxoadipyl‑CoA/3‑oxo‑5,6‑dehydrosuberyl‑CoA thiolase
Reads supporting (aligned to +/- strand):  ref base A (22/28);  new base T (3/2);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

GATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCAT  >  NC_002947/3715254‑3715550
                                                                                                                                                   |                                                                                                                                                     
gATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaatc                                                                                                                                                     >  2:805973/1‑150 (MQ=255)
  tttCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaa                                                                                                                                                       <  2:1199105/146‑1 (MQ=255)
  tttCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCatcgaa                                                                                                                                                       >  1:1199105/1‑146 (MQ=255)
      gACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGcc                                                                                                                                               >  2:821661/1‑150 (MQ=255)
       aCCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGATTCGATGCCg                                                                                                                                              <  2:619578/150‑1 (MQ=255)
       aCCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCg                                                                                                                                              <  1:805973/150‑1 (MQ=255)
              cGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAggcgg                                                                                                                                       <  2:1391294/150‑1 (MQ=255)
                aTTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGcc                                                                                                                                     <  2:29986/150‑1 (MQ=255)
                 tttCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCt                                                                                                                                    >  1:596758/1‑150 (MQ=255)
                       gcgcCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTtcatca                                                                                                                              <  1:768705/150‑1 (MQ=255)
                         gcCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGc                                                                                                                            <  1:1376904/150‑1 (MQ=255)
                           cAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTtcattttgag                                                                                                                          >  1:288906/1‑144 (MQ=255)
                           cAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCgg                                                                                                                          <  1:266584/150‑1 (MQ=255)
                              gcggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTTGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGtt                                                                                                                       <  1:216904/150‑1 (MQ=255)
                              gcggcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGtt                                                                                                                       <  2:655575/150‑1 (MQ=255)
                                 gcCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGAc                                                                                                                    <  2:1237477/150‑1 (MQ=255)
                                         ccTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGcc                                                                                                            >  1:568086/1‑150 (MQ=255)
                                          cTGCGCTGCTGCGGCTTTGTGCTGGCTGCGCATGGCAAAGGCATCCGGGTCGGCTCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCa                                                                                                           <  2:470271/150‑1 (MQ=255)
                                               ctgctgCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGa                                                                                                      <  1:821661/150‑1 (MQ=255)
                                                tgctgcGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGAt                                                                                                     <  1:1230604/150‑1 (MQ=255)
                                                        cTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGgtgt                                                                                             >  1:140878/1‑150 (MQ=255)
                                                                     gcgcAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTGTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGc                                                                                <  2:754871/150‑1 (MQ=255)
                                                                       gcAGGGCAAAGGCCTCCTGGTCGGCGAGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGAATCGAATCGATGCCGTAGGCGGCCTTCATCACCGGGTTGACGAAACGCAAGCCGATGGTGGTGTCGAACCGCTcctccg                                                                              >  2:585468/1‑145 (MQ=255)
                                                                       gcAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCg                                                                              >  1:588610/1‑150 (MQ=255)
                                                                           ggCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCa                                                                                                           >  2:841974/1‑117 (MQ=255)
                                                                           ggCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCa                                                                                                           <  1:841974/117‑1 (MQ=255)
                                                                             cAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGcc                                                                        >  2:884515/1‑150 (MQ=255)
                                                                                             ggCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGa                                                        >  1:1213950/1‑150 (MQ=255)
                                                                                                   cgAGATGCCGAACTGTTCGTCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                                <  1:1081094/88‑1 (MQ=255)
                                                                                                   cgAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGaaa                                                                                                                >  2:1081094/1‑88 (MQ=255)
                                                                                                     agaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCa                                                <  2:1406745/150‑1 (MQ=255)
                                                                                                      gaTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCCTCAGCCGGCCTATGAAACGCCACGCCATGGTTGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCAt                                               <  1:884515/150‑1 (MQ=255)
                                                                                                       tggCCGAACAGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATc                                              >  2:1352446/3‑150 (MQ=255)
                                                                                                       aTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCgg                                                                                                                          <  2:154049/74‑1 (MQ=255)
                                                                                                       aTGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCgg                                                                                                                          >  1:154049/1‑74 (MQ=255)
                                                                                                         gCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATcac                                            >  2:3964/1‑150 (MQ=255)
                                                                                                          ccGAACTGTTCGGCCACGTTTTCATCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATcaca                                           >  1:404558/1‑150 (MQ=255)
                                                                                                                tGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCGTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCAt                                               <  2:871468/140‑1 (MQ=255)
                                                                                                                tGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCAt                                               >  1:871468/1‑140 (MQ=255)
                                                                                                                    cGGCCACGTCTGCAGCCGGTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAc                                 <  1:1092125/150‑1 (MQ=255)
                                                                                                                        cACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGa                             <  2:1213950/150‑1 (MQ=255)
                                                                                                                         aCGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGAc                            >  1:768610/1‑150 (MQ=255)
                                                                                                                           gTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACAt                          <  1:999579/150‑1 (MQ=255)
                                                                                                                            ttttCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATg                         <  2:596758/150‑1 (MQ=255)
                                                                                                                              ttCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGa                       <  2:803820/150‑1 (MQ=255)
                                                                                                                                 aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATggtggt                                                                                               >  2:1335156/1‑75 (MQ=255)
                                                                                                                                 aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATggtggt                                                                                               <  1:1335156/75‑1 (MQ=255)
                                                                                                                                 aGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTc                    >  1:624446/1‑150 (MQ=255)
                                                                                                                                 aGCCGTTTCCGGCATAGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTc                    >  2:1378240/1‑150 (MQ=255)
                                                                                                                                   ccGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCca                  >  2:641078/1‑150 (MQ=255)
                                                                                                                                    cGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCcac                 <  2:624446/150‑1 (MQ=255)
                                                                                                                                     gTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCcaca                <  1:641078/150‑1 (MQ=255)
                                                                                                                                       ctacGGCATCGATTCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAc                                                   >  2:1182389/4‑113 (MQ=255)
                                                                                                                                       ctacGGCATCGATTCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTAc                                                   <  1:1182389/110‑1 (MQ=255)
                                                                                                                                       ctacGGCATCGATTCGATGCCGGAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCAcacc              >  1:63569/4‑150 (MQ=255)
                                                                                                                                           ggCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCg          <  1:520497/150‑1 (MQ=255)
                                                                                                                                           ggCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCg          >  1:1204455/1‑150 (MQ=255)
                                                                                                                                             cATCGATTCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGc        >  2:500543/1‑150 (MQ=255)
                                                                                                                                               tcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgc                               >  1:967642/1‑125 (MQ=255)
                                                                                                                                               tcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCAcgc                               <  2:967642/125‑1 (MQ=255)
                                                                                                                                               tcgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGGTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCa      >  1:756671/1‑150 (MQ=255)
                                                                                                                                                cgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACa                                                                                       >  2:1204926/1‑68 (MQ=255)
                                                                                                                                                cgaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACa                                                                                       <  1:1204926/68‑1 (MQ=255)
                                                                                                                                                 gaatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGc    >  1:701429/1‑150 (MQ=255)
                                                                                                                                                  aatcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCa   >  2:577995/1‑150 (MQ=255)
                                                                                                                                                   atcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCAt  <  2:404558/150‑1 (MQ=255)
                                                                                                                                                   atcgaTGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCAt  >  2:320850/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATTTCGACCGGCACGATTTCCCGCGCCAGGCGGCCGCGGGCCTGAGCTGCTGCGGCTTTGTGCTGGCTGCGCAGGGCAAAGGCATCCTGGTCGGCGCGCGAGATGCCGAACTGTTCGGCCACGTTTTCAGCCGTTTCCGGCATCGAATCGATGCCGTAGGCGGCCTTCATCAGCGGGTTGACGAAACGCCAGCCGATGGTGGTGTCGAACAGCTCGGCCGCACGCCCGAATGCCTGCTCCGACTTACCCATCACAAACGGTGCACGCGACATGGACTCCACACCACCGGCCAGCAT  >  NC_002947/3715254‑3715550

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