Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,133,804 T→C H106R (CAC→CGC)  wzb ← putative protein‑tyrosine‑phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,133,8040TC99.6% 781.7 / ‑4.9 246H106R (CAC→CGC) wzbputative protein‑tyrosine‑phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (109/136);  total (110/136)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

Only 100 of 260 total aligned reads displayed.

GGGCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCG  >  NC_000913/2133715‑2133893
                                                                                         |                                                                                         
gggCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGGTATCCCAGCGAc                                                                                        >  1:1043577/1‑93 (MQ=255)
 ggCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGAcc                                                                                       >  1:4589798/1‑93 (MQ=255)
 ggCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGAcc                                                                                       >  1:3455996/1‑93 (MQ=255)
  gCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCa                                                                                      >  2:1217659/1‑93 (MQ=255)
   cAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCaa                                                                                     >  1:4153306/1‑93 (MQ=255)
    aGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCcaa                                                                                    >  1:7076110/1‑93 (MQ=255)
    aGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCaaa                                                                                    <  2:2639727/93‑1 (MQ=255)
     gACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAAc                                                                                   >  1:6989422/1‑93 (MQ=255)
     gACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAAc                                                                                   <  2:7961460/93‑1 (MQ=255)
     gACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCAGGGATTTCACATTCGTTATCCCAGCGACCAAAc                                                                                   <  1:678651/93‑1 (MQ=255)
       ccGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACac                                                                                 >  2:9766539/1‑92 (MQ=255)
       ccGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAg                                                                                 >  1:4957254/1‑93 (MQ=255)
        cGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGc                                                                                >  2:5993631/1‑93 (MQ=255)
        cGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGc                                                                                <  1:3634809/93‑1 (MQ=255)
         gTTCAAGTGATGTGTACGCCGCTGCAAACGTTTCCCGGCTTATACGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCa                                                                               <  2:8979890/93‑1 (MQ=255)
          ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              <  1:1880625/93‑1 (MQ=255)
          ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              <  2:7657741/93‑1 (MQ=255)
          ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTCATCCCAGCGACCAAACAGCAt                                                                              >  2:2801680/1‑93 (MQ=255)
           tCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATc                                                                             >  2:3638675/1‑93 (MQ=255)
            cAAGAAATGTGTCCACCGCGGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCTTTTTCCCAGCGACCAACAGCATCAc                                                                           <  2:4542118/93‑1 (MQ=255)
              aGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACt                                                                          >  2:6575947/1‑93 (MQ=255)
                tAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACttt                                                                        >  1:8269876/1‑93 (MQ=255)
                tAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACttt                                                                        <  2:4370664/93‑1 (MQ=255)
                  aTGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGc                                                                      <  2:354275/93‑1 (MQ=255)
                   tgtgTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGcc                                                                     >  2:254052/1‑93 (MQ=255)
                    gtgtACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCcg                                                                    >  2:1865675/1‑93 (MQ=255)
                    gtgtACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGACTTGACATTCGCTATCCCAGCGACCAAACAGCATAACTTTGctg                                                                    >  2:4960544/1‑91 (MQ=255)
                        acacCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCa                                                                <  2:6739492/93‑1 (MQ=255)
                        acacCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTAGCCGCGCa                                                                <  2:324888/93‑1 (MQ=255)
                           ccGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtct                                                             >  2:5252720/1‑93 (MQ=255)
                           ccGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtct                                                             <  1:1237050/93‑1 (MQ=255)
                           ccGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGGTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtct                                                             <  1:3282870/93‑1 (MQ=255)
                            cGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtctc                                                            <  1:2606356/93‑1 (MQ=255)
                            cGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtctc                                                            <  2:8080275/93‑1 (MQ=255)
                              cTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCg                                                          <  2:1439219/93‑1 (MQ=255)
                              cTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCg                                                          >  1:9079793/1‑93 (MQ=255)
                              cTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCg                                                          >  1:9993467/1‑93 (MQ=255)
                              cTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCg                                                          <  1:5957647/93‑1 (MQ=255)
                              cTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCg                                                          <  1:6616361/93‑1 (MQ=255)
                               tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         <  2:7406806/93‑1 (MQ=255)
                               tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         <  1:7739341/93‑1 (MQ=255)
                               tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         <  1:6262028/93‑1 (MQ=255)
                               tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         <  1:8820196/93‑1 (MQ=255)
                                 cAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTg                                                       >  1:3234965/1‑93 (MQ=255)
                                 cAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTg                                                       >  1:3146159/1‑93 (MQ=255)
                                 cAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTg                                                       >  2:4037002/1‑93 (MQ=255)
                                 cAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTg                                                       >  2:789076/1‑93 (MQ=255)
                                   aaCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGcc                                                     >  2:3332373/1‑93 (MQ=255)
                                       tttCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtct                                                 >  1:8701319/1‑93 (MQ=255)
                                        ttCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtctc                                                <  2:8371290/93‑1 (MQ=255)
                                           ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCa                                             <  1:2482226/93‑1 (MQ=255)
                                           ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCa                                             >  2:4822262/1‑93 (MQ=255)
                                             ggCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATa                                           >  1:4435520/1‑93 (MQ=255)
                                              gCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATaa                                          >  2:2451539/1‑93 (MQ=255)
                                                ttttGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGc                                        >  2:7258696/1‑93 (MQ=255)
                                                ttttGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGc                                        >  1:9789997/1‑93 (MQ=255)
                                                   tGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGtt                                     <  1:4801752/93‑1 (MQ=255)
                                                   tGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGtt                                     <  1:9016951/93‑1 (MQ=255)
                                                    gCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTc                                    >  2:9602862/1‑93 (MQ=255)
                                                    gCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTc                                    <  1:9705831/93‑1 (MQ=255)
                                                     cGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCCGTGCCATCTCGCATAAGCGTTCg                                   >  2:5524005/1‑93 (MQ=255)
                                                      gATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGa                                  >  1:1962445/1‑93 (MQ=255)
                                                      gATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGa                                  <  2:6443069/93‑1 (MQ=255)
                                                       aTACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGat                                 >  1:8305097/1‑93 (MQ=255)
                                                        tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  2:2764840/93‑1 (MQ=255)
                                                        tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  1:4043704/93‑1 (MQ=255)
                                                        tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  1:5740487/93‑1 (MQ=255)
                                                        tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  1:7099856/93‑1 (MQ=255)
                                                        tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  2:8809774/93‑1 (MQ=255)
                                                        tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  2:1785082/93‑1 (MQ=255)
                                                         aCGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGatat                               >  2:4160702/1‑93 (MQ=255)
                                                          cGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATg                              <  1:4743981/93‑1 (MQ=255)
                                                           ggATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATgg                             <  1:4830580/93‑1 (MQ=255)
                                                            gATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgc                            >  2:3521649/1‑93 (MQ=255)
                                                            gATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgc                            <  1:7943963/93‑1 (MQ=255)
                                                             aTCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgcg                           >  1:6305420/1‑93 (MQ=255)
                                                             aTCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgcg                           >  2:5070659/1‑93 (MQ=255)
                                                             aTCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgcg                           <  1:3957801/93‑1 (MQ=255)
                                                              tCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgcgc                          <  1:5464646/93‑1 (MQ=255)
                                                               cGGTGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCt                         <  2:3069387/93‑1 (MQ=255)
                                                                ggggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCtt                        >  1:2981216/1‑93 (MQ=255)
                                                                 gggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCttt                       <  2:6400636/93‑1 (MQ=255)
                                                                 gggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGCCGCttt                       >  2:8572987/1‑93 (MQ=255)
                                                                    aTTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTcc                    >  2:913198/1‑93 (MQ=255)
                                                                     tttCACATTCGTTATCCCAGCGACCAAACAGCATCCCTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGAGTTCGATATGGCGCTTTTCCa                   <  1:1808467/93‑1 (MQ=255)
                                                                     tttCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCa                   <  1:6289617/93‑1 (MQ=255)
                                                                     tttCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCa                   <  1:9461945/93‑1 (MQ=255)
                                                                      ttCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCAt                  >  2:1067138/1‑93 (MQ=255)
                                                                       tCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATg                 <  1:9158458/93‑1 (MQ=255)
                                                                           aTTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCa             >  2:3917340/1‑93 (MQ=255)
                                                                              cGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCaaaa          >  1:72638/1‑93 (MQ=255)
                                                                              cGTTATCCCAGCGACCAAACAGAATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCaaaa          <  1:7911563/93‑1 (MQ=255)
                                                                                ttATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATc        >  2:20661/1‑93 (MQ=255)
                                                                                 tATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCa       >  1:1532363/1‑93 (MQ=255)
                                                                                 tATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCAGGGTCAAAATCa       <  1:284729/93‑1 (MQ=255)
                                                                                   tCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGGGCCATCTCGCATAAGCGTTCGATATGGCGCTTCTCCATGGTCAAAATCAgg     >  1:5353666/1‑93 (MQ=255)
                                                                                    cccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGt    <  1:579708/93‑1 (MQ=255)
                                                                                    cccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGt    <  2:4657849/93‑1 (MQ=255)
                                                                                    cccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGt    >  1:1425964/1‑93 (MQ=255)
                                                                                      cAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCg  <  2:563796/93‑1 (MQ=255)
                                                                                         |                                                                                         
GGGCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCG  >  NC_000913/2133715‑2133893

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