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)
         tCAACTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGt                                                                                       >  1:1204415/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCAc                                                                                    >  2:460943/1‑99 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACt                                                                                   >  1:237828/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACt                                                                                   >  2:272385/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACt                                                                                   >  2:439706/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACt                                                                                   >  1:199954/1‑100 (MQ=255)
             cTTTGATGAACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACt                                                                                   >  2:1718248/1‑100 (MQ=255)
               ttccgatctggCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTc                                                                                 <  1:119956/91‑1 (MQ=255)
                 gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCa                                                                                <  2:148338/99‑1 (MQ=255)
                 gatgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCAt                                                                               >  2:1582931/1‑100 (MQ=255)
                  atgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATg                                                                              >  1:1695078/1‑100 (MQ=255)
                   tgaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATCTTTTTCGGGGAAAcatcat                                                                               >  1:1387871/1‑99 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgc                                                                            <  1:1132713/100‑1 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgc                                                                            <  1:1485789/100‑1 (MQ=255)
                    gaACGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgc                                                                            >  2:1769012/1‑100 (MQ=255)
                     aaCGGCGTATTCGCCAGCTTAGCCAGACGACGGGCTATTTCAGTTTTACCGACACCGGGCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATgg                                                                             <  2:430591/98‑1 (MQ=255)
                     aaCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgcg                                                                           <  1:1717081/100‑1 (MQ=255)
                      aCGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGgcgc                                                                          >  2:1685137/1‑100 (MQ=255)
                       cGGCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCa                                                                         <  1:1331055/100‑1 (MQ=255)
                        ggCGCATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCagatcggaa                                                                        >  2:119955/1‑92 (MQ=255)
                           gcATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGctc                                                                     <  2:1769072/100‑1 (MQ=255)
                            cATTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TGTTCGGGGTCACTTCATGGCGCAGctct                                                                    >  1:2040396/1‑100 (MQ=255)
                             aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   <  1:1190026/100‑1 (MQ=255)
                             aTTCGCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCtt                                                                   >  2:599276/1‑100 (MQ=255)
                                 gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                               >  2:557951/1‑100 (MQ=255)
                                 gCCAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGtt                                                               >  2:111938/1‑100 (MQ=255)
                                   cAGCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGa                                                             <  2:1974828/100‑1 (MQ=255)
                                     gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGctcaacgaagagc                                                           <  2:247801/100‑11 (MQ=255)
                                     gCTTAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGc                                                           <  1:373054/100‑1 (MQ=255)
                                       ttAGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTg                                                         <  1:447350/100‑1 (MQ=255)
                                         aGCCAGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  1:1639086/100‑1 (MQ=255)
                                           ccAGGCGACGGGCGATTACAGTTTTCCCGACACCGGTCGGGCCGAGCATCAGGATA‑TTGTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATg                                                     <  2:700812/100‑1 (MQ=255)
                                             aGACGACGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCg                                                   >  1:1007488/1‑100 (MQ=255)
                                              gacgatGGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCa                                                       <  2:999384/96‑1 (MQ=255)
                                              gacgacGGGCGATTTCAGTTTTCCCGACACCGGGCGGGCCGATCATCAGGATATTGCCCGGGG‑CACTTCATGGCGCAGCTCAACGAAGAGCTGcgcca                                                    <  2:723371/98‑5 (MQ=255)
                                               acgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  2:1167932/100‑1 (MQ=255)
                                               acgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  1:46588/100‑1 (MQ=255)
                                                cgacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGAc                                                 <  1:1987901/99‑1 (MQ=255)
                                                 gacgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGc                                               <  1:1879156/100‑1 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGcc                                              >  1:1421528/1‑100 (MQ=255)
                                                  acgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGAAGCTCTTCGTTGAGCTGCATGCGACGc                                               >  1:542905/1‑99 (MQ=255)
                                                   cgGGCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCa                                             >  1:2076161/1‑100 (MQ=255)
                                                     ggCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCACGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGc                                           >  2:756213/1‑100 (MQ=255)
                                                      gCGATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCg                                          >  1:1441461/1‑100 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        >  1:806347/1‑100 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        <  1:2299544/100‑1 (MQ=255)
                                                        gATTTCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGt                                        >  2:1620001/1‑100 (MQ=255)
                                                            tCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                    >  1:977266/1‑100 (MQ=255)
                                                            tCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACg                                    >  1:977154/1‑100 (MQ=255)
                                                            tCAGTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTGGAGCTGCATGCGTCGCCAGCGGTTACg                                    >  2:2166456/1‑100 (MQ=255)
                                                               gTTTTACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCag                                 <  2:2225534/100‑1 (MQ=255)
                                                                ttttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCaga                                <  2:237826/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               <  1:1812114/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               <  1:357066/100‑1 (MQ=255)
                                                                 tttACCGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCagag                               >  2:2065197/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATATTTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCa                             >  1:1645302/1‑99 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCGCTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:1295734/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:317495/1‑100 (MQ=255)
                                                                     ccGACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAAt                           >  2:836303/1‑100 (MQ=255)
                                                                       gACACCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCg                         >  2:2374864/1‑100 (MQ=255)
                                                                        acacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                        >  2:585350/1‑100 (MQ=255)
                                                                        acacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                        <  2:1332634/100‑1 (MQ=255)
                                                                        acacCGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGc                        <  2:1723788/100‑1 (MQ=255)
                                                                           ccGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcac                     >  1:927791/1‑100 (MQ=255)
                                                                            cGGTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCcaca                    <  2:622510/100‑1 (MQ=255)
                                                                              gTCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  >  2:1518798/1‑100 (MQ=255)
                                                                               tCGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGa                  <  2:1573623/99‑1 (MQ=255)
                                                                                cGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAc                <  1:2095448/100‑1 (MQ=255)
                                                                                cGGGCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAAc                >  2:1322209/1‑100 (MQ=255)
                                                                                  ggCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc              >  1:1605861/1‑100 (MQ=255)
                                                                                  ggCCGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGc              >  2:1922981/1‑100 (MQ=255)
                                                                                    ccGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCtt            >  1:1089289/1‑100 (MQ=255)
                                                                                     cGATCATCAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTg           >  1:1416816/1‑100 (MQ=255)
                                                                                       atcatcAGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGc         <  2:2066170/100‑1 (MQ=255)
                                                                                          atcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGGGCAATCGCCACAGAACGCTTGGCGtt      <  2:2155988/100‑1 (MQ=255)
                                                                                           tcaGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTg     <  1:2060840/100‑1 (MQ=255)
                                                                                             aGGATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTc   <  1:1620013/100‑1 (MQ=255)
                                                                                              ggATA‑TTTTCGGGGTCACTTCATGGCGCAGCTCTTCGTTGAGCTGCATGCGACGCCAGCGGTTACGCAGAGCAATCGCCACAGAACGCTTGGCGTTGTcc  >  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: ‑