Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,121,597 (T)5→4 coding (151/1332 nt) hslU ← molecular chaperone and ATPase component of HslUV protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,121,5930T.96.5% 624.7 / 13.5 145coding (155/1332 nt)hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base . (71/69);  total (72/73)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 159 total aligned reads displayed.

TGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCC  >  NC_000913/4121495‑4121688
                                                                                                  |                                                                                               
  gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtttt                                                                                              >  1:135585/1‑100 (MQ=255)
   tCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttc                                                                                             >  2:2590855/1‑99 (MQ=255)
     ncTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           >  2:1985775/2‑97 (MQ=255)
     ncTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           >  2:636903/2‑97 (MQ=255)
     gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           >  1:2467326/1‑97 (MQ=255)
     gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           <  2:2012707/100‑4 (MQ=255)
     gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           <  2:3192982/100‑4 (MQ=255)
      cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtttttggg                                                                                          >  2:1361208/1‑97 (MQ=255)
      cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  2:897149/100‑5 (MQ=255)
      cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           >  2:1051035/1‑96 (MQ=255)
       ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCgggg                                                                                        >  1:1720208/1‑100 (MQ=255)
        tCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGt                                                                                       >  2:2393253/1‑100 (MQ=255)
         cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTc                                                                                      <  2:1346104/100‑1 (MQ=255)
         cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTc                                                                                      >  2:863624/1‑100 (MQ=255)
         cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTGTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  2:2863508/97‑5 (MQ=255)
          aaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCa                                                                                     >  2:1278348/1‑100 (MQ=255)
            cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACt                                                                                   >  2:2722013/1‑100 (MQ=255)
             tttGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACtt                                                                                  <  1:1873184/100‑1 (MQ=255)
                 atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATg                                                                              <  1:686241/100‑1 (MQ=255)
                 atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATg                                                                              >  2:1727719/1‑100 (MQ=255)
                 atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATg                                                                              >  2:2904063/1‑100 (MQ=255)
                     aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCCATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCGCTTGATGgcgc                                                                          <  2:1765656/100‑1 (MQ=255)
                      cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCa                                                                         >  1:744715/1‑100 (MQ=255)
                        ncgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGc                                                                       >  2:747378/2‑100 (MQ=255)
                        gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGc                                                                       <  2:2656967/100‑1 (MQ=255)
                        gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGc                                                                       >  1:259874/1‑100 (MQ=255)
                         cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGct                                                                      >  1:2280097/1‑100 (MQ=255)
                         cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGct                                                                      <  1:2970738/100‑1 (MQ=255)
                            aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   >  1:3112289/1‑100 (MQ=255)
                             ttCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTc                                                                  >  2:673546/1‑100 (MQ=255)
                               cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                                >  2:633458/1‑100 (MQ=255)
                                gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                               >  1:604677/1‑100 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              >  2:1227103/1‑100 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:2818884/100‑1 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:1218171/100‑1 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:1992464/100‑1 (MQ=255)
                                  cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  1:2935488/100‑1 (MQ=255)
                                  cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  1:369930/100‑1 (MQ=255)
                                    gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                           <  2:2306037/100‑1 (MQ=255)
                                    gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                           <  2:201936/100‑1 (MQ=255)
                                      ttAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTg                                                         <  2:1080877/100‑1 (MQ=255)
                                        aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  1:1840983/100‑1 (MQ=255)
                                         gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  2:2802761/1‑100 (MQ=255)
                                         gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  1:3146100/1‑100 (MQ=255)
                                         gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  2:1117885/1‑100 (MQ=255)
                                            aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  2:1781365/1‑100 (MQ=255)
                                            aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  2:2339225/1‑100 (MQ=255)
                                             gacgacGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGa                                                  <  2:555768/100‑1 (MQ=255)
                                               cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACg                                                <  1:1841497/100‑1 (MQ=255)
                                                 acgGGCGATTTCAGTTTTAACGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCACGGTGCAGCTCAACGAAGAGCTGcgccatgaagt                                              <  1:1108443/100‑11 (MQ=21)
                                                    ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                           >  1:2578119/1‑100 (MQ=255)
                                                     gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCg                                          <  1:1727724/100‑1 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         >  2:376566/1‑100 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         <  1:1117886/100‑1 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         <  1:2641977/100‑1 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         <  2:1720203/100‑1 (MQ=255)
                                                       gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        <  2:2565969/100‑1 (MQ=255)
                                                        aTTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGtt                                       <  1:2813236/100‑1 (MQ=255)
                                                         tttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTa                                      >  2:1122350/1‑100 (MQ=255)
                                                          ttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTa                                      <  2:2208933/100‑1 (MQ=255)
                                                          ttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTAc                                     >  2:314954/1‑100 (MQ=255)
                                                             ccgatctccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                  <  1:648260/93‑1 (MQ=255)
                                                             aGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                  <  1:2393264/100‑1 (MQ=255)
                                                              gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  2:2952934/1‑100 (MQ=255)
                                                              gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 <  1:1060487/100‑1 (MQ=255)
                                                              gTTTTACCGACACCCGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 <  2:2280086/100‑1 (MQ=255)
                                                                tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               >  2:2777026/1‑100 (MQ=255)
                                                                 tctccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                              <  2:3277669/97‑1 (MQ=255)
                                                                  tACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCa                             >  2:3235041/1‑100 (MQ=255)
                                                                    ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCaagatcgg                           >  2:648259/1‑93 (MQ=255)
                                                                    ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCaga                           >  1:3277687/1‑98 (MQ=255)
                                                                    ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1837922/1‑100 (MQ=255)
                                                                        cacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       >  2:1925765/1‑100 (MQ=255)
                                                                        cacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       <  1:2165586/100‑1 (MQ=255)
                                                                        cacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       <  1:625894/100‑1 (MQ=255)
                                                                          ccGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                     >  2:2572747/1‑100 (MQ=255)
                                                                            ggTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                   <  2:207634/100‑1 (MQ=255)
                                                                            ggTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                   >  2:975084/1‑100 (MQ=255)
                                                                             gTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  >  2:46060/1‑100 (MQ=255)
                                                                              tCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGaa                 >  1:2858831/1‑100 (MQ=255)
                                                                               cGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAc                <  2:1276839/100‑1 (MQ=255)
                                                                                gggCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACg               >  1:1233476/1‑100 (MQ=255)
                                                                                   ccGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCtt            <  1:1399403/100‑1 (MQ=255)
                                                                                    cGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  1:789479/1‑100 (MQ=255)
                                                                                    cGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  1:1961844/100‑1 (MQ=255)
                                                                                    cGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  2:2367068/1‑100 (MQ=255)
                                                                                     gATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTgg          <  1:1053283/100‑1 (MQ=255)
                                                                                     gATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTgg          <  1:3079811/100‑1 (MQ=255)
                                                                                      atcatcAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc         >  2:2029516/1‑100 (MQ=255)
                                                                                      atcatcAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc         <  1:863626/100‑1 (MQ=255)
                                                                                       tcatcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  2:2722599/100‑1 (MQ=255)
                                                                                       tcatcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  1:2722026/100‑1 (MQ=255)
                                                                                       tcatcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  1:3096391/100‑1 (MQ=255)
                                                                                       tcatcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  2:2601805/100‑1 (MQ=255)
                                                                                          tcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:3202062/100‑1 (MQ=255)
                                                                                           caGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGt    >  1:1439963/1‑100 (MQ=255)
                                                                                            aGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc   >  2:1565547/1‑100 (MQ=255)
                                                                                             ggATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  1:2563201/1‑100 (MQ=255)
                                                                                             ggATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  1:1805554/1‑100 (MQ=255)
                                                                                             ggATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  <  1:314955/100‑1 (MQ=255)
                                                                                                  |                                                                                               
TGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCC  >  NC_000913/4121495‑4121688

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