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)
       ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCgggg                                                                                        >  1:1720208/1‑100 (MQ=255)
        tCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGt                                                                                       >  2:2393253/1‑100 (MQ=255)
         cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTc                                                                                      <  2:1346104/100‑1 (MQ=255)
         cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTc                                                                                      >  2:863624/1‑100 (MQ=255)
         cAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTGTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  2:2863508/97‑5 (MQ=255)
          aaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCa                                                                                     >  2:1278348/1‑100 (MQ=255)
            cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACt                                                                                   >  2:2722013/1‑100 (MQ=255)
             tttGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACtt                                                                                  <  1:1873184/100‑1 (MQ=255)
                 atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATg                                                                              <  1:686241/100‑1 (MQ=255)
                 atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATg                                                                              >  2:1727719/1‑100 (MQ=255)
                 atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATg                                                                              >  2:2904063/1‑100 (MQ=255)
                     aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCCATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCGCTTGATGgcgc                                                                          <  2:1765656/100‑1 (MQ=255)
                      cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCa                                                                         >  1:744715/1‑100 (MQ=255)
                        ncgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGc                                                                       >  2:747378/2‑100 (MQ=255)
                        gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGc                                                                       <  2:2656967/100‑1 (MQ=255)
                        gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGc                                                                       >  1:259874/1‑100 (MQ=255)
                         cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGct                                                                      >  1:2280097/1‑100 (MQ=255)
                         cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGct                                                                      <  1:2970738/100‑1 (MQ=255)
                            aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   >  1:3112289/1‑100 (MQ=255)
                             ttCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTc                                                                  >  2:673546/1‑100 (MQ=255)
                               cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                                >  2:633458/1‑100 (MQ=255)
                                gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                               >  1:604677/1‑100 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              >  2:1227103/1‑100 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:2818884/100‑1 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:1218171/100‑1 (MQ=255)
                                 ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:1992464/100‑1 (MQ=255)
                                  cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  1:2935488/100‑1 (MQ=255)
                                  cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  1:369930/100‑1 (MQ=255)
                                    gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                           <  2:2306037/100‑1 (MQ=255)
                                    gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                           <  2:201936/100‑1 (MQ=255)
                                      ttAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTg                                                         <  2:1080877/100‑1 (MQ=255)
                                        aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  1:1840983/100‑1 (MQ=255)
                                         gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  2:2802761/1‑100 (MQ=255)
                                         gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  1:3146100/1‑100 (MQ=255)
                                         gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  2:1117885/1‑100 (MQ=255)
                                            aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  2:1781365/1‑100 (MQ=255)
                                            aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  2:2339225/1‑100 (MQ=255)
                                             gacgacGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGa                                                  <  2:555768/100‑1 (MQ=255)
                                               cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACg                                                <  1:1841497/100‑1 (MQ=255)
                                                 acgGGCGATTTCAGTTTTAACGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCACGGTGCAGCTCAACGAAGAGCTGcgccatgaagt                                              <  1:1108443/100‑11 (MQ=21)
                                                    ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                           >  1:2578119/1‑100 (MQ=255)
                                                     gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCg                                          <  1:1727724/100‑1 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         >  2:376566/1‑100 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         <  1:1117886/100‑1 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         <  1:2641977/100‑1 (MQ=255)
                                                      cGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCgg                                         <  2:1720203/100‑1 (MQ=255)
                                                       gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        <  2:2565969/100‑1 (MQ=255)
                                                        aTTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGtt                                       <  1:2813236/100‑1 (MQ=255)
                                                         tttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTa                                      >  2:1122350/1‑100 (MQ=255)
                                                          ttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTa                                      <  2:2208933/100‑1 (MQ=255)
                                                          ttCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTAc                                     >  2:314954/1‑100 (MQ=255)
                                                             ccgatctccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                  <  1:648260/93‑1 (MQ=255)
                                                             aGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                  <  1:2393264/100‑1 (MQ=255)
                                                              gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  2:2952934/1‑100 (MQ=255)
                                                              gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 <  1:1060487/100‑1 (MQ=255)
                                                              gTTTTACCGACACCCGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 <  2:2280086/100‑1 (MQ=255)
                                                                tttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               >  2:2777026/1‑100 (MQ=255)
                                                                 tctccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                              <  2:3277669/97‑1 (MQ=255)
                                                                  tACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCa                             >  2:3235041/1‑100 (MQ=255)
                                                                    ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCaagatcgg                           >  2:648259/1‑93 (MQ=255)
                                                                    ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCaga                           >  1:3277687/1‑98 (MQ=255)
                                                                    ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1837922/1‑100 (MQ=255)
                                                                        cacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       >  2:1925765/1‑100 (MQ=255)
                                                                        cacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       <  1:2165586/100‑1 (MQ=255)
                                                                        cacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       <  1:625894/100‑1 (MQ=255)
                                                                          ccGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                     >  2:2572747/1‑100 (MQ=255)
                                                                            ggTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                   <  2:207634/100‑1 (MQ=255)
                                                                            ggTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAg                   >  2:975084/1‑100 (MQ=255)
                                                                             gTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  >  2:46060/1‑100 (MQ=255)
                                                                              tCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGaa                 >  1:2858831/1‑100 (MQ=255)
                                                                               cGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAc                <  2:1276839/100‑1 (MQ=255)
                                                                                gggCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACg               >  1:1233476/1‑100 (MQ=255)
                                                                                   ccGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCtt            <  1:1399403/100‑1 (MQ=255)
                                                                                    cGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  1:789479/1‑100 (MQ=255)
                                                                                    cGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  1:1961844/100‑1 (MQ=255)
                                                                                    cGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  2:2367068/1‑100 (MQ=255)
                                                                                     gATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTgg          <  1:1053283/100‑1 (MQ=255)
                                                                                     gATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTgg          <  1:3079811/100‑1 (MQ=255)
                                                                                      atcatcAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc         >  2:2029516/1‑100 (MQ=255)
                                                                                      atcatcAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc         <  1:863626/100‑1 (MQ=255)
                                                                                       tcatcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  2:2722599/100‑1 (MQ=255)
                                                                                       tcatcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  1:2722026/100‑1 (MQ=255)
                                                                                       tcatcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  1:3096391/100‑1 (MQ=255)
                                                                                       tcatcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCg        <  2:2601805/100‑1 (MQ=255)
                                                                                          tcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:3202062/100‑1 (MQ=255)
                                                                                           caGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGt    >  1:1439963/1‑100 (MQ=255)
                                                                                            aGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc   >  2:1565547/1‑100 (MQ=255)
                                                                                             ggATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  1:2563201/1‑100 (MQ=255)
                                                                                             ggATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  1:1805554/1‑100 (MQ=255)
                                                                                             ggATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  <  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: ‑