Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,930,115 G→A G308G (GGC→GGT fucO ← L‑1,2‑propanediol oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,930,1150GA96.6% 1157.0 / 11.0 355G308G (GGC→GGTfucOL‑1,2‑propanediol oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (4/8);  new base A (167/175);  total (172/183)
Fisher's exact test for biased strand distribution p-value = 3.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

Only 100 of 368 total aligned reads displayed.

TGGCGGAATACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACGCCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCGCAAC  >  NC_000913/2930025‑2930204
                                                                                          |                                                                                         
  gCGGAATACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCa                                                                                       <  1:5652213/93‑1 (MQ=255)
   cGGAATACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCAt                                                                                      <  1:2230170/93‑1 (MQ=255)
     gAATACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATaa                                                                                    >  1:5610101/1‑93 (MQ=255)
      aaTACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAAc                                                                                   >  2:2780267/1‑93 (MQ=255)
       aTACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACCCCCATAAcc                                                                                  >  2:7224122/1‑92 (MQ=255)
       aTACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAAcg                                                                                  >  2:6722029/1‑93 (MQ=255)
        tACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAAcgc                                                                                 >  1:7694987/1‑93 (MQ=255)
        tACCGACATCACGGTTGAGAGCAAAAACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAAcgc                                                                                 >  2:2936569/1‑93 (MQ=255)
          ccGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAAcgcgc                                                                               <  2:9764180/93‑1 (MQ=255)
          ccGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAAcgcgc                                                                               >  1:6894166/1‑93 (MQ=255)
            gACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAAcgcgcgc                                                                             <  2:7910785/93‑1 (MQ=255)
              cATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCTCGa                                                                           >  1:1927629/1‑93 (MQ=255)
               aTCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGat                                                                          >  2:9580606/1‑93 (MQ=255)
                tCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACGCCCATAACGCGCGCGata                                                                         <  2:7286712/93‑1 (MQ=255)
                tCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGata                                                                         >  1:469139/1‑93 (MQ=255)
                 ctccgTTGAGAGCAAACACCGCTTCAAAAGCGGCATTACGCTCATTTTCCAGGCTCATCCCTTCCCCTTTCACCCCCATAGCGCGCGCGatat                                                                        <  2:8139995/89‑1 (MQ=255)
                 cACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGatat                                                                        <  1:3419291/93‑1 (MQ=255)
                 cACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGatat                                                                        <  2:4648060/93‑1 (MQ=255)
                  aCGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATc                                                                       <  2:473555/93‑1 (MQ=255)
                     gTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATcgcg                                                                    >  2:8233800/1‑93 (MQ=255)
                     gTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATcgcg                                                                    >  2:8729323/1‑93 (MQ=255)
                       tGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGt                                                                  <  1:5520608/93‑1 (MQ=255)
                         agagCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTAc                                                                <  1:7568702/93‑1 (MQ=255)
                         agagCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTAc                                                                <  1:5344599/93‑1 (MQ=255)
                          gagCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACt                                                               <  2:6164962/93‑1 (MQ=255)
                          gagCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACt                                                               >  2:314653/1‑93 (MQ=255)
                            gCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCAATAACGCGCGCGATATCGCGGTACTtc                                                             >  1:2744749/1‑93 (MQ=255)
                             cAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATATCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTtct                                                            >  1:1314474/1‑93 (MQ=255)
                             cAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTtct                                                            <  1:7794438/93‑1 (MQ=255)
                             cAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTtct                                                            >  1:3842259/1‑93 (MQ=255)
                             cAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTtct                                                            <  1:885278/93‑1 (MQ=255)
                              aaaCACCGCTTGAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACGCCCATAACGCGCGCGATATCGCGGTACTtctc                                                           <  1:8036801/93‑1 (MQ=255)
                              aaaCACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTtctc                                                           <  1:7852265/93‑1 (MQ=255)
                              aaaCACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTtctc                                                           >  2:2549650/1‑93 (MQ=255)
                              aaaCACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTtctc                                                           >  2:4846179/1‑93 (MQ=255)
                               aaCACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCa                                                          >  2:1712245/1‑93 (MQ=255)
                               aaCACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCa                                                          >  1:6392986/1‑93 (MQ=255)
                               aaCACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCa                                                          >  2:8564992/1‑93 (MQ=255)
                                acacCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCAc                                                         >  2:1515861/1‑93 (MQ=255)
                                acacCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCAc                                                         <  1:1767355/93‑1 (MQ=255)
                                 cacCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCAcc                                                        <  1:8332222/93‑1 (MQ=255)
                                   ccGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCgg                                                      <  1:5013430/93‑1 (MQ=255)
                                     gCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTa                                                    >  1:3332503/1‑93 (MQ=255)
                                     ccTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTa                                                    <  2:1958232/92‑1 (MQ=255)
                                       ttCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACGCCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTaaa                                                  <  2:9091770/93‑1 (MQ=255)
                                       ttCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTaaa                                                  <  2:3930085/93‑1 (MQ=255)
                                       ttCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTaaa                                                  <  1:8381238/93‑1 (MQ=255)
                                       ttCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTaaa                                                  <  2:7796201/93‑1 (MQ=255)
                                        tCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAg                                                 >  2:4247249/1‑93 (MQ=255)
                                        tCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAg                                                 <  2:1301149/93‑1 (MQ=255)
                                          aaCAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTc                                               >  2:4182827/1‑93 (MQ=255)
                                          aaCAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTc                                               >  2:5546417/1‑93 (MQ=255)
                                           aCAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCa                                              <  1:4856122/93‑1 (MQ=255)
                                           aCAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCa                                              <  1:3018844/93‑1 (MQ=255)
                                            cAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAg                                             >  2:1153871/1‑93 (MQ=255)
                                             cgCGGCATTACGCGCCTCTTCCAGGCTCATTCCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGc                                            <  1:8646266/92‑1 (MQ=255)
                                                ggCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGtt                                         <  2:7662433/93‑1 (MQ=255)
                                                   aTTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTACAGTCAGCGTtata                                      >  1:1299817/1‑93 (MQ=255)
                                                   aTTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTtata                                      >  1:2442420/1‑93 (MQ=255)
                                                    ttACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTCCTTCTCACCGGTAAAGTCAGCGTTATaa                                     <  2:3791876/93‑1 (MQ=255)
                                                    ttACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATaa                                     <  1:9564964/93‑1 (MQ=255)
                                                    ttACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATaa                                     <  2:6350144/93‑1 (MQ=255)
                                                          gcCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATg                               >  2:6815823/1‑93 (MQ=255)
                                                          gcCTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATg                               >  1:4022424/1‑93 (MQ=255)
                                                           ccTCTTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGa                              >  1:9764409/1‑93 (MQ=255)
                                                             tctTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACa                            <  1:3065508/93‑1 (MQ=255)
                                                             tctTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACa                            <  1:6230235/93‑1 (MQ=255)
                                                             cctTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACa                            <  2:1378273/92‑1 (MQ=255)
                                                              ctTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACAt                           <  2:6576696/93‑1 (MQ=255)
                                                              ctTCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACAt                           <  2:5784800/93‑1 (MQ=255)
                                                               ttCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATg                          >  2:3341319/1‑93 (MQ=255)
                                                               ttCCAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATg                          >  2:1810469/1‑93 (MQ=255)
                                                                 ccAGGCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCg                        >  1:1414415/1‑93 (MQ=255)
                                                                     gCTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTaa                    >  1:8558106/1‑93 (MQ=255)
                                                                      cTCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAAc                   <  1:2307284/93‑1 (MQ=255)
                                                                       tCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACa                  >  1:6827666/1‑93 (MQ=255)
                                                                       tCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACa                  >  1:2078982/1‑93 (MQ=255)
                                                                       tCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACa                  >  1:8584459/1‑93 (MQ=255)
                                                                       tCATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACa                  >  2:2400558/1‑93 (MQ=255)
                                                                        cATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAg                 <  2:1006389/93‑1 (MQ=255)
                                                                        cATACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGCACTTCTCACCGGTAAAGTCAGCGCTATAACGCATGACATGCGGTAACAg                 >  1:2626447/1‑93 (MQ=255)
                                                                         aTACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAgg                <  2:3606840/93‑1 (MQ=255)
                                                                          tACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGa               >  2:2249863/1‑93 (MQ=255)
                                                                          tACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGa               >  2:8851914/1‑93 (MQ=255)
                                                                          tACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGa               <  1:4143109/93‑1 (MQ=255)
                                                                          tACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGa               >  1:7094018/1‑93 (MQ=255)
                                                                          tACCTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGAGGTAACAGGa               >  2:4459692/1‑93 (MQ=255)
                                                                           cccTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGAt              <  1:4425850/92‑1 (MQ=255)
                                                                            ccTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGCAACAGGATg             >  1:7363709/1‑93 (MQ=255)
                                                                             cTTCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATgg            <  1:1161733/93‑1 (MQ=255)
                                                                               tCCACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCg          >  2:5693636/1‑93 (MQ=255)
                                                                                ccACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGt         >  2:7482373/1‑93 (MQ=255)
                                                                                 cACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGtt        <  1:7284003/93‑1 (MQ=255)
                                                                                 cACTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACAGGCGGTAACAGGATGGCGtt        <  1:6009167/93‑1 (MQ=255)
                                                                                   cTTTCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCg      >  1:8985423/1‑93 (MQ=255)
                                                                                    tttCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCGc     >  2:9847506/1‑93 (MQ=255)
                                                                                      tCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCGCaa   <  2:6028202/93‑1 (MQ=255)
                                                                                      tCACACCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCGCaa   <  1:9967000/93‑1 (MQ=255)
                                                                                       cacaCCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCGCAac  >  2:9291199/1‑93 (MQ=255)
                                                                                       cacaCCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCGCAac  <  1:931392/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGGCGGAATACCGACATCACGGTTGAGAGCAAACACCGCTTCAACAGCGGCATTACGCGCCTCTTCCAGGCTCATACCTTCCACTTTCACGCCCATAACGCGCGCGATATCGCGGTACTTCTCACCGGTAAAGTCAGCGTTATAACGCATGACATGCGGTAACAGGATGGCGTTCGCAAC  >  NC_000913/2930025‑2930204

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: