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.4% 478.0 / ‑4.6 154H106R (CAC→CGC) wzbputative protein‑tyrosine‑phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (73/80);  total (74/80)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 158 total aligned reads displayed.

GGGCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCG  >  NC_000913/2133715‑2133893
                                                                                         |                                                                                         
gggCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGAc                                                                                        >  1:2742078/1‑93 (MQ=255)
 ggCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGAcc                                                                                       >  1:1400093/1‑93 (MQ=255)
 ggCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGAcc                                                                                       >  2:2660520/1‑93 (MQ=255)
    aGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCaaa                                                                                    >  2:4575305/1‑93 (MQ=255)
     gACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGGTTTCACATTCGTTATCCCAGCGAGCAAAc                                                                                   >  2:4102044/1‑93 (MQ=255)
     gACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAAc                                                                                   <  1:1807031/93‑1 (MQ=255)
      aCCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACa                                                                                  <  1:2865591/93‑1 (MQ=255)
      aCCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACa                                                                                  <  2:2247967/93‑1 (MQ=255)
        cGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGc                                                                                >  2:2123166/1‑93 (MQ=255)
          ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              <  2:676770/93‑1 (MQ=255)
          ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              <  1:4277932/93‑1 (MQ=255)
          ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              >  2:516359/1‑93 (MQ=255)
             aaGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCAc                                                                           >  2:2939608/1‑93 (MQ=255)
             aaGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCAc                                                                           >  2:1937290/1‑93 (MQ=255)
             aaGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCAc                                                                           <  2:1007127/93‑1 (MQ=255)
              aGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACt                                                                          <  2:2495254/93‑1 (MQ=255)
               gTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACtt                                                                         <  2:1353189/93‑1 (MQ=255)
               atAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACtt                                                                         <  1:3525048/92‑1 (MQ=255)
                tAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACttt                                                                        >  2:2606842/1‑93 (MQ=255)
                tAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACttt                                                                        >  2:66152/1‑93 (MQ=255)
                tAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACttt                                                                        >  2:934462/1‑93 (MQ=255)
                tAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACttt                                                                        <  1:1166359/93‑1 (MQ=255)
                tAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACttt                                                                        <  1:1481327/93‑1 (MQ=255)
                   tgtgTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGcc                                                                     >  1:1305946/1‑93 (MQ=255)
                    gtgtACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCcg                                                                    >  1:1636232/1‑93 (MQ=255)
                     tgtACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCcgc                                                                   <  1:1721148/93‑1 (MQ=255)
                       tACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCgcgc                                                                 <  2:4050128/93‑1 (MQ=255)
                        acacCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCa                                                                <  2:3778918/93‑1 (MQ=255)
                          cccGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtc                                                              <  1:2408806/92‑1 (MQ=255)
                          acCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtc                                                              <  1:1076160/93‑1 (MQ=255)
                          acCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtc                                                              <  2:3987477/93‑1 (MQ=255)
                           ccGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtct                                                             <  2:318594/93‑1 (MQ=255)
                             gCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTcc                                                           <  2:1691943/93‑1 (MQ=255)
                             gCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTcc                                                           <  2:179155/93‑1 (MQ=255)
                             gCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTcc                                                           <  1:769765/93‑1 (MQ=255)
                             gCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTcc                                                           >  1:594632/1‑93 (MQ=255)
                               tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         <  1:4076228/93‑1 (MQ=255)
                               tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         >  2:3450874/1‑93 (MQ=255)
                               tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         >  1:502997/1‑93 (MQ=255)
                                gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        <  2:2554475/93‑1 (MQ=255)
                                gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        >  2:3441490/1‑93 (MQ=255)
                                gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        <  2:53257/93‑1 (MQ=255)
                                   aaCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGcc                                                     <  1:2029235/93‑1 (MQ=255)
                                    aCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGGCa                                                    >  1:3294903/1‑93 (MQ=255)
                                    aCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCa                                                    <  2:552687/93‑1 (MQ=255)
                                      gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATTCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtc                                                  >  1:3734671/1‑93 (MQ=255)
                                      gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtc                                                  <  2:3161698/93‑1 (MQ=255)
                                      gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtc                                                  >  2:2891579/1‑93 (MQ=255)
                                       tttGCCGGGTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtct                                                 <  1:375577/93‑1 (MQ=255)
                                       tttCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtct                                                 >  2:3302334/1‑93 (MQ=255)
                                           ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCa                                             <  1:2237649/93‑1 (MQ=255)
                                           ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGACATCTCGCa                                             >  1:4207567/1‑93 (MQ=255)
                                           ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCCGTGCCATCTCGCa                                             >  2:2885797/1‑93 (MQ=255)
                                             ggCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATa                                           >  2:1735409/1‑93 (MQ=255)
                                              gCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATaa                                          >  2:1739264/1‑93 (MQ=255)
                                               cTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAg                                         <  2:38849/93‑1 (MQ=255)
                                               cTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAg                                         <  1:793021/93‑1 (MQ=255)
                                                ttttGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGc                                        >  2:2092783/1‑93 (MQ=255)
                                                   tGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGtt                                     <  2:3211466/93‑1 (MQ=255)
                                                   tGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGtt                                     <  2:2027324/93‑1 (MQ=255)
                                                   tGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGtt                                     <  1:35165/93‑1 (MQ=255)
                                                    gCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTc                                    >  1:124671/1‑93 (MQ=255)
                                                      gATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGa                                  <  1:111535/93‑1 (MQ=255)
                                                        tACTGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  2:1844425/93‑1 (MQ=255)
                                                        tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  1:3899527/93‑1 (MQ=255)
                                                         aCGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGatat                               >  1:611576/1‑93 (MQ=255)
                                                          cGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATg                              <  1:200348/93‑1 (MQ=255)
                                                          cGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGATGACATCTCGCATAAGCGTTCGATATg                              >  1:3058939/1‑93 (MQ=255)
                                                           ggATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATgg                             >  1:2942518/1‑93 (MQ=255)
                                                            gATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgc                            <  1:210760/93‑1 (MQ=255)
                                                             aTCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgcg                           >  2:4266655/1‑93 (MQ=255)
                                                               cGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCTTCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCt                         <  1:1870007/93‑1 (MQ=255)
                                                               cGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCt                         >  1:4686432/1‑93 (MQ=255)
                                                                ggggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCtt                        >  1:87796/1‑93 (MQ=255)
                                                                ggggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCtt                        >  1:2709805/1‑93 (MQ=255)
                                                                 gggATTTCACATTGTTTATCCGAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCttt                       <  2:2699687/93‑1 (MQ=255)
                                                                 gggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCttt                       >  2:3891587/1‑93 (MQ=255)
                                                                 gggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCttt                       >  1:2072899/1‑93 (MQ=255)
                                                                 gggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCttt                       >  2:1675667/1‑93 (MQ=255)
                                                                  ggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCtttt                      <  1:3254696/93‑1 (MQ=255)
                                                                   gATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTc                     >  1:3414685/1‑93 (MQ=255)
                                                                     tttCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCa                   <  1:3676604/93‑1 (MQ=255)
                                                                      ttCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCAt                  <  2:4441234/93‑1 (MQ=255)
                                                                        cacaTTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATgg                <  2:3310800/93‑1 (MQ=255)
                                                                          caTTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTc              >  2:3925114/1‑93 (MQ=255)
                                                                          caTTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTc              >  2:1984332/1‑93 (MQ=255)
                                                                           aTTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCa             >  1:4165657/1‑93 (MQ=255)
                                                                            ttCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCaa            >  2:1128684/1‑93 (MQ=255)
                                                                            ttCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCaa            <  1:1985844/93‑1 (MQ=255)
                                                                              cGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCaaaa          <  2:3579970/93‑1 (MQ=255)
                                                                              cGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCaaaa          >  1:2292544/1‑93 (MQ=255)
                                                                                ttATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATc        <  1:696156/93‑1 (MQ=255)
                                                                                ttATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATc        >  2:1481941/1‑93 (MQ=255)
                                                                                  aTCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAg      <  1:1952184/93‑1 (MQ=255)
                                                                                   tCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAgg     >  1:3416961/1‑93 (MQ=255)
                                                                                   tCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAgg     >  2:4323813/1‑93 (MQ=255)
                                                                                    cccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGt    <  2:3060693/93‑1 (MQ=255)
                                                                                     ccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTc   >  2:2564344/1‑93 (MQ=255)
                                                                                     ccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTc   >  1:241932/1‑93 (MQ=255)
                                                                                      cAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCg  <  2:2486075/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 < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: