Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 163,973 T→A 5.0% intergenic (‑140/+47) scpC ← / ← pntB propionyl‑CoA:succinate CoA transferase/pyridine nucleotide transhydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947163,9730TA 5.0% 304.3 / 4.8 99intergenic (‑140/+47)scpC/pntBpropionyl‑CoA:succinate CoA transferase/pyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (51/43);  new base A (2/3);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.31e-01

TGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATA  >  NC_002947/163833‑164121
                                                                                                                                             |                                                                                                                                                    
tGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCg                                                                                                                                              <  4:8137/149‑1 (MQ=255)
tGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCg                                                                                                                                              >  7:188640/1‑149 (MQ=255)
  gATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTgccgg                                                                                                                                              <  7:332002/148‑5 (MQ=255)
  gATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTGCCggg                                                                                                                                            <  7:139287/149‑1 (MQ=255)
  gATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCggg                                                                                                                                            <  7:195215/149‑1 (MQ=255)
  gATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCgg                                                                                                                                             <  6:54475/148‑1 (MQ=255)
  gATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTAGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCggg                                                                                                                                            <  1:44034/149‑1 (MQ=255)
       ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGGTTCTGACCGGGGCTTc                                                                                                                                       <  2:118588/149‑1 (MQ=255)
       ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTc                                                                                                                                       >  7:63292/1‑149 (MQ=255)
        tATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTcc                                                                                                                                      >  3:108590/1‑149 (MQ=255)
              gggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTGCCGGGGCttttctt                                                                                                                                   <  3:277467/146‑6 (MQ=255)
              gggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTtat                                                                                                                                 <  6:91642/148‑1 (MQ=255)
              gggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTtat                                                                                                                                 <  3:383001/148‑1 (MQ=255)
               ggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTatat                                                                                                                               <  8:63292/149‑1 (MQ=255)
                gCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTatat                                                                                                                               >  2:350093/1‑148 (MQ=255)
                  tACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCtt                                                                                                                                    >  7:72170/1‑141 (MQ=255)
                  tACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCtt                                                                                                                                    <  8:72170/141‑1 (MQ=255)
                  tACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTatat                                                                                                                               >  5:356241/1‑146 (MQ=255)
                   aCTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGc                                                                                                                           >  2:223691/1‑149 (MQ=255)
                    cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCATTAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCt                                                                                                                          >  4:26765/1‑149 (MQ=255)
                    cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCt                                                                                                                          >  7:361704/1‑149 (MQ=255)
                    cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGCTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCGTTTATATCAGCt                                                                                                                          >  7:169058/1‑149 (MQ=255)
                       aaaCTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACAGGGGCTTCCTTTATATCAGCTGGt                                                                                                                       >  5:204880/1‑149 (MQ=255)
                         actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGttt                                                                                                                     <  8:65548/149‑1 (MQ=255)
                           tactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTg                                                                                                                    >  8:232784/1‑148 (MQ=255)
                              tGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGc                                                                                                                >  6:252220/1‑149 (MQ=255)
                                 aGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGc                                                                                                             >  3:390264/1‑149 (MQ=255)
                                    gcAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTgg                                                                                                          <  7:30472/149‑1 (MQ=255)
                                     cAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                         <  6:356241/149‑1 (MQ=255)
                                     cAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                         >  3:337250/1‑149 (MQ=255)
                                      aGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgtg                                                                                                        >  5:399948/1‑149 (MQ=255)
                                      aGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgtg                                                                                                        >  1:178048/1‑149 (MQ=255)
                                        tCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgtgt                                                                                                       >  1:144204/1‑148 (MQ=255)
                                          tAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCAc                                                                                                    >  3:270235/1‑149 (MQ=255)
                                           agagATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGtt                                                                                                                      <  8:143517/130‑1 (MQ=255)
                                           agagATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGtt                                                                                                                      >  7:143517/1‑130 (MQ=255)
                                               aTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGCCCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCAc                                                                                               >  8:275100/1‑149 (MQ=255)
                                                tttCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCAc                                                                                               <  4:270235/148‑1 (MQ=255)
                                                   cGCCCATTGCCGGCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGc                                                                                            <  8:269776/148‑1 (MQ=255)
                                                    gCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGcc                                                                                           <  4:390264/148‑1 (MQ=255)
                                                      ccATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCtt                                                                                         >  5:13352/1‑148 (MQ=255)
                                                      ccATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCtt                                                                                         >  4:116069/1‑148 (MQ=255)
                                                      ccATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCCCCTTg                                                                                        >  6:326366/1‑149 (MQ=255)
                                                      ccATTGCCGTCCCACCAACATGCTCGTATGAGAAGACTTGTATTTGAGACGTTGATCCGGATCAATAAAAGAAAACCCCCGCCACATTCTGACCGGCCGATCCATTATATCAGAAGCTTTGCACCGGCAGGTCTCACTCCACCGcc                                                                                           >  6:344843/1‑146 (MQ=255)
                                                        aTTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGAc                                                                                      <  4:108590/149‑1 (MQ=255)
                                                            ccGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCAt                                                                                   <  8:361704/148‑1 (MQ=255)
                                                            ccGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATg                                                                                  >  5:265528/1‑149 (MQ=255)
                                                              gTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCgg                                                                                                                                             <  5:349325/88‑1 (MQ=255)
                                                              gTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCgg                                                                                                                                             >  6:349325/1‑88 (MQ=255)
                                                               tCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCTACCGCCTTGACCATGTcc                                                                               <  8:188640/149‑1 (MQ=255)
                                                                  cACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGcc                                                                                           <  4:347690/134‑1 (MQ=255)
                                                                  cACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTc                                                                             >  8:321133/1‑148 (MQ=255)
                                                                     gACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCTCTCCACCGCCTTGACCATGTCCTCGATg                                                                         <  1:223691/149‑1 (MQ=255)
                                                                      aCCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCAACGCCTTGACCATGTCCTCGATGa                                                                        >  8:153572/1‑149 (MQ=255)
                                                                        cATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGAcc                                                                      <  3:116069/149‑1 (MQ=255)
                                                                           ggTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACcttc                                                                   >  2:104410/1‑149 (MQ=255)
                                                                           ggTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGAGGCTGGTGTCACTCCACCGCCTTGACAATGTCCTCGATGACctcg                                                                   >  6:81897/1‑147 (MQ=255)
                                                                               gTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTgg                                                               >  4:224685/1‑149 (MQ=255)
                                                                               gTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTg                                                                <  5:252220/148‑1 (MQ=255)
                                                                                tATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGc                                                              <  6:265528/149‑1 (MQ=255)
                                                                                   gagaAGCCTTGTATTTTAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCAc                                                                                                    <  7:173742/108‑1 (MQ=255)
                                                                                   gagaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCACCGGCAGATTCTGACCGGGGCTTCCTTTATATAAGCTGGTTTGCAGCGGATGGTTTCAc                                                                                                    >  8:173742/1‑108 (MQ=255)
                                                                                     gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCa                                                          <  6:204880/148‑1 (MQ=255)
                                                                                       aGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCg                                                       <  8:169058/149‑1 (MQ=255)
                                                                                           ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCgggg                                                                                                                                           >  5:18694/1‑61 (MQ=255)
                                                                                           ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCgggg                                                                                                                                           <  6:18694/61‑1 (MQ=255)
                                                                                           ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTg                                                                                        >  4:51270/1‑112 (MQ=255)
                                                                                           ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTg                                                                                        <  3:51270/112‑1 (MQ=255)
                                                                                           ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAaca                                                   >  2:60409/1‑149 (MQ=255)
                                                                                                tttGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcat                                               >  7:98035/1‑148 (MQ=255)
                                                                                                  tGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcatca                                             <  2:178048/148‑1 (MQ=255)
                                                                                                   gAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATg                                           >  1:321330/1‑149 (MQ=255)
                                                                                                      tCGTTGATCCGGATCAATAAAAGAAAAGCCGCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                        <  5:326366/149‑1 (MQ=255)
                                                                                                      tCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                        <  7:232784/149‑1 (MQ=255)
                                                                                                         ttGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTg                                     >  5:90392/1‑149 (MQ=255)
                                                                                                           gATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTc                                   >  6:83561/1‑149 (MQ=255)
                                                                                                             tCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTcc                                  <  4:337250/148‑1 (MQ=255)
                                                                                                               cggtgcAATAAAAGAAAAGCCCAGGGAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATCACCTTCTTGGCGTCACCGAACAGCATCATGGTCTTGTCCAGa                               <  6:396622/144‑1 (MQ=255)
                                                                                                                ggttCAATAAAAGAAAAGCGCCGGCTGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACAGCATCATGGTCTTGTCCAGAt                              <  6:13352/146‑1 (MQ=255)
                                                                                                                ggATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGa                               >  3:63142/1‑148 (MQ=255)
                                                                                                                     aaTAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGaa                          <  3:223309/148‑1 (MQ=255)
                                                                                                                     aaTAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGaa                          <  6:399948/148‑1 (MQ=255)
                                                                                                                       tAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAg                       >  2:383439/1‑149 (MQ=255)
                                                                                                                       tAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTACCTTTATATCAGAAGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGACCTCGATGACCTTCTTGGCGTCACCCAACACCATCATGGTCTTGTCCAGATAGAACAg                       >  1:244011/1‑149 (MQ=255)
                                                                                                                        aaaagaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCGCCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                      >  1:376588/1‑149 (MQ=255)
                                                                                                                        aaaagaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                      >  6:354275/1‑149 (MQ=255)
                                                                                                                         aaagaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                      <  1:350093/148‑1 (MQ=255)
                                                                                                                          aagaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTc                    >  5:150462/1‑149 (MQ=255)
                                                                                                                           agaaaagCCCCGGCAGATTCTGTCCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCg                   <  7:321133/149‑1 (MQ=255)
                                                                                                                           agaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCg                   >  3:269053/1‑149 (MQ=255)
                                                                                                                            gaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGt                  >  5:49417/1‑149 (MQ=255)
                                                                                                                            gaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGt                  >  5:185446/1‑149 (MQ=255)
                                                                                                                                 gCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTcc             <  7:275100/149‑1 (MQ=255)
                                                                                                                                 gCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTc              >  6:239427/1‑148 (MQ=255)
                                                                                                                                  ccccGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTc              >  7:14999/1‑147 (MQ=255)
                                                                                                                                   cccggtcAGAATCTGCCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAg           >  7:304888/7‑149 (MQ=255)
                                                                                                                                    ccggtcAGAATCTGCCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAgg          >  2:89504/6‑149 (MQ=255)
                                                                                                                                         cAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCgg      <  2:144204/149‑1 (MQ=255)
                                                                                                                                          aGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGc     <  5:83561/149‑1 (MQ=255)
                                                                                                                                             ttCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATa  <  3:26765/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCCGGCATA  >  NC_002947/163833‑164121

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