Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 14,116 A→G 17.9% intergenic (+37/‑52) dnaK → / → dnaJ chaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091314,1160AG17.9% 406.8 / 84.9 199intergenic (+37/‑52)dnaK/dnaJchaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone
Reads supporting (aligned to +/- strand):  ref base A (87/74);  new base G (18/17);  total (107/92)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

Only 100 of 206 total aligned reads displayed.

AACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCG  >  NC_000913/14017‑14215
                                                                                                   |                                                                                                   
  caacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAg                                                                                                   <  2:39792/100‑1 (MQ=255)
  caacGCGAAAGATGACGATGTTGCCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGaa                                                                                                    >  1:1731580/1‑99 (MQ=255)
   aacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAgg                                                                                                  <  1:2447266/100‑1 (MQ=255)
    acGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAggg                                                                                                 <  2:271684/100‑1 (MQ=255)
    acGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAggg                                                                                                 <  1:1668631/100‑1 (MQ=255)
     cgcgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAggg                                                                                                 <  1:4391723/99‑1 (MQ=255)
      gcgAAAGATGAAGATGTTGTCGAAGCTGAATTTGAAGAAGTAAAAGACAAAAAATAATCGCCCTATAAACGGAAAATTATACTGACACGGGCGAAGGGGa                                                                                               >  1:1514398/1‑100 (MQ=255)
       cgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGaa                                                                                              >  2:4231165/1‑100 (MQ=255)
       cgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGaa                                                                                              >  1:1414121/1‑100 (MQ=255)
       cgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGaa                                                                                              <  1:1583417/100‑1 (MQ=255)
          aaGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAttt                                                                                           >  2:94868/1‑100 (MQ=255)
            gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcn                                                                                         <  1:2189790/100‑2 (MQ=255)
            gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                         <  1:3609364/100‑1 (MQ=255)
            gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                         <  1:3160624/100‑1 (MQ=255)
             aTGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCct                                                                                        >  1:4166040/1‑100 (MQ=255)
              tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                         <  1:2712355/98‑1 (MQ=255)
              tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctc                                                                                       >  2:1887657/1‑100 (MQ=255)
                aCGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctct                                                                                      >  2:4056101/1‑99 (MQ=255)
                 cGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTcc                                                                                    >  2:3795059/1‑100 (MQ=255)
                    tgttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                 <  1:725769/100‑1 (MQ=38)
                     gttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCACTCCGcc                                                                                 >  1:3991414/1‑99 (MQ=255)
                       tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                              <  2:2896255/100‑1 (MQ=37)
                       tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                              <  2:3746913/100‑1 (MQ=37)
                       tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                              <  1:1074200/100‑1 (MQ=37)
                       tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGt                                                                              <  1:3583813/100‑1 (MQ=255)
                       tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACAGGCGAAGGGGAATTTCCTCTCAGCCCGt                                                                              >  1:1883833/1‑100 (MQ=255)
                       ngtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGt                                                                              >  1:1181412/2‑100 (MQ=255)
                        gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGTg                                                                             <  1:4231135/100‑1 (MQ=37)
                        gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                             <  2:3257721/100‑1 (MQ=255)
                        gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                             <  1:325368/100‑1 (MQ=255)
                              gCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                       <  2:4153548/100‑1 (MQ=255)
                              gCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                       <  1:2991950/100‑1 (MQ=255)
                               cTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCAt                                                                      >  2:1056445/1‑100 (MQ=255)
                                tGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATc                                                                     <  2:678520/100‑1 (MQ=255)
                                 gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCt                                                                    >  1:1328173/1‑100 (MQ=255)
                                  aaTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTa                                                                   >  2:3781312/1‑100 (MQ=255)
                                   aTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAg                                                                  >  1:3101922/1‑100 (MQ=255)
                                     ttGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAggg                                                                >  1:1935466/1‑100 (MQ=255)
                                      tGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGCCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAgggg                                                               >  1:63935/1‑100 (MQ=255)
                                       gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGc                                                              >  2:2236161/1‑100 (MQ=255)
                                           aagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAtt                                                          <  1:1022951/100‑1 (MQ=255)
                                            agTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAttt                                                         >  1:2071457/1‑100 (MQ=255)
                                             gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTa                                                        >  2:2438033/1‑100 (MQ=255)
                                             gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTa                                                        >  2:1425498/1‑100 (MQ=255)
                                              tCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaa                                                       <  2:1843552/100‑1 (MQ=255)
                                               cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                      <  2:4003244/100‑1 (MQ=255)
                                               cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                      <  2:3334414/100‑1 (MQ=255)
                                               cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCCCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                      <  1:1633827/100‑1 (MQ=255)
                                               cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGAGAAATTCCTCTCCGCCCGTGCATTCCTCTAGGGGCAATTTaaa                                                      >  1:2820867/1‑100 (MQ=255)
                                                 aaGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaa                                                    <  2:1400209/100‑1 (MQ=255)
                                                    aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                 <  2:2996696/100‑1 (MQ=255)
                                                    aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                 >  2:2409712/1‑100 (MQ=255)
                                                    aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                 >  2:3461312/1‑100 (MQ=255)
                                                     cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGCAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                >  1:265309/1‑100 (MQ=255)
                                                     cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                >  2:312886/1‑100 (MQ=255)
                                                     cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                >  2:2762227/1‑100 (MQ=255)
                                                     cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                >  1:3216666/1‑100 (MQ=255)
                                                            taatCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAg                                         <  1:3936202/100‑1 (MQ=255)
                                                              atCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCa                                       >  1:856093/1‑100 (MQ=255)
                                                              atCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCa                                       >  2:253885/1‑100 (MQ=255)
                                                               tCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCaa                                      <  2:3514365/100‑1 (MQ=255)
                                                                cGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAg                                     >  1:2359327/1‑100 (MQ=255)
                                                                cGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAg                                     >  1:384464/1‑100 (MQ=255)
                                                                  cccTATAAACGGGTAATTATACTGACACGGACGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                   >  1:2891708/1‑100 (MQ=255)
                                                                    cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCACTTTAAAAAAGATGGCTAAGCAAGatta                                 >  1:3791489/1‑100 (MQ=255)
                                                                    cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatta                                 <  2:3889684/100‑1 (MQ=255)
                                                                     tataAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGattat                                <  2:1259177/100‑1 (MQ=255)
                                                                      ataAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCACGattatt                               >  2:876224/1‑100 (MQ=255)
                                                                       taAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAttatta                              >  1:2188767/1‑100 (MQ=255)
                                                                         aaCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACg                            <  1:2206827/100‑1 (MQ=255)
                                                                          aCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGAAATTTAAAAAAGATGGCTAAGCGAGATTATTACga                           >  1:3192736/1‑100 (MQ=255)
                                                                             ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAt                        >  2:144403/1‑100 (MQ=255)
                                                                                aagtataCTGAGAGGGGCGAAGGGGAATGTCCTCTCCGGCGGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAtttt                     <  1:3689837/97‑1 (MQ=255)
                                                                                 aTTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                    <  2:3593828/100‑1 (MQ=255)
                                                                                  ttATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAg                   <  1:491535/100‑1 (MQ=255)
                                                                                   tataCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                  <  2:1782988/100‑1 (MQ=255)
                                                                                   tataCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                  >  1:338582/1‑100 (MQ=255)
                                                                                   tataCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                  <  1:2075883/100‑1 (MQ=255)
                                                                                      aCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGt               >  2:35792/1‑100 (MQ=255)
                                                                                       cTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGtt              >  2:536361/1‑100 (MQ=255)
                                                                                       cTGACACGGGCGAAGGGGAATTTCCCCTCCGCCCGTGCATGCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCg                >  1:1235785/1‑98 (MQ=255)
                                                                                       aatgcacGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCg                >  1:4381221/5‑98 (MQ=25)
                                                                                        tGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGttt             >  2:1376112/1‑100 (MQ=255)
                                                                                           cccGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTcc           <  2:3694938/97‑1 (MQ=255)
                                                                                           cacGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa          >  1:1844400/1‑100 (MQ=37)
                                                                                           cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa          >  2:119010/1‑100 (MQ=255)
                                                                                           cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa          <  2:4361305/100‑1 (MQ=255)
                                                                                           cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa          <  1:3719100/100‑1 (MQ=255)
                                                                                             cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa        <  1:2409701/100‑1 (MQ=255)
                                                                                             cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa        <  2:75647/100‑1 (MQ=255)
                                                                                              ggggGAAGGGGAATTGCCTCTCCGCCGGTGAATTCATCTAGGGGCAATTTAAAAAAGAGGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaa         <  1:969618/98‑1 (MQ=255)
                                                                                              gggCGGAGAGGACATTCCCCTTCGCCCGTGCAATCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTGTCCaaaa       >  1:1262985/1‑100 (MQ=37)
                                                                                              gggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa       >  2:2926349/1‑100 (MQ=255)
                                                                                               ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc      >  1:3476891/1‑100 (MQ=38)
                                                                                               ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc      >  2:1625029/1‑100 (MQ=38)
                                                                                               ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc      >  2:2142612/1‑100 (MQ=38)
                                                                                               ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc      >  2:3536120/1‑100 (MQ=255)
                                                                                               ggCGAAGGAGAATTTCCTCTCCGCCCGTGCATTCATCTAGAGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa        <  1:3551193/98‑1 (MQ=255)
                                                                                                gCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACa     <  1:1193409/100‑1 (MQ=255)
                                                                                                  gAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGc   <  2:3967594/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
AACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCG  >  NC_000913/14017‑14215

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: