Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,457,953 T→C 64.5% K403E (AAA→GAA)  fadJ ← enoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,457,9530TC64.5% 156.1 / 145.3 155K403E (AAA→GAA) fadJenoyl‑CoA hydratase/epimerase and isomerase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (24/31);  new base C (58/42);  total (82/73)
Fisher's exact test for biased strand distribution p-value = 9.52e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.17e-01

Only 100 of 156 total aligned reads displayed.

GTGAGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTGTTTGTCAC  >  NC_000913/2457855‑2458051
                                                                                                  |                                                                                                  
gtgAGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTc                                                                                                   <  1:2423501/100‑1 (MQ=255)
 tgAGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCa                                                                                                  >  2:2591271/1‑100 (MQ=255)
 tgAGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCa                                                                                                  >  1:1548461/1‑100 (MQ=255)
  gAGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCaa                                                                                                 <  2:11012/100‑1 (MQ=255)
  gAGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCaa                                                                                                 <  2:3544985/100‑1 (MQ=255)
   aGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAAt                                                                                                >  2:1737571/1‑100 (MQ=255)
   aGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAAt                                                                                                >  1:2947451/1‑100 (MQ=255)
      ggcgATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCg                                                                                             >  1:613331/1‑100 (MQ=255)
      ggcgATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCg                                                                                             <  1:458675/100‑1 (MQ=255)
       gcgATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCga                                                                                            >  1:785398/1‑100 (MQ=255)
         gATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCgaga                                                                                          >  1:544228/1‑100 (MQ=255)
         gATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCgaga                                                                                          >  1:2699452/1‑100 (MQ=255)
           tatCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGAtt                                                                                        >  2:1172955/1‑100 (MQ=255)
             tCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCAACCATCTGTTGTTCCAATTCGAGAtttt                                                                                      >  1:2997121/1‑100 (MQ=255)
              cACCAATCGGTAAAGATGACGTCTTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTc                                                                                     >  2:2399646/1‑100 (MQ=255)
              cACCAATAGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTGCCGCCACCATCTGTTGTTCCAATTCGAGATTTTc                                                                                     >  2:2174147/1‑100 (MQ=255)
               aCCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCa                                                                                    >  1:3570955/1‑100 (MQ=255)
               aCCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCACTTCGAGAttttga                                                                                    <  1:2809775/100‑3 (MQ=255)
                 cAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCaaa                                                                                  <  2:1544951/100‑1 (MQ=255)
                      ggTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCg                                                                             >  2:1095069/1‑100 (MQ=255)
                      ggTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCg                                                                             >  1:3856123/1‑100 (MQ=255)
                       gTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGc                                                                            <  1:3636709/100‑1 (MQ=255)
                        tAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCt                                                                           >  2:104624/1‑100 (MQ=255)
                         aaaGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCtt                                                                          >  2:1411663/1‑100 (MQ=255)
                            gATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCaa                                                                       >  1:1953668/1‑100 (MQ=255)
                              tGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCaata                                                                     <  1:3568266/100‑1 (MQ=255)
                               gACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCaataa                                                                    >  2:2858142/1‑100 (MQ=255)
                                 cGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATc                                                                  <  2:3330754/100‑1 (MQ=255)
                                 cGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATc                                                                  <  2:2043760/100‑1 (MQ=255)
                                  gTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCa                                                                 <  1:3040403/100‑1 (MQ=255)
                                  gTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCa                                                                 >  2:2343827/1‑100 (MQ=255)
                                  gTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCa                                                                 <  1:2527503/100‑1 (MQ=255)
                                   tATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAg                                                                >  2:1067681/1‑100 (MQ=255)
                                    aTTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGa                                                               >  1:3907584/1‑100 (MQ=255)
                                    aTTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAg                                                                >  2:279041/1‑99 (MQ=255)
                                    aTTCGAAGCAAAGATGCTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCCTCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCACTAATCAGa                                                               >  2:830961/1‑100 (MQ=255)
                                     ttCAAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAAATTCCGCCACCATCTGTTGTTCCAATTCGAGAATTTCAAACACCGCTTCAATAATCAGa                                                               >  2:1898012/1‑99 (MQ=255)
                                        gAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATcgc                                                           >  1:597083/1‑100 (MQ=255)
                                        gAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATcgc                                                           >  2:2254183/1‑100 (MQ=255)
                                         aaGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATcgcg                                                          >  1:974904/1‑100 (MQ=255)
                                          aGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGa                                                         >  2:1460129/1‑100 (MQ=255)
                                            cAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATg                                                       >  1:1641097/1‑100 (MQ=255)
                                            cAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATg                                                       >  1:3598246/1‑100 (MQ=255)
                                            cAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATg                                                       >  2:2731405/1‑100 (MQ=255)
                                            cAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATg                                                       >  2:2647025/1‑100 (MQ=255)
                                             aaaGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATgg                                                      <  2:785392/100‑1 (MQ=255)
                                               aGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGc                                                    <  1:3239039/100‑1 (MQ=255)
                                               aGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGc                                                    <  1:682337/100‑1 (MQ=255)
                                                gATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCa                                                   >  1:3553835/1‑100 (MQ=255)
                                                 aTGGTATGAGCGTCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCa                                                   <  1:1424044/99‑1 (MQ=255)
                                                     tATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGcc                                              >  2:2496115/1‑100 (MQ=255)
                                                      aTGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCcg                                             >  1:717545/1‑100 (MQ=255)
                                                          gcggcgCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGAt                                         <  1:369946/100‑1 (MQ=255)
                                                          gcggcgCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGAt                                         <  1:1278183/100‑1 (MQ=255)
                                                           cggcgCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATa                                        <  1:2388277/100‑1 (MQ=255)
                                                             gcgcAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGt                                      >  1:1730442/1‑100 (MQ=255)
                                                              cgcAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTc                                     <  2:1756070/100‑1 (MQ=255)
                                                              cgcAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTc                                     <  2:1669021/100‑1 (MQ=255)
                                                              cgcAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTc                                     <  2:4134658/100‑1 (MQ=255)
                                                              cgcAATTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTc                                     <  1:1319150/100‑1 (MQ=255)
                                                               gcAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTcc                                    >  1:2183778/1‑100 (MQ=255)
                                                                 aaTTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCcgt                                  <  2:1848849/100‑1 (MQ=255)
                                                                 aaTTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCcgt                                  <  1:2255495/100‑1 (MQ=255)
                                                                 aaTTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTc                                     <  2:1386379/97‑1 (MQ=255)
                                                                  aTTTTGCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCcgtc                                 <  1:110288/100‑1 (MQ=255)
                                                                      tGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTc                             <  2:573838/100‑1 (MQ=255)
                                                                       gCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTcc                            >  2:3159609/1‑100 (MQ=255)
                                                                       gCTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTcc                            >  2:339793/1‑100 (MQ=255)
                                                                        cTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCg                           >  1:3504614/1‑100 (MQ=255)
                                                                        cTCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCg                           >  1:823020/1‑100 (MQ=255)
                                                                         tCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCgg                          <  2:3540286/100‑1 (MQ=255)
                                                                         tCAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCgg                          <  1:1023039/100‑1 (MQ=255)
                                                                         ncAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCgg                          >  1:3654848/2‑100 (MQ=255)
                                                                          cAACTTCCGCCACCATCTGTTGTTTCAATTCGAGACTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGga                         >  2:376837/1‑100 (MQ=255)
                                                                          cAACTTCCGCCACCATCTGTTGTTCCAATTCGCGACTTTCAAACACCGCCTCCATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTcc                            >  2:170857/1‑97 (MQ=255)
                                                                          cAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGga                         >  2:704034/1‑100 (MQ=255)
                                                                          cAACTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCAATGGGCAAAGCCGCGATAGTCCGTCGTTCCgg                          >  2:1061098/1‑99 (MQ=255)
                                                                          cAACTTCCGCCACCATCTGTCGTTTCAATCCGAGATTTTCAACCACCGCTTCAATAATCACATCGCGCTGGGCAACGCCGCGATAGTCCCTCCTTCCgg                          >  2:799841/1‑99 (MQ=255)
                                                                             cTTCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGAt                      <  1:1092245/100‑1 (MQ=255)
                                                                              ttCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGAtt                     <  1:4083650/100‑1 (MQ=255)
                                                                              ttCCGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGAtt                     <  2:590528/100‑1 (MQ=255)
                                                                                ccGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTaa                   <  2:2934785/100‑1 (MQ=255)
                                                                                 cGCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGCGATTaa                   >  1:3447686/1‑99 (MQ=255)
                                                                                  gCCACCATCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGCCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATg                 >  2:2278643/1‑100 (MQ=255)
                                                                                   ccaccaTCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGc                >  1:335672/1‑100 (MQ=255)
                                                                                    ctccaTCTGTTGTTCCAATTCGAGATTTTCTAACAGCGCTTTAATAATCAGGTGGCGATGGACAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGcc               <  2:3127610/98‑1 (MQ=255)
                                                                                    caccaTCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGcc               <  1:328428/100‑1 (MQ=255)
                                                                                      ccaTCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAg             <  1:893886/100‑1 (MQ=255)
                                                                                      ccaTCTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAg             <  2:3614847/100‑1 (MQ=255)
                                                                                       caTCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGc            <  2:1743483/100‑1 (MQ=255)
                                                                                       caTCTGTTGTTCCAATTCGAGATTTTCAAACACCGCGTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGc            >  1:2913957/1‑100 (MQ=255)
                                                                                         tCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTg          <  1:104625/100‑1 (MQ=255)
                                                                                         tCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTg          <  1:1461113/100‑1 (MQ=255)
                                                                                          cTGTTGTTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTGt         <  1:1999237/100‑1 (MQ=255)
                                                                                           tgttgtTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCTTCGTTCCGGAGATTAATGCCAGCTGtt        >  2:3665286/1‑100 (MQ=255)
                                                                                           tgttgtTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTGtt        >  2:696766/1‑100 (MQ=255)
                                                                                           tgttgtTCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTGtt        >  1:3943575/1‑100 (MQ=255)
                                                                                            gttgttCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTACGTCGTTCCGGAGATTAATGCCAGCTGttt       >  2:2767941/1‑100 (MQ=255)
                                                                                             ttgttCCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTGTTTg      <  1:3159628/100‑1 (MQ=255)
                                                                                               gttCCAATTCGAGATTTTCAAACACCGCTTCAAGAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTAGTTCCGGAGATTAATGCCAGCTGTTTGTc    >  2:2490825/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
GTGAGCGGCGATATCACCAATCGGTAAAGATGACGTATTCGAAGCAAAGATGGTATGAGCGGCGCAATTTTGCTCAACTTCCGCCACCATCTGTTGTTTCAATTCGAGATTTTCAAACACCGCTTCAATAATCAGATCGCGATGGGCAAAGCCGCGATAGTCCGTCGTTCCGGAGATTAATGCCAGCTGTTTGTCAC  >  NC_000913/2457855‑2458051

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