Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,121,597 (T)5→4 100% 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.97.2% 758.8 / 8.2 173coding (155/1332 nt)hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base . (87/81);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 190 total aligned reads displayed.

TTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCC  >  NC_000913/4121494‑4121688
                                                                                                   |                                                                                               
ttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCGGATCATCAGGatat                                                                                                 >  2:2604341/1‑100 (MQ=255)
  ggTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttt                                                                                               >  1:1334860/1‑100 (MQ=255)
   gTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtttt                                                                                              >  1:2093508/1‑100 (MQ=255)
     cGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcg                                                                                            <  1:1406023/100‑3 (MQ=255)
      gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGCTATTTTCCg                                                                                           >  2:3107518/1‑100 (MQ=255)
      gCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGGCGGGCCGATGATCAGGATAttttcgg                                                                                           >  1:1944624/1‑97 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  1:1176435/100‑5 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          >  2:463605/1‑96 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          >  1:3211156/1‑96 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  1:2015206/100‑5 (MQ=255)
        ttCAACTTTGATGCCCGGCGCCTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCgggg                                                                                        <  1:1127337/100‑1 (MQ=255)
        ttCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCgggg                                                                                        <  1:2382041/100‑1 (MQ=255)
         tccaCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcgg                                                                                           >  2:3615082/4‑94 (MQ=255)
         tCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGt                                                                                       >  1:1400203/1‑100 (MQ=255)
            aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCAc                                                                                    <  1:2665068/100‑1 (MQ=255)
            aCTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCAc                                                                                    >  2:372808/1‑100 (MQ=255)
                tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCa                                                                                >  2:1785146/1‑100 (MQ=255)
                tgatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCa                                                                                >  2:3700557/1‑100 (MQ=255)
                 gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCAt                                                                               <  2:2425823/100‑1 (MQ=255)
                  atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCCCTTCAt                                                                               >  1:440589/1‑99 (MQ=255)
                   tgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATgg                                                                             >  2:2869158/1‑100 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgc                                                                            >  2:575485/1‑100 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgc                                                                            >  1:1333878/1‑100 (MQ=255)
                     aaCGGCGCATTCGCCAGCTTAGCCAGACGACGGTCGAATTAAGTTTTACCGACACCGGTCGGGTCTATCATCAGGATA‑TTTTCGGGGTCACTTCATgg                                                                             <  2:2688109/98‑1 (MQ=255)
                      aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCCGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgcgc                                                                          <  1:3043680/100‑1 (MQ=255)
                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGCTTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCTGGGTCACTTCATGGcgccgc                                                                       >  2:1059571/1‑100 (MQ=255)
                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGc                                                                       <  2:1467508/100‑1 (MQ=255)
                         gcgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGc                                                                       >  2:1080067/1‑100 (MQ=255)
                          cgcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGc                                                                       <  2:1224274/99‑1 (MQ=255)
                           gcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGctc                                                                     >  2:946744/1‑100 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  1:1083650/1‑100 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  1:1033251/1‑100 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTGTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  2:2033638/1‑100 (MQ=255)
                             aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   >  2:886786/1‑100 (MQ=255)
                               tCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCg                                                                 <  1:1778630/100‑1 (MQ=255)
                                cGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGCTCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGt                                                                >  1:3234742/1‑100 (MQ=255)
                                  ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  1:186508/100‑1 (MQ=255)
                                  ccAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTg                                                              <  2:152931/100‑1 (MQ=255)
                                   cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  1:1547850/100‑1 (MQ=255)
                                   cAGCTTAGCCAGACGACGCGCTATTTCATTTTTCCCGACACCGGTCGGGCCGAGCATCAGGATA‑TTTACGGGGTCACTTCATGGCGCAGCTCTTCGt                                                                <  2:3225751/97‑1 (MQ=255)
                                    aGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAg                                                            >  1:3196032/1‑100 (MQ=255)
                                    aGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAg                                                            >  1:1292860/1‑100 (MQ=255)
                                         aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  1:3328076/100‑1 (MQ=255)
                                          gCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCAt                                                      >  1:754738/1‑100 (MQ=255)
                                          cccAGACGACGGGCGCTTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGAT‑CGTTTCGGGGTCACTTCATGCCGCAGCTCTTCGTTCAGCTGcct                                                      >  1:1761774/2‑98 (MQ=255)
                                           ccAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                     <  2:2788040/100‑1 (MQ=255)
                                            cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGc                                                    >  1:3095584/1‑100 (MQ=255)
                                            cAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGc                                                    >  1:2296523/1‑100 (MQ=255)
                                             aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  2:1448417/1‑100 (MQ=255)
                                                cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACg                                                >  1:2486032/1‑100 (MQ=255)
                                                 gacgGGCGATTTCAGTTTTACCGGCACCGGTCGGGCCGATCATCAGGAT‑GTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  2:361903/98‑1 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGcc                                              >  2:2494640/1‑100 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGc                                               <  2:1869380/99‑1 (MQ=255)
                                                   cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCa                                             >  1:3408909/1‑100 (MQ=255)
                                                   cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCa                                             >  2:938921/1‑100 (MQ=255)
                                                    gggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                            <  1:2853873/100‑1 (MQ=255)
                                                    gggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAg                                            >  2:521318/1‑100 (MQ=255)
                                                     ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGCCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                           >  1:675166/1‑100 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        <  2:3403080/100‑1 (MQ=255)
                                                          tttCAGTTTTGCCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTGCATGGCGCGGCTCTTCGTTGAGCTTCATGGGACGCCAGGGGTTa                                      >  1:3627734/1‑100 (MQ=255)
                                                             cAGTTTTACCGACACCGGTCGGGTCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                    <  1:2986609/99‑1 (MQ=255)
                                                             cAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                   <  2:1373767/100‑1 (MQ=255)
                                                              aGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCa                                  <  1:3675327/100‑1 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGc                                   <  2:3599710/98‑1 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  1:1427368/1‑100 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  1:3340762/1‑100 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 >  1:3136457/1‑100 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               <  1:902540/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               >  2:2381821/1‑100 (MQ=255)
                                                                  ttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                              <  1:3410567/100‑1 (MQ=255)
                                                                  ttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGc                              <  1:968863/100‑1 (MQ=255)
                                                                   tACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCa                             <  1:2712171/100‑1 (MQ=255)
                                                                    aCCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCaa                            >  1:3312472/1‑100 (MQ=255)
                                                                    aCCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCaa                            >  2:2893496/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  1:414622/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1988447/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  1:498698/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1201352/1‑100 (MQ=255)
                                                                      cGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATc                          >  2:2715248/1‑100 (MQ=255)
                                                                       gACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                         >  2:2256666/1‑100 (MQ=255)
                                                                         cacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       >  2:381787/1‑100 (MQ=255)
                                                                         cacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGcc                       >  1:3175449/1‑100 (MQ=255)
                                                                          acCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCca                      >  2:343478/1‑100 (MQ=255)
                                                                            cGGTCGGGCCGATCATCAGGATA‑TTTGCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                    <  2:2995072/100‑1 (MQ=255)
                                                                             ggTCGGGCAGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCACATCGTTGAGATGAACGCTACGCCATCGGTTACGCAGAGCAATCGCCCCAg                   >  1:1616668/1‑100 (MQ=255)
                                                                              gTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  <  1:1047704/100‑1 (MQ=255)
                                                                               tCGGGCCGATCATCAGGTTATTTTTCGGGGTCACTTCATGGCGCTGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                     <  2:3185825/97‑1 (MQ=255)
                                                                                  ggCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc              <  2:1794753/100‑1 (MQ=255)
                                                                                     cGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  2:1334851/100‑1 (MQ=255)
                                                                                     cGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  1:3511953/100‑1 (MQ=255)
                                                                                     cGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           <  1:209191/100‑1 (MQ=255)
                                                                                      gATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGCCGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  2:2712001/1‑99 (MQ=255)
                                                                                         catcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGt       >  1:3053430/1‑100 (MQ=255)
                                                                                          atcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGtt      >  2:75171/1‑100 (MQ=255)
                                                                                           tcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGCGCAATCGCCACAGAACGCTTcgcg        >  2:1267291/1‑97 (MQ=255)
                                                                                           tcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:2031550/100‑1 (MQ=255)
                                                                                           tcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:1042778/100‑1 (MQ=255)
                                                                                             aGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc   <  1:847069/100‑1 (MQ=255)
                                                                                              ggATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  1:60689/1‑100 (MQ=255)
                                                                                              ggATA‑TTTTCGGGGTCACGTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  2:1389466/1‑100 (MQ=255)
                                                                                                   |                                                                                               
TTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCC  >  NC_000913/4121494‑4121688

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