Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 163,973 T→A 6.5% 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 6.5% 185.6 / 2.6 61intergenic (‑140/+47)scpC/pntBpropionyl‑CoA:succinate CoA transferase/pyridine nucleotide transhydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base A (2/2);  total (27/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 = 9.64e-01

ACATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCC  >  NC_002947/163830‑164115
                                                                                                                                                |                                                                                                                                               
aCATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTg                                                                                                                                             <  7:263588/148‑1 (MQ=255)
aCATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTg                                                                                                                                             >  8:263588/1‑148 (MQ=255)
  aTGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGAcc                                                                                                                                          >  7:251539/1‑149 (MQ=255)
      aTTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTAGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCgggg                                                                                                                                      >  5:297837/1‑149 (MQ=255)
         gTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTGCCGGGGCtt                                                                                                                                  >  5:129122/1‑149 (MQ=255)
          ttATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTGCCGGGGCttt                                                                                                                                 <  3:302507/149‑2 (MQ=255)
              ctcGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCtt                                                                                                                              >  2:90352/1‑148 (MQ=255)
                 gggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGTCAGAATCTGCCGGGGCttttcttt                                                                                                                            <  6:129122/147‑7 (MQ=255)
                 gggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCtt                                                                                                                              <  2:249890/145‑1 (MQ=255)
                  ggCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCAGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTAatat                                                                                                                         >  6:221203/1‑149 (MQ=255)
                       cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAg                                                                                                                      <  7:214797/147‑1 (MQ=255)
                       cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCt                                                                                                                    >  6:49216/1‑149 (MQ=255)
                       cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCt                                                                                                                    <  5:240222/149‑1 (MQ=255)
                       cTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCt                                                                                                                    >  1:181841/1‑149 (MQ=255)
                         gAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTgg                                                                                                                  >  6:240636/1‑149 (MQ=255)
                            actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACTATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGttt                                                                                                               >  2:215832/1‑149 (MQ=255)
                            actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACTATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATAAGCTGGttt                                                                                                               >  2:211510/1‑149 (MQ=255)
                            actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGttt                                                                                                               >  3:277271/1‑149 (MQ=255)
                            actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGttt                                                                                                               <  4:186562/149‑1 (MQ=255)
                            actactGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGtt                                                                                                                >  5:55034/1‑148 (MQ=255)
                                   cAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCg                                                                                                         <  7:280784/148‑1 (MQ=255)
                                        cAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                   <  6:184921/149‑1 (MQ=255)
                                        cAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGgt                                                                                                   <  5:49216/149‑1 (MQ=255)
                                             tAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCAc                                                                                              >  5:2649/1‑149 (MQ=255)
                                                agaTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                            <  8:281912/148‑1 (MQ=255)
                                                agaTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTc                                                                                            >  7:281912/1‑148 (MQ=255)
                                                 gaTTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCa                                                                                          <  4:53946/149‑1 (MQ=255)
                                                   tttCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTTGTGTCACTCCAcc                                                                                        <  6:297837/149‑1 (MQ=255)
                                                            ttGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTTCAGCGGCTGGTGTCACTCCACCGCCTTGAcc                                                                               <  6:55034/149‑1 (MQ=255)
                                                               ccGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATg                                                                            >  7:15660/1‑149 (MQ=255)
                                                                                  gTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTg                                                          <  2:168803/148‑1 (MQ=255)
                                                                                   tATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTatat                                                                                                                         >  3:154650/1‑84 (MQ=255)
                                                                                   tATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTatat                                                                                                                         <  4:154650/84‑1 (MQ=255)
                                                                                        gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCAc                                                   <  2:181841/149‑1 (MQ=255)
                                                                                        gaAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCAc                                                   >  7:150774/1‑149 (MQ=255)
                                                                                             cTTGTATTTGAGTCGTTGATCCGTATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGTTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAac                                              <  8:251539/149‑1 (MQ=255)
                                                                                              ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAaca                                             >  2:250279/1‑149 (MQ=255)
                                                                                              ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAaca                                             >  2:265306/1‑149 (MQ=255)
                                                                                              ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAaca                                             >  8:186145/1‑149 (MQ=255)
                                                                                              ttGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACAGGGGCTTCCTTTATATAAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCAt                                                                             >  3:162412/1‑117 (MQ=255)
                                                                                                   tttGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATGAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcat                                         >  3:287162/1‑148 (MQ=255)
                                                                                                        gTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGt                                   <  3:302889/149‑1 (MQ=255)
                                                                                                         tCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                  <  5:240636/149‑1 (MQ=255)
                                                                                                         tCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                  <  2:154219/149‑1 (MQ=255)
                                                                                                          cGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTc                                  <  8:150774/148‑1 (MQ=255)
                                                                                                           gTTTATCCTGATCAATAAAAGAATAGCCCGGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCTCGGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCtt                                <  5:221203/149‑1 (MQ=255)
                                                                                                           gTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCtt                                <  6:2649/149‑1 (MQ=255)
                                                                                                                  cGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGa                         <  1:285028/149‑1 (MQ=255)
                                                                                                                   ggATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCTGAt                        <  2:126173/149‑1 (MQ=255)
                                                                                                                   ggATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGAt                        <  7:186145/149‑1 (MQ=255)
                                                                                                                   ggATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGAt                        <  4:277271/149‑1 (MQ=255)
                                                                                                                    gATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATa                       <  1:211510/149‑1 (MQ=255)
                                                                                                                    gATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGAAGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATa                       <  1:215832/149‑1 (MQ=255)
                                                                                                                         aTAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGaa                    >  5:11587/1‑147 (MQ=255)
                                                                                                                           aaaagaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGt                >  1:145024/1‑149 (MQ=255)
                                                                                                                            aaagaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACcatc                                        <  2:106209/124‑1 (MQ=255)
                                                                                                                              agaaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCg             <  2:145024/149‑1 (MQ=255)
                                                                                                                                aaaagCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGtt           <  2:98836/149‑1 (MQ=255)
                                                                                                                                   agCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTc        <  1:90352/149‑1 (MQ=255)
                                                                                                                                     ccccGGTCAGAATCTGCCGGGGCTTTTCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTc        >  8:40468/1‑148 (MQ=255)
                                                                                                                                         ggCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGcc  >  1:179367/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                               
ACATGGATTGTTATCTCGGGCTACTGAAACTACTGCAGCGCAGTCTAGAGATTTCGCCCATTGCCGTCCCACGACCATGGTCGTATGAGAAGCCTTGTATTTGAGTCGTTGATCCGGATCAATAAAAGAAAAGCCCCGGCAGATTCTGACCGGGGCTTCCTTTATATCAGCTGGTTTGCAGCGGCTGGTGTCACTCCACCGCCTTGACCATGTCCTCGATGACCTTCTTGGCGTCACCGAACACCATCATGGTCTTGTCCAGATAGAACAGTTCGTTGTCCAGGCC  >  NC_002947/163830‑164115

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: