Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 353,717 C→G 6.1% R338G (CGC→GGC)  prpE → propionate‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913353,7170CG 6.1% 213.4 / 6.1 66R338G (CGC→GGC) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (25/36);  new base G (2/2);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

TTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACAACTA  >  NC_000913/353572‑353862
                                                                                                                                                 |                                                                                                                                                 
ttCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTcgcg                                                                                                                                                <  3:155682/149‑1 (MQ=255)
   gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTg                                                                                                                                              <  6:50013/148‑1 (MQ=255)
   gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGc                                                                                                                                             >  7:266817/1‑149 (MQ=255)
   gTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGc                                                                                                                                             >  3:56129/1‑149 (MQ=255)
     atatCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCt                                                                                                                                            >  2:78757/1‑148 (MQ=255)
      tatCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTg                                                                                                                                           >  3:159341/1‑148 (MQ=255)
      tatCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTg                                                                                                                                           >  3:95692/1‑148 (MQ=255)
          gTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGaaaa                                                                                                                                       >  6:6357/1‑148 (MQ=255)
           tttACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGaaaaaa                                                                                                                                     >  7:326048/1‑149 (MQ=255)
              aCGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTc                                                                                                                                  >  8:118593/1‑149 (MQ=255)
              aCGCGCCGCTGCTGGCGGGGATGGCGACTAACTTATACGAAGTATTGCCGACATGGCCGGACTGCGGCGTGAGGTGGACAATCGTCGAGAAATATCAGGTTAGCCCGATGATCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAAtta                                                                                                                                  >  4:271449/1‑148 (MQ=255)
                gcgcCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                                                >  4:207110/1‑149 (MQ=255)
                gcgcCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTAGCGTGCTGAAAAAATTccc                                                                                                                                >  3:208743/1‑149 (MQ=255)
                gcgcCGCTGCTGGCGGGGATGGCGACTAGCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTccc                                                                                                                                <  8:92207/149‑1 (MQ=255)
                   ccgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGCCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGTTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTAc                                                                                                                             >  6:290491/1‑149 (MQ=255)
                   ccgcTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTAc                                                                                                                             <  5:217140/149‑1 (MQ=255)
                              ggggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTc                                                                                                                  <  3:102359/149‑1 (MQ=255)
                               gggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCg                                                                                                                 <  4:56129/149‑1 (MQ=255)
                               gggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCg                                                                                                                 >  4:33467/1‑149 (MQ=255)
                                ggATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGc                                                                                                                <  1:31561/149‑1 (MQ=255)
                                     gCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAaca                                                                                                           >  5:81834/1‑149 (MQ=255)
                                          tATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtc                                                                                                      <  3:126249/149‑1 (MQ=255)
                                            tCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctc                                                                                                    >  4:168778/1‑149 (MQ=255)
                                            tCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctc                                                                                                    >  4:164544/1‑149 (MQ=255)
                                              gTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctctc                                                                                                  >  2:7277/1‑149 (MQ=255)
                                              gTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctct                                                                                                   <  7:98332/148‑1 (MQ=255)
                                                 tACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCtcgt                                                                                                <  3:207110/148‑1 (MQ=255)
                                                   cGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAac                                                                                                            >  2:50180/1‑134 (MQ=255)
                                                   cGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAac                                                                                                            <  1:50180/134‑1 (MQ=255)
                                                    gAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTc                                                                                                                  >  5:319936/1‑127 (MQ=255)
                                                    gAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTc                                                                                                                  <  6:319936/127‑1 (MQ=255)
                                                    gAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCt                                                                                            >  7:40207/1‑149 (MQ=255)
                                                        gATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGaa                                                                                        <  1:239030/149‑1 (MQ=255)
                                                             ccGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGct                                                                                   <  1:286802/149‑1 (MQ=255)
                                                              cGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGctc                                                                                  >  4:67043/1‑149 (MQ=255)
                                                                  cTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATc                                                                              <  5:11582/149‑1 (MQ=255)
                                                                       cGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctggct                                                                         >  7:15287/1‑149 (MQ=255)
                                                                        ggACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctggctg                                                                        <  8:9691/149‑1 (MQ=255)
                                                                          aCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGga                                                                      <  3:168778/149‑1 (MQ=255)
                                                                          aCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTGGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGga                                                                      <  3:164544/149‑1 (MQ=255)
                                                                          aCTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATctggctg                                                                        >  3:120914/1‑147 (MQ=255)
                                                                           cTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGgag                                                                     <  3:271449/149‑1 (MQ=255)
                                                                           cTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGgag                                                                     <  7:118593/149‑1 (MQ=255)
                                                                              cggcgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGaa                                                                   <  1:78757/148‑1 (MQ=255)
                                                                              cggcgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAAc                                                                  >  1:56782/1‑149 (MQ=255)
                                                                                gcgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCg                                                                <  8:326048/149‑1 (MQ=255)
                                                                                 cgTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGc                                                               <  4:159341/149‑1 (MQ=255)
                                                                                         ggACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAg                                                       >  1:130052/1‑149 (MQ=255)
                                                                                               tcgtcgAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGAcc                                                 <  6:49102/149‑1 (MQ=255)
                                                                                                 gtcgAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCt                                                                                            <  6:276071/104‑1 (MQ=255)
                                                                                                    gagaAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCa                                             <  2:130052/148‑1 (MQ=255)
                                                                                                         atatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTggg                                       <  1:39419/149‑1 (MQ=255)
                                                                                                         atatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCGCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTgg                                        <  3:33467/148‑1 (MQ=255)
                                                                                                          tatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAtctc                                                                                                    <  2:143217/87‑1 (MQ=255)
                                                                                                          tatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGt                                      <  1:7277/149‑1 (MQ=255)
                                                                                                          tatCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGt                                      <  8:266817/149‑1 (MQ=255)
                                                                                                            tCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGt                                      >  3:178241/1‑147 (MQ=255)
                                                                                                              aGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAg                                   <  4:120914/148‑1 (MQ=255)
                                                                                                                   aGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATa                              <  3:67043/148‑1 (MQ=255)
                                                                                                                     ccGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATGCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGc                           <  5:6357/149‑1 (MQ=255)
                                                                                                                         aTGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATGCGCAAACACGATCTCTCGTCGCGGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGa                       <  4:95692/149‑1 (MQ=255)
                                                                                                                              ctcAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGAACTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGATAACCGCTGGACGAGCCGACCGCCACTTGGATGAGCAATACGCTGGATGTGc                  >  4:184332/1‑149 (MQ=255)
                                                                                                                                 aGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGc                  <  7:37558/146‑1 (MQ=255)
                                                                                                                                      cGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATc          <  7:72484/149‑1 (MQ=255)
                                                                                                                                         ccGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGAc       <  8:40207/149‑1 (MQ=255)
                                                                                                                                           gCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACaa     <  3:205415/149‑1 (MQ=255)
                                                                                                                                              aTTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACAacta  >  5:217880/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
TTCGTATATCGTTTACGCGCCGCTGCTGGCGGGGATGGCGACTATCGTTTACGAAGGATTGCCGACCTGGCCGGACTGCGGCGTGTGGTGGACAATCGTCGAGAAATATCAGGTTAGCCGGATGTTCTCAGCGCCGACCGCCATTCGCGTGCTGAAAAAATTCCCTACCGCTGAAATTCGCAAACACGATCTCTCGTCGCTGGAAGTGCTCTATCTGGCTGGAGAACCGCTGGACGAGCCGACCGCCAGTTGGGTGAGCAATACGCTGGATGTGCCGGTCATCGACAACTA  >  NC_000913/353572‑353862

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: