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.94.9% 320.6 / 10.5 78coding (155/1332 nt)hslUmolecular chaperone and ATPase component of HslUV protease
Reads supporting (aligned to +/- strand):  ref base T (2/2);  new base . (41/33);  total (43/35)
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.88e-01

TTGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTCC  >  NC_000913/4121494‑4121688
                                                                                                   |                                                                                               
ttGGTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGatat                                                                                                 <  1:635123/100‑1 (MQ=255)
 tGGTCGCTTCAACATTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAtt                                                                                                >  1:1621549/1‑100 (MQ=255)
  ggTCGCTTCAACTTTGATGAACGGCGCATTCGCCCGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGCCCCCGCTCCGGCCGATCATCCGGATAttt                                                                                               >  2:1001098/1‑100 (MQ=255)
  ggTCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttt                                                                                               >  1:160898/1‑100 (MQ=255)
    tCGCTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttc                                                                                             >  2:1758712/1‑99 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  1:268333/100‑5 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          >  1:175690/1‑96 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  2:198178/100‑5 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  1:234861/100‑5 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          <  1:65359/100‑5 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          >  1:1373216/1‑96 (MQ=255)
       cTTCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATAttttcggg                                                                                          >  1:706713/1‑96 (MQ=255)
         tCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGt                                                                                       >  1:1204415/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCAc                                                                                    >  2:460943/1‑99 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACt                                                                                   >  1:237828/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACt                                                                                   >  2:272385/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACt                                                                                   >  2:439706/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACt                                                                                   >  1:199954/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACt                                                                                   >  2:1718248/1‑100 (MQ=255)
               ttccgatctggCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTc                                                                                 <  1:119956/91‑1 (MQ=255)
                 gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCa                                                                                <  2:148338/99‑1 (MQ=255)
                 gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCAt                                                                               >  2:1582931/1‑100 (MQ=255)
                  atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATg                                                                              >  1:1695078/1‑100 (MQ=255)
                   tgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATCTTTTTCGGGGAAAcatcat                                                                               >  1:1387871/1‑99 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgc                                                                            <  1:1132713/100‑1 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgc                                                                            <  1:1485789/100‑1 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgc                                                                            >  2:1769012/1‑100 (MQ=255)
                     aaCGGCGTATTCGCCAGCTTAGCCAGACGACGGGCTATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATgg                                                                             <  2:430591/98‑1 (MQ=255)
                     aaCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgcg                                                                           <  1:1717081/100‑1 (MQ=255)
                      aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGgcgc                                                                          >  2:1685137/1‑100 (MQ=255)
                       cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCa                                                                         <  1:1331055/100‑1 (MQ=255)
                        ggCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCagatcggaa                                                                        >  2:119955/1‑92 (MQ=255)
                           gcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGctc                                                                     <  2:1769072/100‑1 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATGTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  1:2040396/1‑100 (MQ=255)
                             aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   <  1:1190026/100‑1 (MQ=255)
                             aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   >  2:599276/1‑100 (MQ=255)
                                 gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                               >  2:557951/1‑100 (MQ=255)
                                 gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                               >  2:111938/1‑100 (MQ=255)
                                   cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  2:1974828/100‑1 (MQ=255)
                                     gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGctcaacgaagagc                                                           <  2:247801/100‑11 (MQ=255)
                                     gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                           <  1:373054/100‑1 (MQ=255)
                                       ttAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTg                                                         <  1:447350/100‑1 (MQ=255)
                                         aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  1:1639086/100‑1 (MQ=255)
                                           ccAGGCGACGGGCGATTACAGTTTTCCCGACACCGGTCGGGCCGAGCATCAGGATATTGTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                     <  2:700812/100‑1 (MQ=255)
                                             aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  1:1007488/1‑100 (MQ=255)
                                              gacgatGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  2:999384/96‑1 (MQ=255)
                                              gacgacGGGCGATTTCAGTTTTCCCGACACCGGGCGGGCCGATCATCAGGATATTGCCCGGGGCACTTCATGGCGCAGCTCAACGAAGAGCTGcgcca                                                    <  2:723371/98‑5 (MQ=255)
                                               acgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  2:1167932/100‑1 (MQ=255)
                                               acgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  1:46588/100‑1 (MQ=255)
                                                cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  1:1987901/99‑1 (MQ=255)
                                                 gacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGc                                               <  1:1879156/100‑1 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGcc                                              >  1:1421528/1‑100 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGAAGCTCTTCGTTGAGCTGCATGCGACGc                                               >  1:542905/1‑99 (MQ=255)
                                                   cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCa                                             >  1:2076161/1‑100 (MQ=255)
                                                     ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCACGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                           >  2:756213/1‑100 (MQ=255)
                                                      gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCg                                          >  1:1441461/1‑100 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        >  1:806347/1‑100 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        <  1:2299544/100‑1 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        >  2:1620001/1‑100 (MQ=255)
                                                            tCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                    >  1:977266/1‑100 (MQ=255)
                                                            tCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                    >  1:977154/1‑100 (MQ=255)
                                                            tCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTGGAGCTGCATGCGTCGCCAGCGGTTACg                                    >  2:2166456/1‑100 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 <  2:2225534/100‑1 (MQ=255)
                                                                ttttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCaga                                <  2:237826/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               <  1:1812114/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               <  1:357066/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               >  2:2065197/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCa                             >  1:1645302/1‑99 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCGCTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1295734/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:317495/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:836303/1‑100 (MQ=255)
                                                                       gACACCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                         >  2:2374864/1‑100 (MQ=255)
                                                                        acacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                        >  2:585350/1‑100 (MQ=255)
                                                                        acacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                        <  2:1332634/100‑1 (MQ=255)
                                                                        acacCGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                        <  2:1723788/100‑1 (MQ=255)
                                                                           ccGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                     >  1:927791/1‑100 (MQ=255)
                                                                            cGGTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                    <  2:622510/100‑1 (MQ=255)
                                                                              gTCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  >  2:1518798/1‑100 (MQ=255)
                                                                               tCGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  <  2:1573623/99‑1 (MQ=255)
                                                                                cGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAc                <  1:2095448/100‑1 (MQ=255)
                                                                                cGGGCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAc                >  2:1322209/1‑100 (MQ=255)
                                                                                  ggCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc              >  1:1605861/1‑100 (MQ=255)
                                                                                  ggCCGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc              >  2:1922981/1‑100 (MQ=255)
                                                                                    ccGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCtt            >  1:1089289/1‑100 (MQ=255)
                                                                                     cGATCATCAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  1:1416816/1‑100 (MQ=255)
                                                                                       atcatcAGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc         <  2:2066170/100‑1 (MQ=255)
                                                                                          atcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGGGCAATCGCCACAGAACGCTTGGCGtt      <  2:2155988/100‑1 (MQ=255)
                                                                                           tcaGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:2060840/100‑1 (MQ=255)
                                                                                             aGGATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc   <  1:1620013/100‑1 (MQ=255)
                                                                                              ggATATTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  2:756792/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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: