Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,794,973 G→C 15.8% M15I (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,9730GC15.8% 521.8 / 87.2 196M15I (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.40e-01

Only 100 of 219 total aligned reads displayed.

TGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATG  >  NC_000913/3794874‑3795072
                                                                                                   |                                                                                                   
tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGGTCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                     >  1:2635858/1‑100 (MQ=255)
tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                     >  2:3315513/1‑100 (MQ=255)
tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                     >  2:2720498/1‑100 (MQ=255)
tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                     >  2:2356475/1‑100 (MQ=255)
tGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatg                                                                                                     >  2:2381681/1‑100 (MQ=255)
 gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                    >  1:893792/1‑100 (MQ=255)
 gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                    >  1:314234/1‑100 (MQ=255)
 gAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACatgatga                                                                                                    >  2:3428085/1‑100 (MQ=255)
   aGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTCGCGACatgatgatg                                                                                                  >  1:786790/1‑100 (MQ=255)
    ggTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGt                                                                                                 <  2:2445656/100‑1 (MQ=255)
        tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                             <  1:3288260/100‑2 (MQ=255)
        tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACCGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                             <  1:3407702/100‑2 (MQ=255)
        tAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCACTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGcc                                                                                             <  2:689817/100‑2 (MQ=255)
         aaCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaa                                                                                            >  2:1913989/1‑100 (MQ=255)
          aCGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCaaa                                                                                           >  1:1571990/1‑100 (MQ=255)
           cGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATCATCATGTCGCcaac                                                                                          <  1:3688883/100‑2 (MQ=255)
             gAATACATGGCCTGGCTGAATCGCGACGCATAAGATCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtc                                                                                        <  1:4238953/100‑1 (MQ=255)
             gAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtc                                                                                        <  1:87180/100‑1 (MQ=255)
             gAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtc                                                                                        <  2:2919572/100‑1 (MQ=255)
               aTACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGtctc                                                                                      <  1:2524725/100‑1 (MQ=255)
                  cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTa                                                                                    <  1:4100984/99‑1 (MQ=255)
                  cATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTAt                                                                                   >  2:3680770/1‑100 (MQ=255)
                    tGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGGCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCg                                                                                 >  1:556527/1‑100 (MQ=255)
                     ggCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAAGACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTCTCGc                                                                                >  1:2407496/1‑100 (MQ=255)
                      gCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCa                                                                               >  1:3185454/1‑100 (MQ=255)
                      gCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGACCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCa                                                                               <  2:2207627/100‑1 (MQ=255)
                       ccTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCAc                                                                              <  2:1991307/100‑1 (MQ=255)
                       ccTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCAc                                                                              <  2:536666/100‑1 (MQ=255)
                        cTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACg                                                                             >  1:1755436/1‑100 (MQ=255)
                         tGGCTGAGTGGTGACGCATAAGAGCTCTGCATGAAAATACTGGGGATCGGCCCGTTTTTTGTTGGCGACATTATGATGTCGCAAAGTCTCTATCg                                                                                 <  1:1450815/95‑1 (MQ=255)
                         tGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGc                                                                            <  1:20661/100‑1 (MQ=255)
                            cTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTcc                                                                         >  1:354896/1‑100 (MQ=255)
                             tGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                        >  1:674959/1‑100 (MQ=255)
                             tGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                        <  1:3400825/100‑1 (MQ=255)
                             tGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCa                                                                        >  2:4194492/1‑100 (MQ=255)
                               aaTCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGCCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAgg                                                                      <  1:2812591/100‑1 (MQ=255)
                                 tCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcg                                                                    <  1:4027690/100‑1 (MQ=255)
                                   gcgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                  <  2:4066843/100‑1 (MQ=255)
                                   gcgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcg                                                                  >  2:1313015/1‑100 (MQ=255)
                                    cgACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGgcgcgc                                                                 <  2:1534129/100‑1 (MQ=255)
                                     gACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCt                                                                >  2:3384096/1‑100 (MQ=255)
                                           tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcccc                                                          <  2:1790432/100‑1 (MQ=255)
                                           tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcccc                                                          >  2:3578040/1‑100 (MQ=255)
                                           tAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATcccc                                                          >  2:1502770/1‑100 (MQ=255)
                                            aaGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCa                                                         >  2:2284125/1‑100 (MQ=255)
                                             agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                        >  2:1112892/1‑100 (MQ=255)
                                             agagCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAg                                                        >  2:4363608/1‑100 (MQ=255)
                                                gCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCg                                                     >  2:1188987/1‑100 (MQ=255)
                                                gCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCg                                                     >  1:1542571/1‑100 (MQ=255)
                                                 ctctGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGa                                                    >  1:1997201/1‑100 (MQ=255)
                                                     gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGGAAAGTCTCTATCGCACGCTCCAGGCGCGGTGTCCCCAGGCGATatt                                                >  2:820150/1‑98 (MQ=255)
                                                     gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAAt                                                >  2:2148727/1‑100 (MQ=255)
                                                     gCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAAt                                                >  1:1006261/1‑100 (MQ=255)
                                                      cATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATc                                               <  2:2995993/100‑1 (MQ=255)
                                                         gAAAATACTGGTGATCGGCCCGTCTTGCGTTGCCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAt                                            <  1:712934/100‑1 (MQ=255)
                                                          aaaaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtg                                           <  1:2735105/100‑1 (MQ=255)
                                                          aaaaTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtg                                           <  2:3892351/100‑1 (MQ=255)
                                                          aaaaTACTGGTGATCGGCCCGGCTTGGGTTGGCGACATGATGATGTCGCAAAGTCGCTAGCGCACGCGCCAGGCGCGCTATCCCCAGGCGATAATCGAtg                                           <  1:665612/100‑1 (MQ=255)
                                                            aaTACTGGTGATCCGCCCGTCTTGTGTGGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGAtgtg                                         <  1:3819336/100‑1 (MQ=255)
                                                             aTACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGa                                        <  2:1743522/100‑1 (MQ=255)
                                                               aCTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATg                                      <  2:1464472/100‑1 (MQ=255)
                                                                cTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATgg                                     <  2:2857155/100‑1 (MQ=255)
                                                                   gTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAc                                  >  2:3913280/1‑100 (MQ=255)
                                                                    tGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                                 <  2:3019235/100‑1 (MQ=255)
                                                                    tGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCAcc                                 >  1:3019218/1‑100 (MQ=255)
                                                                          gCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATg                           <  2:1952773/100‑1 (MQ=255)
                                                                           cccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATgg                          <  1:2334611/100‑1 (MQ=255)
                                                                           cccGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATgg                          <  2:579496/100‑1 (MQ=255)
                                                                             cGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTg                        <  1:3214270/100‑1 (MQ=255)
                                                                              gTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGc                       >  1:2473765/1‑100 (MQ=255)
                                                                                  tGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTc                   <  2:1358779/100‑1 (MQ=255)
                                                                                  tGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTc                   <  2:185441/100‑1 (MQ=255)
                                                                                   gactttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTcc                  >  2:3459452/7‑100 (MQ=255)
                                                                                    ggTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCa                 <  1:226894/100‑1 (MQ=255)
                                                                                    ggTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCa                 <  2:4164926/100‑1 (MQ=255)
                                                                                     ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                >  2:4346822/5‑100 (MQ=255)
                                                                                     ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                >  2:1026318/5‑100 (MQ=255)
                                                                                     ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                >  2:3081697/5‑100 (MQ=255)
                                                                                     ctttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                >  2:1626332/5‑100 (MQ=255)
                                                                                     cttagCGACATCATCATGGCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATATGAAGACACAGGCAAAGTAACGTCCa                 >  1:107436/5‑99 (MQ=255)
                                                                                      tttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCat                >  2:2258901/4‑99 (MQ=255)
                                                                                       ttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatta              >  1:2772031/3‑100 (MQ=255)
                                                                                       ttgCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatta              >  1:935005/3‑100 (MQ=255)
                                                                                       ttgCGACATCATCATGTCGCAAAGTCACTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatta              >  2:2686465/3‑100 (MQ=255)
                                                                                       tGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCTAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatta              >  1:4386062/1‑100 (MQ=255)
                                                                                       tGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCatta              <  2:2497618/100‑1 (MQ=255)
                                                                                        ggCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCattat             >  1:3310633/1‑100 (MQ=255)
                                                                                         gCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCattatt            >  1:2175682/1‑100 (MQ=255)
                                                                                         gCGACATCATCATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCattatt            >  2:2659911/1‑100 (MQ=255)
                                                                                          cGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCAttatta           <  2:3760686/100‑1 (MQ=255)
                                                                                           gACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATtattat          >  2:2042460/1‑100 (MQ=255)
                                                                                           gACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCATCCATtattat          >  1:3889491/1‑100 (MQ=255)
                                                                                               tgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATcgcg      <  2:2774108/100‑1 (MQ=255)
                                                                                               tgatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATcgcg      <  2:83420/100‑1 (MQ=255)
                                                                                                gatgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCgg     >  1:1363087/1‑100 (MQ=255)
                                                                                                 atgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGa    >  2:3859807/1‑100 (MQ=255)
                                                                                                  tgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGAt   >  1:2233926/1‑100 (MQ=255)
                                                                                                  tgatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGCCCATTATTATCGCGGAt   >  1:1812479/1‑100 (MQ=255)
                                                                                                   gatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATg  <  2:3487100/100‑1 (MQ=255)
                                                                                                   gatgTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATACCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATtattatattggatg  >  1:2575309/1‑92 (MQ=255)
                                                                                                   |                                                                                                   
TGAAGGTGTAACGGAATACATGGCCTGGCTGAATCGCGACGCATAAGAGCTCTGCATGAAAATACTGGTGATCGGCCCGTCTTGGGTTGGCGACATGATGATGTCGCAAAGTCTCTATCGCACGCTCCAGGCGCGCTATCCCCAGGCGATAATCGATGTGATGGCACCGGCATGGTGCCGTCCATTATTATCGCGGATG  >  NC_000913/3794874‑3795072

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: