Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,027,595 C→T 6.3% intergenic (+367/+14) fre → / ← fadA NAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase (thiolase I)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,027,5950CT 6.3% 200.8 / 2.5 63intergenic (+367/+14)fre/fadANAD(P)H‑flavin reductase/3‑ketoacyl‑CoA thiolase (thiolase I)
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base T (2/2);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

GTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGC  >  NC_000913/4027448‑4027741
                                                                                                                                                   |                                                                                                                                                  
gTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGACTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCa                                                                                                                                                   >  6:207483/1‑149 (MQ=255)
    cTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATGCGGCGTGAACGCCTTATCCg                                                                                                                                               <  6:226765/149‑1 (MQ=18)
      aCAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                              >  6:118255/1‑148 (MQ=255)
       cAAATAGTCCGAACCCTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGAAGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATGCGGCGTGAACGCCTTATCCggtc                                                                                                                                            <  1:21298/149‑3 (MQ=18)
        aaaTAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                             >  3:164941/1‑147 (MQ=255)
            aGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTTGTGCATGATGCCTGATGCGACGCTGGCTATTCTTATAAGGCCTACATTGGTGCAGGATCTTTAGACAGGATAAGGCGTTCACGCCGCATCCGGCAAGTgg                                                                                                                                       >  4:41391/1‑149 (MQ=18)
               ccGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGtt                                                                                                                                     >  7:86033/1‑148 (MQ=255)
                 gAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTaaa                                                                                                                                  <  7:187131/149‑1 (MQ=255)
                  caCCGTAGGCCGGATACGGCGTTCACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAAc                                                                                                                                 <  3:171403/148‑1 (MQ=255)
                  aaCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAAc                                                                                                                                 <  4:126479/149‑1 (MQ=255)
                  aaCCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTAGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAAc                                                                                                                                 <  2:162106/149‑1 (MQ=255)
                      gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGATGGCGCGTCTTATCAGGCATACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGc                                                                                                                             >  4:28238/1‑149 (MQ=255)
                       tAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCt                                                                                                                            >  2:66404/1‑149 (MQ=255)
                          gCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCa                                                                                                                          >  3:122889/1‑148 (MQ=255)
                           ccGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCaaa                                                                                                                        >  1:70419/1‑149 (MQ=255)
                               aTAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACAcc                                                                                                                    <  3:120802/149‑1 (MQ=255)
                                  aGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTc                                                                                                                 <  1:86316/149‑1 (MQ=255)
                                   ggCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCg                                                                                                                >  6:79340/1‑149 (MQ=255)
                                   ggCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCg                                                                                                                >  8:82219/1‑149 (MQ=255)
                                      gTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCaa                                                                                                             <  7:224032/149‑1 (MQ=255)
                                       tttACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCaa                                                                                                             >  6:165886/1‑148 (MQ=255)
                                          aCGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATAcc                                                                                                         >  7:210113/1‑149 (MQ=255)
                                              cgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAACCACCGTCGCAATACCCTGa                                                                                                     >  1:190772/1‑149 (MQ=255)
                                              cgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGa                                                                                                     <  7:99557/149‑1 (MQ=255)
                                              cgcATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCAGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGa                                                                                                     >  6:102501/1‑149 (MQ=255)
                                                 aTCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGAcc                                                                                                   >  3:29213/1‑148 (MQ=255)
                                                  tCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCAGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCa                                                                                                 <  5:102501/149‑1 (MQ=255)
                                                    cGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGa                                                                                               <  4:29213/149‑1 (MQ=255)
                                                     ggCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGAc                                                                                              <  3:63295/149‑1 (MQ=255)
                                                     ggCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGAc                                                                                              <  7:82219/149‑1 (MQ=255)
                                                         aTTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGAt                                                                                          <  6:212114/149‑1 (MQ=255)
                                                            ggTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATaca                                                                                       >  8:60840/1‑149 (MQ=255)
                                                                cATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAtc                                                                                   <  8:210113/149‑1 (MQ=255)
                                                                  tgatgCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGa                                                                                           >  6:46824/1‑139 (MQ=255)
                                                                  tgatgCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGa                                                                                           <  5:46824/139‑1 (MQ=255)
                                                                     tgCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGc                                                                              <  4:164941/149‑1 (MQ=255)
                                                                       ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCa                                                                            <  2:70419/149‑1 (MQ=255)
                                                                        cTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAtcgt                                                                                 <  3:10603/143‑1 (MQ=255)
                                                                        cTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAtcgt                                                                                 >  4:10603/1‑143 (MQ=255)
                                                                               gACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCCAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGCCCAAAc                                                                    >  4:23620/1‑149 (MQ=255)
                                                                               gACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACAGGATAATGCTTTCACGCCGCATCCGGCAAGTGGTTAAACCAGCTCAAACACATTCGCAATACCATGACCCAGACCGATCCACATCGTCGCAAGACAAAAc                                                                    >  2:62194/1‑149 (MQ=38)
                                                                                   cTGGTGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGaa                                                                >  4:6321/1‑149 (MQ=255)
                                                                                          gTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCttt                                                         <  5:79340/149‑1 (MQ=255)
                                                                                          gTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCttt                                                         <  3:6321/149‑1 (MQ=255)
                                                                                                tCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTc                                                   <  7:60840/149‑1 (MQ=255)
                                                                                                    gCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAt                                                                                    >  2:22010/1‑112 (MQ=255)
                                                                                                    gCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAt                                                                                    <  1:22010/112‑1 (MQ=255)
                                                                                                    gCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACACAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCAt                                                >  7:85278/1‑148 (MQ=255)
                                                                                                     ccTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCa                                              <  5:118255/149‑1 (MQ=255)
                                                                                                        aCATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGAt                                           >  4:159964/1‑149 (MQ=255)
                                                                                                         cATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTTACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGAtt                                          <  3:41391/149‑1 (MQ=255)
                                                                                                                 ccGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCg                                  >  4:46382/1‑149 (MQ=255)
                                                                                                                        ggTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGc                            <  3:28238/148‑1 (MQ=255)
                                                                                                                         gTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCt                           >  8:106075/1‑148 (MQ=255)
                                                                                                                           ttgccGGATGCGGCGTGAACGCCTTATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGa                         >  6:11114/4‑148 (MQ=255)
                                                                                                                                ggATAAGGAGTTCACGCAGCATCCGCCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGCTAAACATCGTCGCCAGACCAAACTGAACGAGTATGCGTTCCATCAGATTCAGCAGCGTGCTGCTGATAcgcg                   >  2:192172/1‑149 (MQ=255)
                                                                                                                                 gATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATAcgcg                   >  7:217586/1‑148 (MQ=255)
                                                                                                                                  aTAAGGCGTTCACGCCGTATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCa                 >  8:72786/1‑149 (MQ=255)
                                                                                                                                    aaGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAtcgt                                                                                 >  6:111960/1‑83 (MQ=255)
                                                                                                                                    aaGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACAtcgt                                                                                 <  5:111960/83‑1 (MQ=255)
                                                                                                                                     aGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCg              >  6:17847/1‑149 (MQ=255)
                                                                                                                                     aGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCg              >  3:96198/1‑149 (MQ=255)
                                                                                                                                        cGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGaa           >  8:233817/1‑149 (MQ=255)
                                                                                                                                               gccgcATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATCCCCTGCCCCAGCCCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAAcc      >  1:90619/1‑147 (MQ=255)
                                                                                                                                                 cgcATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGc  >  5:74521/1‑149 (MQ=255)
                                                                                                                                                 cgcATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGc  <  3:46382/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTGCCTACAAATAGTCCGAACCGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATTGGTGCATGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACATTGGTGCCGGATCGGTAGACCGGATAAGGCGTTCACGCCGCATCCGGCAAGTGGTTAAACCCGCTCAAACACCGTCGCAATACCCTGACCCAGACCGATACACATCGTCGCCAGACCAAACTGAACGTCTTTGCGTTCCATCAGATTCAGCAGCGTGGTGCTGATACGCGCACCGGAACAACCCAGC  >  NC_000913/4027448‑4027741

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: