Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 825,545 T→G V188V (GTA→GTC ybhR ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913825,5450TG97.4% 650.8 / 4.5 189V188V (GTA→GTCybhRputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (76/108);  total (78/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 189 total aligned reads displayed.

CGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATTACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGCTCGCTGT  >  NC_000913/825448‑825642
                                                                                                 |                                                                                                 
cGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGAc                                                                                                 <  2:5418319/100‑1 (MQ=255)
   tGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGcc                                                                                              <  2:2841478/100‑1 (MQ=255)
   tGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGcc                                                                                              <  2:3464530/100‑1 (MQ=255)
   tGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACCATCATGACGcc                                                                                              <  2:3031083/100‑1 (MQ=255)
    gAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCg                                                                                             >  2:3151379/1‑100 (MQ=255)
      aGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGAt                                                                                           >  2:328602/1‑100 (MQ=255)
       gATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATa                                                                                          >  1:4330574/1‑100 (MQ=255)
       gATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATa                                                                                          >  2:499038/1‑100 (MQ=255)
        aTCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAg                                                                                         <  2:2339496/100‑1 (MQ=255)
         tCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgt                                                                                        >  2:201272/1‑100 (MQ=255)
           tGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtgg                                                                                      <  2:1408973/100‑1 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATTACGCCGATAgtggt                                                                                     >  2:3693867/1‑100 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggt                                                                                     >  1:2325273/1‑100 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggt                                                                                     <  2:4830886/100‑1 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggt                                                                                     >  1:3896679/1‑100 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggt                                                                                     <  1:5224722/100‑1 (MQ=255)
            gCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAgtggt                                                                                     >  1:2845315/1‑100 (MQ=255)
              cAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGa                                                                                   >  2:2375188/1‑100 (MQ=255)
              cAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGa                                                                                   <  2:2579730/100‑1 (MQ=255)
               aGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGAt                                                                                  >  2:1702760/1‑100 (MQ=255)
                ggtggtGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATc                                                                                 <  2:4645096/100‑1 (MQ=255)
                ggtggtGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATc                                                                                 <  1:789723/100‑1 (MQ=255)
                  tggtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCAt                                                                               >  1:580228/1‑100 (MQ=255)
                  tggtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCAt                                                                               <  1:5159624/100‑1 (MQ=255)
                   ggtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATg                                                                              <  1:2006593/100‑1 (MQ=255)
                    gtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATgg                                                                             >  1:3125922/1‑100 (MQ=255)
                    gtgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATgg                                                                             >  1:4091794/1‑100 (MQ=255)
                     tgAGCGGCGAAGCCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGAAGGAAAGTGAAGGGACGATCATGACGCCCATAGTGGTGATCATGGc                                                                            <  2:1466931/100‑1 (MQ=255)
                     tgAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGc                                                                            <  2:3940906/100‑1 (MQ=255)
                      gAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCg                                                                           <  1:1070685/100‑1 (MQ=255)
                      gAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGGGACGATCATGACGCCGATAGTGGTGATCATGGCg                                                                           <  1:4130895/100‑1 (MQ=255)
                        gcggcgAAACCTGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGag                                                                         >  2:5307052/1‑99 (MQ=255)
                        gcggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGAt                                                                         <  2:707663/100‑1 (MQ=255)
                         cggcgAAACCAGTAGCTGATCGAGTGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATc                                                                        <  2:1273233/100‑1 (MQ=255)
                         cggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATc                                                                        <  2:879862/100‑1 (MQ=255)
                         cggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATc                                                                        <  1:1902525/100‑1 (MQ=255)
                          ggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCa                                                                       >  2:222085/1‑100 (MQ=255)
                          ggcgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCa                                                                       >  2:4939860/1‑100 (MQ=255)
                            cgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGt                                                                     >  1:4262808/1‑100 (MQ=255)
                            cgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGt                                                                     <  1:5479966/100‑1 (MQ=255)
                            cgAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGt                                                                     <  1:5376476/100‑1 (MQ=255)
                               aaCCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATTAGTGAc                                                                  <  1:4644398/100‑1 (MQ=255)
                               aaCCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGAc                                                                  <  1:5032598/100‑1 (MQ=255)
                                 ccAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACgg                                                                >  2:4698308/1‑100 (MQ=255)
                                 ccAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACgg                                                                <  2:3408896/100‑1 (MQ=255)
                                  cAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGc                                                               >  1:3095869/1‑100 (MQ=255)
                                    gTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGAGAGTGGTGATCATGGCGATCAGTGACGGcac                                                             >  2:3810218/1‑100 (MQ=255)
                                     tAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcacc                                                            <  2:1868893/100‑1 (MQ=255)
                                     tAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcacc                                                            <  1:2609336/100‑1 (MQ=255)
                                       gCTGATCGAGCGTACGTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcaccac                                                          >  2:4170921/1‑100 (MQ=255)
                                       gCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcaccac                                                          >  2:4076360/1‑100 (MQ=255)
                                       gCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcaccac                                                          >  1:3523803/1‑100 (MQ=255)
                                       gCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGcaccac                                                          >  2:1748743/1‑100 (MQ=255)
                                         tGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACaa                                                        >  1:5014776/1‑100 (MQ=255)
                                           aTCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAc                                                      <  2:2449819/100‑1 (MQ=255)
                                            tCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAcc                                                     <  2:1413277/100‑1 (MQ=255)
                                            tCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAcc                                                     <  1:2439570/100‑1 (MQ=255)
                                              gAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCAt                                                   >  2:2679635/1‑100 (MQ=255)
                                              gAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCAt                                                   >  2:2685413/1‑100 (MQ=255)
                                               aGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCAtt                                                  <  2:145452/100‑1 (MQ=255)
                                                 cGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTg                                                <  1:1833945/100‑1 (MQ=255)
                                                 cGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTg                                                <  2:4040759/100‑1 (MQ=255)
                                                  gTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGt                                               <  2:3946286/100‑1 (MQ=255)
                                                  gTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGt                                               <  1:4109944/100‑1 (MQ=255)
                                                  gTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGt                                               <  1:3739807/100‑1 (MQ=255)
                                                   tACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTa                                              <  2:1168530/100‑1 (MQ=255)
                                                    aCCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAAGCATTTGTAg                                             >  2:4770669/1‑100 (MQ=255)
                                                    aCCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAg                                             <  1:4926479/100‑1 (MQ=255)
                                                        tGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCga                                         >  1:875897/1‑100 (MQ=255)
                                                              cGTTCGCGGGCGACGGAAAGTGAAGTGTCGATCATGACGCCGATAGTGGTGATCGTGGCCATCAGTGACAGCACCACAAACCATCTGTAGTCGAGATTCg                                   <  1:4839037/100‑1 (MQ=255)
                                                                ttCGCGGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCggg                                 <  1:614238/100‑1 (MQ=255)
                                                                  cgcgGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGtt                               <  2:948776/100‑1 (MQ=255)
                                                                  cgcgGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGtt                               <  1:3016698/100‑1 (MQ=255)
                                                                   gcgGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGATGGCACCACAAACCATTTGTAGTCGAGATTCGGGTta                              <  1:2167015/100‑1 (MQ=255)
                                                                    cgGGCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTtat                             <  2:1442814/100‑1 (MQ=255)
                                                                     gggCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTtata                            <  2:230050/100‑1 (MQ=255)
                                                                      ggCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATAc                           >  2:381584/1‑100 (MQ=255)
                                                                      ggCGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATAc                           <  1:1933159/100‑1 (MQ=255)
                                                                        ngACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCa                         >  1:137722/2‑100 (MQ=255)
                                                                        cGACGGAAAGTGAAGTGACTATCATGACGCCGATAGTGGTGAGCATGGCGATCAGTGACTGCACCACAAACCAGTTGTAGTCGAGATTCGGGTTATACCa                         >  2:797760/1‑100 (MQ=255)
                                                                        cGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCa                         >  2:2178158/1‑100 (MQ=255)
                                                                        cGACGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCa                         >  1:576951/1‑100 (MQ=255)
                                                                           cGGAAAGTGAAGTGCCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGtt                      <  1:257528/100‑1 (MQ=255)
                                                                           cGGAAAGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGtt                      <  2:1051466/100‑1 (MQ=255)
                                                                              aaaGTGAAGTGACGATCATTACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCg                   >  1:4343759/1‑100 (MQ=255)
                                                                              aaaGTGAAGTGACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCg                   <  2:1691657/100‑1 (MQ=255)
                                                                               aagtgaagtgACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGt                  >  1:3227020/1‑100 (MQ=255)
                                                                               aagtgaagtgACGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGt                  <  1:1117105/100‑1 (MQ=255)
                                                                                atctaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTa                 <  2:4297153/96‑1 (MQ=255)
                                                                                agtgaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTa                 >  1:2012944/1‑100 (MQ=255)
                                                                                agtgaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTa                 >  1:2012263/1‑100 (MQ=255)
                                                                                  tgaagtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTacc               >  2:43174/1‑100 (MQ=255)
                                                                                     agtgaCGATCATGACGTCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTaccacc            <  1:4590192/100‑1 (MQ=255)
                                                                                     agtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTaccacc            <  1:331720/100‑1 (MQ=255)
                                                                                     agtgaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTaccacc            >  1:1733808/1‑100 (MQ=255)
                                                                                      gtgaCGATCATGACGCCGATAGTGGTGGTCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTAccacca           <  2:306448/100‑1 (MQ=255)
                                                                                        gaCGATCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGc         >  1:1704529/1‑100 (MQ=255)
                                                                                         aCGACCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCATCAGCt        <  1:4999556/100‑1 (MQ=255)
                                                                                             tCATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGCTCGCt    <  1:5130224/100‑1 (MQ=255)
                                                                                              cATGACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGCTCGCtg   >  1:3675302/1‑100 (MQ=255)
                                                                                                 |                                                                                                 
CGATGAAGATCTGCCAGGTGGTGAGCGGCGAAACCAGTAGCTGATCGAGCGTACCTTGTTCACGTTCGCGGGCGACGGAAAGTGAAGTGACGATCATTACGCCGATAGTGGTGATCATGGCGATCAGTGACGGCACCACAAACCATTTGTAGTCGAGATTCGGGTTATACCAGTTGCGTACCACCAGCTCGCTGT  >  NC_000913/825448‑825642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: