Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 14,124 T→A 17.4% intergenic (+45/‑44) 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,1240TA17.4% 464.3 / 83.4 195intergenic (+45/‑44)dnaK/dnaJchaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone
Reads supporting (aligned to +/- strand):  ref base T (86/74);  new base A (17/17);  total (104/91)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

Only 100 of 199 total aligned reads displayed.

AAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCG  >  NC_000913/14027‑14223
                                                                                                 |                                                                                                   
  gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                                 <  2:3676554/100‑1 (MQ=255)
  gATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                                 <  1:3160624/100‑1 (MQ=255)
    tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                                 <  1:2712355/98‑1 (MQ=255)
    tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctc                                                                                               >  1:3210391/1‑100 (MQ=255)
    tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctc                                                                                               >  1:2378774/1‑100 (MQ=255)
      aCGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctct                                                                                              >  2:4056101/1‑99 (MQ=255)
       cGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTcc                                                                                            >  2:3795059/1‑100 (MQ=255)
         aTGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGc                                                                                          <  1:4384025/100‑1 (MQ=38)
          tgttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                         <  1:725769/100‑1 (MQ=38)
          tgttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGcc                                                                                         >  2:4307633/1‑100 (MQ=255)
           gttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCACTCCGcc                                                                                         >  1:3991414/1‑99 (MQ=255)
            ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCg                                                                                       <  2:2449625/100‑1 (MQ=38)
            ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCg                                                                                       <  2:611137/100‑1 (MQ=255)
             tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                      <  2:3255686/100‑1 (MQ=37)
             tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                      <  2:3746913/100‑1 (MQ=37)
             tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACAGGCGAAGGGGAATTTCCTCTCAGCCCGt                                                                                      >  1:1883833/1‑100 (MQ=255)
              gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                                     <  2:2853130/100‑1 (MQ=255)
              gtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                                     >  2:3082821/1‑100 (MQ=255)
               tCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGTg                                                                                     <  1:3538000/99‑1 (MQ=37)
               tCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTg                                                                                     >  1:3576748/1‑99 (MQ=255)
                ccaCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCa                                                                                   >  2:2991966/3‑100 (MQ=255)
                cGACACTGAATTTTAAGAAGTCAAAGACAAAAAACAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGTAGGGGAATTTCCTCTCCGCCCGTgtc                                                                                   <  1:1185452/100‑3 (MQ=255)
                   cGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGTgtcag                                                                                 <  1:1223850/99‑5 (MQ=25)
                    gCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                               <  1:2991950/100‑1 (MQ=255)
                    gCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                               <  2:4153548/100‑1 (MQ=255)
                     cTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCAt                                                                              >  2:1056445/1‑100 (MQ=255)
                      tGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATc                                                                             <  1:3815155/100‑1 (MQ=255)
                       gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCt                                                                            >  1:1328173/1‑100 (MQ=255)
                        aaTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTa                                                                           >  2:3781312/1‑100 (MQ=255)
                        aaTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTa                                                                           <  2:450797/100‑1 (MQ=255)
                         aTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAg                                                                          >  1:3101922/1‑100 (MQ=255)
                           ttGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAggg                                                                        >  1:1935466/1‑100 (MQ=255)
                            tGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGCAGGGGAATGTCATCTCCGTCCGTACATCGATGTAgggg                                                                       >  1:3730806/1‑100 (MQ=255)
                            tGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGCCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAgggg                                                                       >  1:63935/1‑100 (MQ=255)
                             gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGc                                                                      <  2:1731586/100‑1 (MQ=255)
                             gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGc                                                                      >  2:1813315/1‑100 (MQ=255)
                             gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGc                                                                      >  2:2720107/1‑100 (MQ=255)
                             gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGc                                                                      >  2:2236161/1‑100 (MQ=255)
                              aagaagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCa                                                                     >  1:2258582/1‑100 (MQ=255)
                                 aagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAtt                                                                  <  1:1022951/100‑1 (MQ=255)
                                   gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTa                                                                >  2:1425498/1‑100 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                              >  1:1026219/1‑100 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                              <  2:915797/100‑1 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                              >  2:1857159/1‑100 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                              <  1:4108651/100‑1 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                              <  2:4003244/100‑1 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                              >  2:2916793/1‑100 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCCCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                              <  1:1633827/100‑1 (MQ=255)
                                     cAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGAGAAATTCCTCTCCGCCCGTGCATTCCTCTAGGGGCAATTTaaa                                                              >  1:2820867/1‑100 (MQ=255)
                                       aaGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaa                                                            <  2:1400209/100‑1 (MQ=255)
                                        aGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaa                                                            >  2:2075891/1‑99 (MQ=255)
                                          ataaaaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTACGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAg                                                          <  1:1436478/97‑1 (MQ=255)
                                          aCAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGa                                                         <  1:2091000/100‑1 (MQ=255)
                                           cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGCAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                        >  1:265309/1‑100 (MQ=255)
                                           cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAg                                                          >  1:2436712/1‑98 (MQ=255)
                                           cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                        >  2:953743/1‑100 (MQ=255)
                                           cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGAt                                                        >  1:3683393/1‑100 (MQ=255)
                                              aaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGc                                                     >  2:2426138/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)
                                                       gCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGa                                            >  1:3593809/1‑100 (MQ=255)
                                                        cccTATAAACGGGTAATTATACTGACACGGACGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGat                                           >  1:2891708/1‑100 (MQ=255)
                                                         ccTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGa                                            >  1:1696070/1‑98 (MQ=255)
                                                          cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatta                                         <  2:3889684/100‑1 (MQ=255)
                                                           tataAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGattat                                        <  2:1259177/100‑1 (MQ=255)
                                                             taAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAttatta                                      <  1:2191181/100‑1 (MQ=255)
                                                             taAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAttatta                                      >  1:2188856/1‑100 (MQ=255)
                                                                   ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAt                                >  2:144403/1‑100 (MQ=255)
                                                                        ttATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAg                           <  1:491535/100‑1 (MQ=255)
                                                                        ttATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAg                           <  1:3218095/100‑1 (MQ=255)
                                                                            aCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGt                       >  2:35792/1‑100 (MQ=255)
                                                                             cTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGtt                      >  2:3696958/1‑100 (MQ=255)
                                                                              tGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGttt                     >  2:1376112/1‑100 (MQ=255)
                                                                                 cacGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa                  >  1:1844400/1‑100 (MQ=37)
                                                                                 cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa                  >  2:119010/1‑100 (MQ=255)
                                                                                  acGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaa                 >  2:3330775/1‑100 (MQ=37)
                                                                                   cGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaa                <  2:75647/100‑1 (MQ=255)
                                                                                    ggggGAAGGGGAATTGCCTCTCCGCCGGTGAATTCATCTAGGGGCAATTTAAAAAAGAGGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaa                 <  1:969618/98‑1 (MQ=255)
                                                                                    gggCGGAGAGGACATTCCCCTTCGCCCGTGCAATCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTGTCCaaaa               >  1:1262985/1‑100 (MQ=37)
                                                                                    gggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCaaaa               >  2:2344016/1‑100 (MQ=37)
                                                                                     ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTACGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc              >  1:3309351/1‑100 (MQ=37)
                                                                                     ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc              >  1:347494/1‑100 (MQ=38)
                                                                                     ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc              >  1:2075629/1‑100 (MQ=38)
                                                                                     ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc              >  1:4190745/1‑100 (MQ=38)
                                                                                     ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc              >  2:3536120/1‑100 (MQ=255)
                                                                                        gAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGc           <  2:3967594/100‑1 (MQ=255)
                                                                                         aaGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCg          <  2:3343511/100‑1 (MQ=255)
                                                                                          aGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCgg         <  1:552155/100‑1 (MQ=255)
                                                                                          aGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCgg         <  2:2790764/100‑1 (MQ=255)
                                                                                           ggggAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGGTGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGa        <  1:1168811/100‑1 (MQ=255)
                                                                                           ggggAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGa        >  1:3719078/1‑100 (MQ=255)
                                                                                           ggggAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGa        >  1:1102678/1‑100 (MQ=255)
                                                                                             ggAATTTCCTCTCCGCCCGTTCATTCATCTAGGGGCAAATTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTATGCGTTTCCAAAACAGCGGaa       <  1:1556864/99‑1 (MQ=255)
                                                                                             ggAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAag      >  1:4265921/1‑100 (MQ=255)
                                                                                               aaTTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAagag    <  2:856097/100‑1 (MQ=255)
                                                                                                 tttCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCg  <  1:3692139/100‑1 (MQ=255)
                                                                                                 tttCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCg  <  1:3189029/100‑1 (MQ=255)
                                                                                                 tttCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCAGAAGAGCg  <  1:2062658/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
AAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCG  >  NC_000913/14027‑14223

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: