Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,265,675 C→A 17.9% N454K (AAC→AAA dnaB → replicative DNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,265,6750CA17.9% 345.4 / 58.1 134N454K (AAC→AAAdnaBreplicative DNA helicase
Reads supporting (aligned to +/- strand):  ref base C (59/51);  new base A (12/12);  total (71/63)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 138 total aligned reads displayed.

TGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGG  >  NC_000913/4265576‑4265774
                                                                                                   |                                                                                                   
tgaCGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAAc                                                                                                     <  1:1596761/100‑1 (MQ=255)
    gtgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTc                                                                                                 <  2:2216102/100‑1 (MQ=255)
     tgtATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGAACGGTACGCCTGACCTTTAACGGTCa                                                                                                >  1:991998/1‑100 (MQ=255)
       tATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCagg                                                                                              <  2:2226860/100‑3 (MQ=255)
        aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATTGGGACGGTACGCCTGACCGTTAAAGGTCaggc                                                                                             <  2:2458452/100‑4 (MQ=255)
        aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggc                                                                                             <  2:1484162/100‑4 (MQ=255)
        aTCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggc                                                                                             <  2:1429741/100‑4 (MQ=255)
         tCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATgg                                                                                            >  2:1643048/1‑100 (MQ=255)
         tCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATgg                                                                                            <  2:2563776/100‑1 (MQ=255)
          cACGAAAACAGTGATTTAAAAGGCATCGCGGGAATTATTTTCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTCTAACGGTCAATGGt                                                                                           <  1:993929/100‑1 (MQ=255)
          cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGt                                                                                           >  2:2084302/1‑100 (MQ=255)
          cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGt                                                                                           >  2:326809/1‑100 (MQ=255)
          cACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggcgt                                                                                           <  1:850346/100‑6 (MQ=255)
           aCGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTc                                                                                          >  2:182407/1‑100 (MQ=255)
           aCGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTc                                                                                          >  2:1009280/1‑100 (MQ=255)
            cGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcg                                                                                         <  2:773892/100‑1 (MQ=255)
            cGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcg                                                                                         >  2:616101/1‑100 (MQ=255)
              aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg                                                                                       >  2:1149914/1‑100 (MQ=255)
              aaaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTcgcg                                                                                       >  2:1149854/1‑100 (MQ=255)
               aaaCAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCGTTAAAGGTCaggcgtaccgt                                                                                      <  2:1277175/100‑11 (MQ=38)
                  cAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGGAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTc                                                                                   >  1:2264472/1‑100 (MQ=255)
                   aGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCg                                                                                  <  1:668066/100‑1 (MQ=255)
                    gTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGa                                                                                 >  1:777664/1‑100 (MQ=255)
                     tGATGTAAAAGGCATCGCGGAAATTGTTATCGGTAAACAAGGTAACGGCGCAATCTGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCGTCGAc                                                                                <  1:741194/100‑1 (MQ=255)
                      gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAAGAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                               >  2:1824357/1‑100 (MQ=255)
                      gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                               >  1:870065/1‑100 (MQ=255)
                      gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                               >  2:581719/1‑100 (MQ=255)
                      gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                               <  1:2245416/100‑1 (MQ=255)
                      gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                               <  2:446791/100‑1 (MQ=255)
                      gATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACa                                                                               >  1:1522319/1‑100 (MQ=255)
                        tttAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTCACCTTTAACGGCCAATGGTCGCGCCTcaacaa                                                                              >  1:2153136/1‑99 (MQ=255)
                        tttAAAAGACATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAAc                                                                             <  2:2343788/100‑1 (MQ=255)
                          tAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTa                                                                           >  1:1942117/1‑100 (MQ=255)
                            aaaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATg                                                                         >  1:467270/1‑100 (MQ=255)
                             aaGGGATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGTGACGGTAGGGCTATCCTTTAACGGACAATGGTCGCGGTTCGACAACTa                                                                           >  1:644353/1‑97 (MQ=255)
                             aaGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGc                                                                        >  2:1223829/1‑100 (MQ=255)
                                gCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCggg                                                                     <  2:2444662/100‑1 (MQ=255)
                                 cATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCgggg                                                                    <  1:951407/100‑1 (MQ=255)
                                   tCGCGGAAATTATTATCGGTAAACAACTTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGcc                                                                  <  1:533664/100‑1 (MQ=255)
                                    cgcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCg                                                                 >  2:888553/1‑100 (MQ=255)
                                     gcgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGc                                                                >  1:2264243/1‑100 (MQ=255)
                                      cgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCa                                                               <  1:185084/100‑1 (MQ=255)
                                      cgGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCa                                                               <  1:2135082/100‑1 (MQ=255)
                                        gAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGt                                                             >  1:2006815/1‑100 (MQ=255)
                                         aaaTTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTa                                                            <  2:1662134/100‑1 (MQ=255)
                                             cttccgatctAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgac                                                        <  2:1698684/91‑1 (MQ=255)
                                                 atCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacgacg                                                    >  2:2371168/1‑100 (MQ=255)
                                                    ggTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGCCCTTTAACGGTCAATGGTCGCGCGTCGACAACGATGCGGGGCCGCAGTATGACGACGaat                                                 <  1:232756/100‑1 (MQ=255)
                                                    ggTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaat                                                 >  2:1140601/1‑100 (MQ=255)
                                                      tAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataa                                               <  2:1027200/100‑1 (MQ=255)
                                                       aaaCAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTacgacagatcggaag                                              >  1:1698675/1‑90 (MQ=255)
                                                       aaaCAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataat                                              >  2:2429139/1‑100 (MQ=255)
                                                        aacaacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataata                                             <  1:1141363/100‑1 (MQ=255)
                                                          caacGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGaataataat                                           <  2:1864553/100‑1 (MQ=255)
                                                           aacGTAACGGCCCAATCGGGACGGTGCGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAatt                                          <  1:1976922/100‑1 (MQ=255)
                                                              gTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAattatt                                       <  2:548561/100‑1 (MQ=255)
                                                              gTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAattatt                                       >  1:325351/1‑100 (MQ=255)
                                                               taccGGCCCAATCGGGACGGTACCCCTGTCCTTTAACGGTCAATGGTCGCGCTTAGAAAACTATGCGGGGCCGCAGTACGACGACGAATAATCattatt                                       <  1:248427/96‑1 (MQ=255)
                                                               tAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTAttt                                      <  2:1353443/100‑1 (MQ=255)
                                                                   ggCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATg                                  <  2:2483302/100‑1 (MQ=255)
                                                                    gCCCAATCGGGACGGTACGCCTGCCCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGa                                 >  1:671420/1‑100 (MQ=255)
                                                                    gCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGa                                 >  2:149343/1‑100 (MQ=255)
                                                                      ccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAt                               >  1:1399015/1‑100 (MQ=255)
                                                                      ccAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCACTCCGACGACGAATACTAATTATTTTATGaa                                >  1:1959453/1‑99 (MQ=255)
                                                                       cAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAAtt                              <  2:1788978/100‑1 (MQ=255)
                                                                        aaTCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGAAAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTa                             <  1:337546/100‑1 (MQ=255)
                                                                          tCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAgg                           <  2:1921686/100‑1 (MQ=255)
                                                                           cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                          <  2:1311979/100‑1 (MQ=255)
                                                                           cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                          <  2:2205209/100‑1 (MQ=255)
                                                                           cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                          >  2:2268317/1‑100 (MQ=255)
                                                                           cGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGt                          <  1:2413827/100‑1 (MQ=255)
                                                                            gggACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTa                         >  1:1484540/1‑100 (MQ=255)
                                                                                  cgaccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTaaa                   >  1:1913111/8‑100 (MQ=255)
                                                                                   gaccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAg                  >  2:2480193/7‑100 (MQ=255)
                                                                                   gaccattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAg                  <  1:1145291/94‑1 (MQ=255)
                                                                                    accattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc                 >  2:506908/6‑100 (MQ=255)
                                                                                    accattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc                 >  2:907603/6‑100 (MQ=255)
                                                                                    accattGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc                 >  2:1722039/6‑100 (MQ=255)
                                                                                    aCGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc                 <  2:1792151/100‑1 (MQ=255)
                                                                                    aCGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc                 <  2:1790624/100‑1 (MQ=255)
                                                                                    aCGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGc                 <  2:540636/100‑1 (MQ=255)
                                                                                     cGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCa                <  2:1900210/100‑1 (MQ=255)
                                                                                     cGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCa                <  2:1831553/100‑1 (MQ=255)
                                                                                      gCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaa               >  2:1972753/1‑100 (MQ=255)
                                                                                       ccTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACAACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaaa              >  1:241062/1‑100 (MQ=255)
                                                                                       attGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGCATAATCATTATTTTCTGAATTAGGTAATTAAAGCaaa              >  1:1744126/3‑100 (MQ=255)
                                                                                       attGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaaa              >  2:861150/3‑100 (MQ=255)
                                                                                       attGACCGTTAAAGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCaaa              >  2:1760654/3‑100 (MQ=255)
                                                                                        cTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAgatcgg             >  2:1970453/1‑95 (MQ=255)
                                                                                         tGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAaca            <  1:423197/100‑1 (MQ=255)
                                                                                          gACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAacac           >  1:851227/1‑100 (MQ=255)
                                                                                           aCCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACt          <  2:2379356/100‑1 (MQ=255)
                                                                                            ccTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACtt         >  1:841171/1‑100 (MQ=255)
                                                                                             cTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTa        >  1:569213/1‑100 (MQ=255)
                                                                                             cTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTa        >  2:761910/1‑100 (MQ=255)
                                                                                             cTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTa        >  2:1220137/1‑100 (MQ=255)
                                                                                              tttAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTAt       >  1:1310676/1‑100 (MQ=255)
                                                                                                   cGGTCAATTGTCGCGCTTCGACAACTATGTGGGGCCGCAGTACGACGACGAATAATAATTACTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAgg  <  1:1595379/100‑1 (MQ=255)
                                                                                                   cGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAgg  >  1:916693/1‑100 (MQ=255)
                                                                                                   cGGTCAATGGTCGCGCTTAGACAACTATGCGGGGACGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCac    <  1:513718/98‑2 (MQ=255)
                                                                                                   |                                                                                                   
TGAGGTGTATCACGAAAACAGTGATTTAAAAGGCATCGCGGAAATTATTATCGGTAAACAACGTAACGGCCCAATCGGGACGGTACGCCTGACCTTTAACGGTCAATGGTCGCGCTTCGACAACTATGCGGGGCCGCAGTACGACGACGAATAATAATTATTTTATGAATTAGGTAATTAAAGCAAACACTTATCAAGG  >  NC_000913/4265576‑4265774

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: