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% 551.9 / ‑5.1 173H106R (CAC→CGC) wzbputative protein‑tyrosine‑phosphatase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (90/82);  total (91/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.08e-01

Only 100 of 177 total aligned reads displayed.

CGGGCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCG  >  NC_000913/2133714‑2133893
                                                                                          |                                                                                         
cGGGCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGa                                                                                         >  1:2544075/1‑93 (MQ=255)
 gggCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGAc                                                                                        >  2:4669843/1‑93 (MQ=255)
   gCAGACCGTTCAAGCAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCCGCGATTTCCAATTCCTTATCCCCGAGACCa                                                                                      >  2:370566/1‑93 (MQ=255)
    cAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCaa                                                                                     <  1:1040396/93‑1 (MQ=255)
    cAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCaa                                                                                     >  1:2099906/1‑93 (MQ=255)
       aCCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACa                                                                                  >  2:4484932/1‑93 (MQ=255)
       aCCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACa                                                                                  <  2:3422481/93‑1 (MQ=255)
        ccGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAg                                                                                 >  2:95350/1‑93 (MQ=255)
        ccGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAg                                                                                 <  2:2270456/93‑1 (MQ=255)
        ccGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAg                                                                                 <  1:3552303/93‑1 (MQ=255)
         cGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGc                                                                                <  1:2508670/93‑1 (MQ=255)
         cGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGc                                                                                >  1:3627005/1‑93 (MQ=255)
          gTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCa                                                                               <  1:3902146/93‑1 (MQ=255)
          gTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCa                                                                               >  1:427496/1‑93 (MQ=255)
           ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCTGCGACCAAACAGCAt                                                                              >  1:1907155/1‑93 (MQ=255)
           ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              <  1:698436/93‑1 (MQ=255)
           ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              >  1:2248163/1‑93 (MQ=255)
           ttCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCAt                                                                              <  2:2502315/93‑1 (MQ=255)
              aaGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCCAACAGGATCAc                                                                           >  1:1643246/1‑93 (MQ=255)
              aaGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCAc                                                                           >  2:1351768/1‑93 (MQ=255)
              aaGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCAc                                                                           >  1:3069734/1‑93 (MQ=255)
                gTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACtt                                                                         <  1:2371333/93‑1 (MQ=255)
                    tgtgTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGcc                                                                     <  1:513244/93‑1 (MQ=255)
                       gtACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCcgcg                                                                  >  2:239324/1‑93 (MQ=255)
                        tACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCgcgc                                                                 >  1:109974/1‑93 (MQ=255)
                        tACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCgcgc                                                                 <  2:1301170/93‑1 (MQ=255)
                          cacCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAt                                                               >  1:11886/1‑93 (MQ=255)
                           acCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtc                                                              >  2:4279974/1‑93 (MQ=255)
                             cGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCAtctc                                                            <  1:3677614/93‑1 (MQ=255)
                              gCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTcc                                                           <  1:1643168/93‑1 (MQ=255)
                              gCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTcc                                                           <  2:2973242/93‑1 (MQ=255)
                                tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         >  1:4656603/1‑93 (MQ=255)
                                tGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCgg                                                         >  1:872995/1‑93 (MQ=255)
                                 gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        <  2:1982653/93‑1 (MQ=255)
                                 gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        >  2:1543711/1‑93 (MQ=255)
                                 gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        >  1:1903189/1‑93 (MQ=255)
                                 gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        <  2:3320923/93‑1 (MQ=255)
                                 gCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        <  2:2939330/93‑1 (MQ=255)
                                 gCAAACGTTGCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGt                                                        <  2:1982654/93‑1 (MQ=255)
                                   aaaCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGc                                                      <  1:4099221/93‑1 (MQ=255)
                                    aaCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGcc                                                     >  1:665303/1‑93 (MQ=255)
                                    aaCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACAAAACAGCATCACTTTGCCGCGCATCTCCGGTGcc                                                     >  1:870227/1‑93 (MQ=255)
                                     aCGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCa                                                    >  2:4348623/1‑93 (MQ=255)
                                       gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtc                                                  >  2:1451185/1‑93 (MQ=255)
                                       gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtc                                                  <  1:2347843/93‑1 (MQ=255)
                                       gTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtc                                                  >  1:3510095/1‑93 (MQ=255)
                                        tttCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtct                                                 >  2:1877966/1‑93 (MQ=255)
                                        tttCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtct                                                 <  1:4649910/93‑1 (MQ=255)
                                        tttCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtct                                                 <  1:3467894/93‑1 (MQ=255)
                                         ttCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCAtctc                                                <  2:1079061/93‑1 (MQ=255)
                                          tCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCg                                               >  2:3582416/1‑93 (MQ=255)
                                            ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCa                                             <  1:3893095/93‑1 (MQ=255)
                                            ccGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCa                                             <  1:4579798/93‑1 (MQ=255)
                                             cGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCAt                                            >  1:1391128/1‑93 (MQ=255)
                                             cGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCAt                                            >  1:1035994/1‑93 (MQ=255)
                                              ggCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATa                                           <  2:4273777/93‑1 (MQ=255)
                                              ggCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATa                                           <  1:4185001/93‑1 (MQ=255)
                                               gCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATaa                                          <  2:1903928/93‑1 (MQ=255)
                                               gCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATaa                                          >  1:1451899/1‑93 (MQ=255)
                                                 ttttTCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGc                                        <  1:2329160/93‑1 (MQ=255)
                                                 ttttGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGc                                        >  2:4652417/1‑93 (MQ=255)
                                                 ttttGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGc                                        >  1:3299644/1‑93 (MQ=255)
                                                    tGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGtt                                     <  1:1698569/93‑1 (MQ=255)
                                                       gATACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGa                                  >  1:620776/1‑93 (MQ=255)
                                                        aTACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGat                                 <  1:1798622/93‑1 (MQ=255)
                                                         tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                >  1:676987/1‑93 (MQ=255)
                                                         tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                >  2:615870/1‑93 (MQ=255)
                                                         tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                <  2:3411612/93‑1 (MQ=255)
                                                         tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                >  1:602166/1‑93 (MQ=255)
                                                         tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGata                                >  2:1389833/1‑93 (MQ=255)
                                                         tACGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAAAAGCATCAGTTTGCCGGGCATCTCCGGTGCCATCTCGCATAAGCGTTAGata                                >  1:853391/1‑93 (MQ=255)
                                                          aCGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGatat                               >  1:541603/1‑93 (MQ=255)
                                                          aCGGATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGatat                               <  1:3208945/93‑1 (MQ=255)
                                                           cGGATTGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATg                              <  1:793282/93‑1 (MQ=255)
                                                            ggATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGCTATgg                             >  1:618374/1‑93 (MQ=255)
                                                            ggATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATgg                             <  2:2552195/93‑1 (MQ=255)
                                                            ggATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATgg                             <  1:90080/93‑1 (MQ=255)
                                                            ggATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATgg                             <  1:2519594/93‑1 (MQ=255)
                                                             gATCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgc                            <  1:3820822/93‑1 (MQ=255)
                                                               tCGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgcgc                          <  1:1247687/93‑1 (MQ=255)
                                                               tCGGGGATTTCACATTCGTTATCCCAGCAACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGgcgc                          <  1:4231579/93‑1 (MQ=255)
                                                                cGGGGATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCt                         <  1:2897045/93‑1 (MQ=255)
                                                                  gggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCttt                       >  2:3933167/1‑93 (MQ=255)
                                                                   ggATTTCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCtttt                      >  1:1728594/1‑93 (MQ=255)
                                                                      tttCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCa                   >  1:3129955/1‑93 (MQ=255)
                                                                      tttCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCa                   >  2:1745448/1‑93 (MQ=255)
                                                                       ttCACATTCGTTATCCCAGCGACCAAACAGCATCACTATGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCAt                  <  1:3229396/93‑1 (MQ=255)
                                                                        tCACATTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATg                 >  1:1278004/1‑93 (MQ=255)
                                                                         cacaTTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATgg                >  2:3037349/1‑93 (MQ=255)
                                                                           caTTCGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTc              >  2:2526638/1‑93 (MQ=255)
                                                                               cGTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCaaaa          <  1:2942754/93‑1 (MQ=255)
                                                                                gTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAAt         <  2:264722/93‑1 (MQ=255)
                                                                                gTTATCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAAt         <  1:3382456/93‑1 (MQ=255)
                                                                                   aTCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAg      >  1:2868349/1‑93 (MQ=255)
                                                                                   aTCCCAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAg      <  1:1987033/93‑1 (MQ=255)
                                                                                     cccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGt    >  2:2042878/1‑93 (MQ=255)
                                                                                      ccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTc   <  1:2792658/93‑1 (MQ=255)
                                                                                      ccAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTc   >  2:4634595/1‑93 (MQ=255)
                                                                                       cAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCg  >  1:4495138/1‑93 (MQ=255)
                                                                                       cAGCGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCg  >  2:1735430/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGGGCAGACCGTTCAAGTAATGTGTACACCGCTGCAAACGTTTCCCGGCTTTTGCGATACGGATCGGGGATTTCACATTCGTTATCCCAGTGACCAAACAGCATCACTTTGCCGCGCATCTCCGGTGCCATCTCGCATAAGCGTTCGATATGGCGCTTTTCCATGGTCAAAATCAGGTCG  >  NC_000913/2133714‑2133893

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