Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 163,722 G→C 6.8% R38G (CGC→GGC)  scpC ← propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947163,7220GC 6.8% 222.7 / 6.2 73R38G (CGC→GGC) scpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (30/38);  new base C (2/3);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.63e-01

CGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAACTACT  >  NC_002947/163580‑163863
                                                                                                                                              |                                                                                                                                             
cggcggGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGa                                                                                                                                          <  6:61067/148‑1 (MQ=255)
    gggCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGAGTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCg                                                                                                                                     >  8:18880/1‑149 (MQ=255)
      gCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCt                                                                                                                                   <  7:93239/149‑1 (MQ=255)
      gCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaag                                                                                                                                   <  6:32683/149‑4 (MQ=255)
      gCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaag                                                                                                                                   <  4:210521/149‑4 (MQ=255)
       ccAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTc                                                                                                                                  <  7:47568/149‑1 (MQ=255)
        cAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGGTGCGCCGGCGCGGGTGAGGCCGCTCa                                                                                                                                 >  2:122614/1‑149 (MQ=255)
         aaCACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCAt                                                                                                                                <  1:197788/149‑1 (MQ=255)
         aaCACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCAt                                                                                                                                >  3:116880/1‑149 (MQ=255)
           cacccCGGCCTCGGTCAGTTGTTTGTCCAGGGCATTGCCCAGGCTGGCGCCGGTCATCAGGCGGATTTTCAGGGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGc                                                                                                                              <  1:122614/149‑1 (MQ=255)
              acCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGa                                                                                                                           >  7:294204/1‑149 (MQ=255)
              acCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGa                                                                                                                           >  5:79566/1‑149 (MQ=255)
                      cGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATg                                                                                                                   <  8:187560/149‑1 (MQ=255)
                           aGTTGTATGTCGAGGTCTTTGCCCAGGCTGTCGCCGGTCATCAGGCTGATTTTCTGTGGCGTTTGGTTGGGAGGTTCGGCCATTGCATGCGGTAGGGGCTTGGCTTGGCCGGGGCGGGTGAAGCCGCTCATGCCGACTGTCATGCCGTc                                                                                                              <  4:71001/149‑1 (MQ=255)
                             ttgtttgtCCAGGTCATTGCACAGGCTGGCTCCGGTCATCAGGCTGATTTGCAGTGGCGATTGCTTGTCACGTTCGGCCAATGCATGCGTTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACATTCATGCAGTCCt                                                                                                            <  7:155079/149‑1 (MQ=255)
                               gtttgtCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTc                                                                                                           >  5:134080/1‑148 (MQ=255)
                                      cAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGa                                                                                                    <  4:116880/148‑1 (MQ=255)
                                      cAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGAc                                                                                                   >  1:18433/1‑149 (MQ=255)
                                       aGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGAc                                                                                                   >  1:338914/1‑148 (MQ=255)
                                        ggTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGAc                                                                                                   >  5:178652/1‑147 (MQ=255)
                                         gTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAg                                                                                                >  6:242579/1‑149 (MQ=255)
                                            aTTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCg                                                                                              <  6:134080/148‑1 (MQ=255)
                                               gCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaa                                                                                                                                    >  2:105122/1‑105 (MQ=255)
                                               gCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCACCCGCGCCGGCGAAGccaa                                                                                                                                    <  1:105122/107‑3 (MQ=255)
                                                cccAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCg                                                                                              >  6:355854/1‑144 (MQ=255)
                                                cccAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCg                                                                                              <  5:355854/144‑1 (MQ=255)
                                                  cAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGAc                                                                                                   <  3:228253/137‑1 (MQ=255)
                                                  cAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGAc                                                                                                   >  4:228253/1‑137 (MQ=255)
                                                   aGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGc                                                                                      <  5:124227/149‑1 (MQ=255)
                                                    ggCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGcc                                                                                     <  2:338914/149‑1 (MQ=255)
                                                        ggCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCa                                                                                  >  2:321763/1‑148 (MQ=255)
                                                         gcgcCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCa                                                                                  <  8:56569/147‑1 (MQ=255)
                                                             cGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCAtt                                                                            <  1:321763/149‑1 (MQ=255)
                                                               gTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTAc                                                                          <  5:242579/149‑1 (MQ=255)
                                                                           aTTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGc                                                                <  2:207386/147‑1 (MQ=255)
                                                                           aTTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAg                                                              >  2:34225/1‑149 (MQ=255)
                                                                             tttCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAgg                                                             >  1:316981/1‑148 (MQ=255)
                                                                              ttCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAggg                                                            <  7:18880/148‑1 (MQ=255)
                                                                              ttCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAggg                                                            <  6:178652/148‑1 (MQ=255)
                                                                               tCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGgga                                                           <  5:182851/148‑1 (MQ=255)
                                                                                     gCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCAc                                                                                 >  5:118599/1‑120 (MQ=255)
                                                                                     gCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCAc                                                                                 <  6:118599/120‑1 (MQ=255)
                                                                                                gCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGAt                                          >  3:96835/1‑148 (MQ=255)
                                                                                                  aCGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGggag                                                          <  1:41615/130‑1 (MQ=255)
                                                                                                  aCGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGggag                                                          >  2:41615/1‑130 (MQ=255)
                                                                                                   cGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACg                                      <  7:222837/149‑1 (MQ=255)
                                                                                                     ttCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGt                                    <  2:18433/149‑1 (MQ=255)
                                                                                                       cGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCTTGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTa                                   <  6:79566/148‑1 (MQ=255)
                                                                                                        ggCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTAc                                  <  1:34225/148‑1 (MQ=255)
                                                                                                         gCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACAt                                <  3:105579/149‑1 (MQ=255)
                                                                                                         gCCAATGCATGCGGTACGACCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACAt                                <  4:235130/149‑1 (MQ=255)
                                                                                                           cAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGAATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGTACAAGCGGATACGATCACGGTACATg                               >  5:254454/1‑148 (MQ=255)
                                                                                                            aaTGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGCTCACGGTACATGGa                             >  4:272599/1‑149 (MQ=255)
                                                                                                               gcatgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGc                                                   >  3:18306/1‑124 (MQ=255)
                                                                                                               gcatgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGc                                                   <  4:18306/124‑1 (MQ=255)
                                                                                                                catgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAATCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGt                         >  1:275589/1‑149 (MQ=255)
                                                                                                                 atgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCAt                                                                             <  6:209145/96‑1 (MQ=255)
                                                                                                                 atgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCAt                                                                             >  5:209145/1‑96 (MQ=255)
                                                                                                                 atgcGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGt                         >  4:305986/1‑148 (MQ=255)
                                                                                                                      gTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTAtctc                   <  4:96835/149‑1 (MQ=255)
                                                                                                                       tACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCg                  <  3:305986/149‑1 (MQ=255)
                                                                                                                         cGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCgg                 <  5:23839/148‑1 (MQ=255)
                                                                                                                            gagcggCTTCACCCGCGCCGGCGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTa             >  6:68527/5‑149 (MQ=255)
                                                                                                                                 gCTTCGCCTGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACAt                                <  6:251805/125‑1 (MQ=255)
                                                                                                                                 gCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACAt                                >  5:251805/1‑125 (MQ=255)
                                                                                                                                  tttCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGaaa       <  7:261971/148‑1 (MQ=255)
                                                                                                                                   ttCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAac      >  8:150974/1‑149 (MQ=255)
                                                                                                                                   ttCGCCGGCGCGGGAGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAac      >  7:306900/1‑149 (MQ=255)
                                                                                                                                    tCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGaaa       >  8:261971/1‑147 (MQ=255)
                                                                                                                                      gCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAactac   >  5:63508/1‑149 (MQ=255)
                                                                                                                                       ccGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAactact  <  8:294204/149‑1 (MQ=255)
                                                                                                                                        cGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAg                                                              >  3:92139/1‑88 (MQ=255)
                                                                                                                                        cGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAg                                                              <  4:92139/88‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                             
CGGCGGGCCAACACACCGGCCTCGGTCAGTTGTTTGTCCAGGTCATTGCCCAGGCTGGCGCCGGTCATCAGGCTGATTTTCAGTGGCGATTGCTTGGCACGTTCGGCCAATGCATGCGGTACGGCCTTGGCTTCGCCGGCGCGGGTGAAGCCGCTCATGCCGACAGTCATGCCGTCCTCGATCAGACCAGCGGCATCAGCCGCACTCATTACCTTGCTGTGCAGGGAGGACAAGCGGATACGATCACGGTACATGGATTGTTATCTCGGGCTACTGAAACTACT  >  NC_002947/163580‑163863

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