Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,794,970 G→C 15.8% M14I (ATG→ATC waaF → ADP‑heptose:LPS heptosyltransferase II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,794,9700GC15.8% 528.4 / 84.9 196M14I (ATG→ATCwaaFADP‑heptose:LPS heptosyltransferase II
Reads supporting (aligned to +/- strand):  ref base G (72/92);  new base C (19/12);  total (91/105)
Fisher's exact test for biased strand distribution p-value = 8.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

Only 100 of 216 total aligned reads displayed.

TGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGG  >  NC_000913/3794871‑3795069
                                                                                                   |                                                                                                   
tGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatg                                                                                                     >  2:3340583/1‑100 (MQ=255)
  cTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgat                                                                                                   >  1:2712771/1‑100 (MQ=255)
   tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGGTCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                  >  1:2635858/1‑100 (MQ=255)
    gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                 >  2:226895/1‑100 (MQ=255)
    gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                 >  2:1857197/1‑100 (MQ=255)
      aagtggAACGAAATACATGGCCTGCCTGAATCGCTACGCATAAGGGCTGGGCAGCAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatgatg                                                                                               <  1:3306593/98‑1 (MQ=255)
          gtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGc                                                                                           <  2:2253062/100‑1 (MQ=255)
          gtAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGc                                                                                           >  2:3298438/1‑100 (MQ=255)
           tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                          <  2:3721068/100‑2 (MQ=255)
           tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                          <  2:1844814/100‑2 (MQ=255)
           tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                          <  1:3288260/100‑2 (MQ=255)
           tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                          <  2:111210/100‑2 (MQ=255)
           tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCACTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                          <  2:689817/100‑2 (MQ=255)
            aaCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaa                                                                                         >  1:4248/1‑100 (MQ=255)
             aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                        >  1:1571990/1‑100 (MQ=255)
              cGGAATACATGGCCTGGCTGAATCGCGACGCATAAGATCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAg                                                                                       <  2:364673/100‑1 (MQ=255)
              cGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGCcaac                                                                                       <  2:3809259/100‑2 (MQ=255)
                gAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtc                                                                                     <  1:87180/100‑1 (MQ=255)
                gAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtc                                                                                     >  1:61096/1‑100 (MQ=255)
                gAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtc                                                                                     <  2:2919572/100‑1 (MQ=255)
                  aTACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtctc                                                                                   <  1:3647546/100‑1 (MQ=255)
                    aCATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTa                                                                                 <  1:3374836/100‑1 (MQ=255)
                     cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTa                                                                                 >  1:3923668/1‑99 (MQ=255)
                     cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTAt                                                                                >  1:3465656/1‑100 (MQ=255)
                     cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGACCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTAt                                                                                >  1:2207619/1‑100 (MQ=255)
                      aTGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGCCATGATGATGTCGCAAAGTCTCTATc                                                                               <  2:2672472/100‑1 (MQ=255)
                       tGGCCTGGCTGAATCGCGACGCATAGGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCg                                                                              <  2:3056722/100‑1 (MQ=255)
                        ggCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAAGACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTCTCGc                                                                             >  1:2407496/1‑100 (MQ=255)
                         gCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGAGCGGCCCGTCTTGGGTTGTCGACATGATGAGGTCCCAAAGTCTCTATCGCa                                                                            >  2:3455904/1‑100 (MQ=255)
                         gCCTGGCTGAATCGCGACGCATAAGAGCGCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCa                                                                            <  2:1443548/100‑1 (MQ=255)
                          ccTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCGCGGCTTGGGTTGGCGACATGATGATGTCTCAAAGTCTCTATCGCAc                                                                           >  1:3211945/1‑100 (MQ=255)
                          ccTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCAc                                                                           >  2:2455029/1‑100 (MQ=255)
                          ccTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCAc                                                                           <  2:536666/100‑1 (MQ=255)
                          ccTGGCTAAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCAc                                                                           >  2:2303840/1‑100 (MQ=255)
                            tGGCTGAGTGGTGACGCATAAGAGCTCTGCATGAAAATACTGGGGATCGGCCCGTTTTTTGTTGGCGACATTATGATGTCGCAAAGTCTCTATCg                                                                              <  1:1450815/95‑1 (MQ=255)
                            tGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGc                                                                         <  2:2390230/100‑1 (MQ=255)
                             ggCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                        >  1:2989387/1‑100 (MQ=255)
                             ggCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                        <  1:3832182/100‑1 (MQ=255)
                             ggCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCt                                                                        >  2:2363753/1‑100 (MQ=255)
                                tGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                     >  2:4194492/1‑100 (MQ=255)
                                  aaTCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGCCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAgg                                                                   <  1:2812591/100‑1 (MQ=255)
                                      gcgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                               >  2:2084282/1‑100 (MQ=255)
                                      gcgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                               <  2:4066843/100‑1 (MQ=255)
                                      gcgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                               >  2:1313015/1‑100 (MQ=255)
                                       cgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcgn                                                              <  2:1140075/100‑2 (MQ=255)
                                       cgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcgc                                                              <  2:1534129/100‑1 (MQ=255)
                                        gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCt                                                             >  2:3384096/1‑100 (MQ=255)
                                         aCGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTa                                                            >  2:2614202/1‑100 (MQ=255)
                                            cATAAGAGCTCTGTATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATc                                                          <  1:3694540/99‑1 (MQ=255)
                                            cATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcc                                                         <  1:782862/100‑1 (MQ=255)
                                            cATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcc                                                         >  2:4014356/1‑100 (MQ=255)
                                              tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcccc                                                       <  1:946083/100‑1 (MQ=255)
                                              tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcccc                                                       >  2:3578040/1‑100 (MQ=255)
                                                agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                     <  2:2697602/100‑1 (MQ=255)
                                                   gCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCg                                                  >  2:1188987/1‑100 (MQ=255)
                                                   gCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCg                                                  >  1:1542571/1‑100 (MQ=255)
                                                     tctGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAt                                                <  2:3477369/100‑1 (MQ=255)
                                                       tGCATGAAAATACTGGTGATCGGCCCGTCGTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGAt                                                <  1:2086417/98‑1 (MQ=255)
                                                       tGCATGAAAATACTGCTGACCGGCCCGTCTTGGGTTGGCGACAGGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGCCGCGCGCTCCCCCCGCGAtca                                              >  1:4123630/1‑98 (MQ=255)
                                                        gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGGAAAGTCTCTATCGCACGCTCCAGGCGCGGTGTCCCCAGGCGATatt                                             >  2:820150/1‑98 (MQ=255)
                                                        gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAAt                                             >  1:1006261/1‑100 (MQ=255)
                                                        gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGAAATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAAt                                             >  1:2233966/1‑100 (MQ=255)
                                                        gCATGAAAATACTGGTGATCGGCCCGCCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATaa                                              <  1:2302750/99‑1 (MQ=255)
                                                         cATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                            >  1:2514105/1‑100 (MQ=255)
                                                         cATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                            <  2:2995993/100‑1 (MQ=255)
                                                         cATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                            <  1:820888/100‑1 (MQ=255)
                                                           tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGa                                          <  2:2168279/100‑1 (MQ=255)
                                                           tGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGa                                          >  2:2735120/1‑100 (MQ=255)
                                                            gAAAATACTGGTTATCGGCCCGTCTTGGGTTGGCGACATGATTATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                         <  1:712491/100‑1 (MQ=255)
                                                             aaaaTACTGGTGATCGGCCGGTCTTGGGTTGCCGACATGATGATGTCGGAAAGTCTCGATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtg                                        <  1:2009455/100‑1 (MQ=255)
                                                             aaaaTACTGGTGATCGGCCCGGCTTGGGTTGGCGACATGATGACGTCGCAAAGTCTCTATCGCACGCTCAAGGTGCGCTATCCCCAGGCGCTAATCGAtg                                        <  1:665753/100‑1 (MQ=255)
                                                                 tACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGAt                                    <  1:185107/100‑1 (MQ=255)
                                                                  aCTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGACGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATg                                   <  1:3029145/100‑1 (MQ=255)
                                                                   cTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATgg                                  <  2:2857155/100‑1 (MQ=255)
                                                                       tGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                              <  2:3019235/100‑1 (MQ=255)
                                                                          tCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGc                           >  2:2002632/1‑100 (MQ=255)
                                                                           cGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCa                          <  2:511510/100‑1 (MQ=255)
                                                                            ggCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCAt                         >  1:1718149/1‑100 (MQ=255)
                                                                            ggCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCAt                         >  2:1471893/1‑100 (MQ=255)
                                                                            ggCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCAt                         >  2:1958147/1‑100 (MQ=255)
                                                                             gCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATg                        <  2:1952456/100‑1 (MQ=255)
                                                                             gCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATg                        <  2:2782067/100‑1 (MQ=255)
                                                                              cccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATgg                       <  2:579496/100‑1 (MQ=255)
                                                                                cGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTg                     >  1:922341/1‑100 (MQ=255)
                                                                                cGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTg                     >  2:1759617/1‑100 (MQ=255)
                                                                                 gTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGc                    >  1:2473765/1‑100 (MQ=255)
                                                                                    ttGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGAGATGGCACCGGCATGGTGCCGt                 <  1:2550757/100‑1 (MQ=255)
                                                                                     tGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTc                >  1:3163676/1‑100 (MQ=255)
                                                                                       ggTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCa              <  2:4164926/100‑1 (MQ=255)
                                                                                        gTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat             >  2:3911289/1‑100 (MQ=255)
                                                                                        ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat             >  2:1626332/5‑100 (MQ=255)
                                                                                        ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat             >  2:1532013/5‑100 (MQ=255)
                                                                                        ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat             >  2:4346822/5‑100 (MQ=255)
                                                                                        cttagCGACATCATCATGGCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATATGAAGACACAGGCAAAGTAACGTCCa              >  1:107436/5‑99 (MQ=255)
                                                                                         tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatt            >  1:3404844/4‑100 (MQ=255)
                                                                                          ttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatta           >  1:2772031/3‑100 (MQ=255)
                                                                                          tGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCTAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatta           >  1:4386062/1‑100 (MQ=255)
                                                                                           ggCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCattat          <  1:4035241/100‑1 (MQ=255)
                                                                                                  tgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATcgcg   <  2:83420/100‑1 (MQ=255)
                                                                                                   gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCgg  >  1:1363087/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TGCTGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGG  >  NC_000913/3794871‑3795069

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: