Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 681,694 A→G 15.2% intergenic (+31/+29) djlC → / ← hscC J domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913681,6940AG15.2% 330.1 / 49.5 138intergenic (+31/+29)djlC/hscCJ domain‑containing HscC co‑chaperone; Hsc56/Hsp70 family chaperone Hsc62; RpoD‑binding transcription inhibitor
Reads supporting (aligned to +/- strand):  ref base A (67/50);  new base G (6/15);  total (73/65)
Fisher's exact test for biased strand distribution p-value = 1.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.28e-01

Only 100 of 144 total aligned reads displayed.

ATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTG  >  NC_000913/681595‑681793
                                                                                                   |                                                                                                   
 tCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAg                                                                                                    <  2:162775/100‑1 (MQ=255)
  cAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGc                                                                                                   <  1:3643372/100‑1 (MQ=255)
  cAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGc                                                                                                   <  2:3416353/100‑1 (MQ=255)
       tCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTa                                                                                              <  1:1895898/100‑1 (MQ=255)
        cAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAc                                                                                             >  2:2827478/1‑100 (MQ=255)
          gTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGCTGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgc                                                                                           >  1:2301857/1‑100 (MQ=255)
          gTTTAACGGTGCTGGTAATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTTCGCCGATGCCGGATGCAGCGTTTAc                                                                                             >  1:575025/1‑98 (MQ=255)
           ttccgatctGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgcc                                                                                          <  1:3635322/92‑1 (MQ=255)
            ttATCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTAcgccg                                                                                         <  1:3619669/100‑1 (MQ=255)
            ttAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGtaaacgctg                                                                                         <  2:3040891/100‑9 (MQ=255)
              aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCa                                                                                       >  2:1097900/1‑100 (MQ=255)
              aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCa                                                                                       <  1:173527/100‑1 (MQ=255)
              aaCGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCa                                                                                       >  1:1096089/1‑100 (MQ=255)
                 ggTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATcc                                                                                    <  2:698187/100‑1 (MQ=38)
                  gTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCg                                                                                   <  2:262974/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  2:2371430/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  2:2081621/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  2:2007993/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  1:222393/100‑1 (MQ=38)
                   tGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCgg                                                                                  <  2:3018082/100‑1 (MQ=38)
                    gCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGTTGCCGGATGCGGCGTAATCGCTGCATCCgcc                                                                                 <  2:2491947/100‑3 (MQ=38)
                    gCTGGTGATAGGAAATAAACACAGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCgg                                                                                  <  1:232255/99‑1 (MQ=255)
                       ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTTCTTTGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGc                                                                              >  1:1240200/1‑100 (MQ=255)
                       ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGCatc                                                                              <  1:2965354/100‑3 (MQ=37)
                       ggTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGc                                                                              >  2:2635643/1‑100 (MQ=255)
                        gTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTGTTGCGCCGATGCCGGATGCGGCGTAAACGCTGCATCCGGCatcg                                                                             <  2:240669/100‑4 (MQ=37)
                         tGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAg                                                                               >  1:781126/1‑97 (MQ=255)
                         tGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGAAGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCg                                                                            >  1:3330400/1‑100 (MQ=255)
                          gATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGAGTTTACGCCGCATCCGGCAGCCgt                                                                           >  2:1470329/1‑100 (MQ=255)
                              ggTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGttt                                                                       >  2:155402/1‑100 (MQ=255)
                               gTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTa                                                                      >  1:2451103/1‑100 (MQ=255)
                               gTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTa                                                                      <  2:3588776/100‑1 (MQ=255)
                                 tATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATg                                                                    <  2:2372718/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                  <  2:152424/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                  <  2:3300477/100‑1 (MQ=255)
                                   tGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATggg                                                                  <  1:876166/100‑1 (MQ=255)
                                   tGCACTCGGCAAAATCCCGATGATGTTTAGCTGACGATGTGTTTTTTGCGCTGATGCCGGATGAAACGTTTAAGCCGCATTATGCAGCCGTttttcacg                                                                   >  1:2478679/1‑95 (MQ=255)
                                    gCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCTTCTTTTCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                 >  1:3548736/1‑100 (MQ=38)
                                    gCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGa                                                                 >  2:886998/1‑100 (MQ=255)
                                        tcGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGAtctc                                                             >  1:2721424/1‑100 (MQ=255)
                                         cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                            <  1:673186/100‑1 (MQ=255)
                                         cGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCa                                                            >  1:448508/1‑100 (MQ=255)
                                           gCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAAt                                                          <  2:3627608/100‑1 (MQ=255)
                                            cAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATg                                                         >  1:1884732/1‑100 (MQ=255)
                                             aaaaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATgg                                                        >  2:3450938/1‑100 (MQ=255)
                                               aaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGca                                                      >  2:2286439/1‑99 (MQ=255)
                                               aaTCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCt                                                      >  1:2612795/1‑100 (MQ=255)
                                                  cgcgATGTTGTTTAGCTGACGATGTGTTTTTTGCTCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGAGTTTATGGGATCTCAATGGATaaa                                                   >  1:2181977/1‑100 (MQ=255)
                                                    cgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGCTGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaata                                                 >  2:30671/1‑100 (MQ=255)
                                                    cgATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaata                                                 >  2:2940034/1‑100 (MQ=255)
                                                    cgATGTTGGTTAGGTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAaata                                                 <  1:2759800/100‑1 (MQ=255)
                                                       tgttgtTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                              <  2:1752317/100‑1 (MQ=255)
                                                       tgttgtTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                              <  1:3468933/100‑1 (MQ=255)
                                                       gttgttggggCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtc                                              <  2:216607/91‑1 (MQ=255)
                                                         ttgttTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgc                                            <  2:1679858/100‑1 (MQ=255)
                                                             ttAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAAtcgctcgc                                        <  2:3272759/100‑1 (MQ=255)
                                                              tAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcc                                       >  1:2990508/1‑100 (MQ=255)
                                                               aGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgccc                                      <  2:1888328/100‑1 (MQ=255)
                                                                gCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccg                                     >  2:787329/1‑100 (MQ=255)
                                                                  tGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTcgcccgcc                                   <  2:3117510/100‑1 (MQ=255)
                                                                    aCGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgc                                 >  2:2723441/1‑100 (MQ=255)
                                                                     cGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgcc                                <  1:1099743/100‑1 (MQ=255)
                                                                      tctgtgTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCccgccgccg                               <  1:3200088/98‑1 (MQ=255)
                                                                       aTGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGa                              >  1:3101940/1‑100 (MQ=255)
                                                                         gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGaga                            >  2:3200112/1‑98 (MQ=255)
                                                                         gtgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAAc                            <  2:51500/100‑1 (MQ=255)
                                                                          tgtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACt                           >  2:55054/1‑100 (MQ=255)
                                                                           gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                          >  1:3231323/1‑100 (MQ=255)
                                                                           gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACtt                          >  2:1218466/1‑100 (MQ=255)
                                                                           gtTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTCATGGGATCTCAATGGCTAAATAATCGCTCGCCCCCCGCCGAAc                            >  1:1385269/1‑98 (MQ=255)
                                                                            tgttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTc                         <  1:3410822/98‑1 (MQ=255)
                                                                             tttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCa                        <  2:585287/100‑1 (MQ=255)
                                                                             tttttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCa                        <  1:2042071/100‑1 (MQ=255)
                                                                               tttGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATc                      <  2:1368610/100‑1 (MQ=255)
                                                                                 tGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCa                    >  2:1312667/1‑100 (MQ=255)
                                                                                 tGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCa                    >  1:2162181/1‑100 (MQ=255)
                                                                                  gcgcCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCAc                   >  1:195305/1‑100 (MQ=255)
                                                                                       gATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCggg              <  1:826354/100‑1 (MQ=255)
                                                                                          tccgatctcAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           <  1:764566/92‑1 (MQ=255)
                                                                                          gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           >  1:1061227/1‑100 (MQ=255)
                                                                                          gCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           <  2:1061233/100‑1 (MQ=255)
                                                                                           ccGGATGCGGCGTAAACGCTGCATGCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGATCGCCCGCCCCCGAACTTCATCCACGCGGGCCat          >  1:1778737/1‑100 (MQ=37)
                                                                                           ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat          <  2:838968/100‑1 (MQ=255)
                                                                                           ccGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat          >  2:1834635/1‑100 (MQ=255)
                                                                                            cGGATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata         >  1:2548433/1‑100 (MQ=38)
                                                                                            cGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATACTCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCa           >  2:2533218/1‑98 (MQ=255)
                                                                                            cGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata         <  1:2940013/100‑1 (MQ=255)
                                                                                            cGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCata         >  2:3030798/1‑100 (MQ=255)
                                                                                             ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        >  2:3489757/1‑100 (MQ=38)
                                                                                             ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        >  2:1457559/1‑100 (MQ=38)
                                                                                             ggATGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCat          >  1:454218/1‑98 (MQ=38)
                                                                                             ggATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCatat        >  2:387634/1‑100 (MQ=255)
                                                                                              gATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATc       <  1:886993/100‑1 (MQ=255)
                                                                                              gATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATc       >  2:226313/1‑100 (MQ=255)
                                                                                               aTGCGGCGTAAACGCTGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCg      >  2:1854779/1‑100 (MQ=39)
                                                                                                 gCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGtt    >  2:698943/1‑100 (MQ=255)
                                                                                                  cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCaagatcgga   >  2:764571/1‑92 (MQ=255)
                                                                                                  cAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGttt   >  2:1198524/1‑100 (MQ=255)
                                                                                                   aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  2:374985/100‑1 (MQ=255)
                                                                                                   aGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTg  <  2:1559322/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ATCAAAATCAGTTTAACGGTGCTGGTGATAGGTTATGCTCTCGGCAAAATCGCGATGTTGTTTAGCTGACGATGTGTTTTTTGCGCCGATGCCGGATGCAGCGTTTACGCCGCATCCGGCAGCCGTGTTTATGGGATCTCAATGGCTAAATAATCGCTCGCCCGCCGCCGAACTTCATCCACGCGGGCCATATCGTTTG  >  NC_000913/681595‑681793

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: