Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 968,406 G→A 6.8% R14Q (CGG→CAG)  lpxK → lipid A 4'kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913968,4060GA 6.8% 267.1 / 4.2 85R14Q (CGG→CAG) lpxKlipid A 4'kinase
Reads supporting (aligned to +/- strand):  ref base G (42/37);  new base A (2/4);  total (44/41)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACA  >  NC_000913/968141‑968679
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                 
tACCACGGCGTTGAAAACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGGGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGACTACCATTGAAAAGGCAGACGAAATCTTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGAACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATtgc                                                                                                                                                                                                                                                                           <  2:381732/274‑1 (MQ=255)
                       acGAACGTGCGATTCAGGCGGCACGGGATGAGTTGCAAAAAAACCGTACATCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGcc                                                                                                                                                                                                                                                   <  2:381831/274‑1 (MQ=255)
                                                      gTTGCAGAAAAACCGTAACTCTCTGGTGATTGCCCACTGCTTGTCTACCATTGAAAAGGCAGACGAACTCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTCCTTGAGCACCTCGGCGTTTACGCGCACCTTAACAAAATGCAGTTTGGCCAATGCTCGAAAAAAGCTGGCCTGGTGAATCACCTTGGTGTCGGCTATTTCTCACACTCTCCTGGCTGTATGGCCTGGTTAGTGGCGCGTGCGGTCTTTGCtctca                                                                                                                                                                                                                    >  4:190874/1‑271 (MQ=255)
                                                               aaaCCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGCTCCGTCTTTGCTCTAAACTAAAACt                                                                                                                                                                                                           >  2:179903/1‑275 (MQ=255)
                                                                         tctctGGTTATTGCCCGCTCCTTGGCGACNTTGAAAAGGGCGGACGAATTCGTGGNGGTCGGGGATGGGGTCATTTTGGAAACCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACANAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGcc                                                                                                                                                                                                 <  3:97385/275‑1 (MQ=255)
                                                                                    tGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAAGGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTggcgg                                                                                                                                                                                                                                                                                  >  3:30329/1‑183 (MQ=255)
                                                                                        cACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTAt                                                                                                                                                                                                                                                                              <  1:30329/183‑1 (MQ=255)
                                                                                         aCCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCAATCCGTCTTTGCTATCAACTAAAACTGAAGCGCTCCTGGCGTTCCCCCGTAc                                                                                                                                                                                 >  2:256393/1‑275 (MQ=255)
                                                                                                        ttGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCATTTGTGGCGGCTATTGCTGCCAATCTGCGGGtt                                                                                                                                                                                                                                                           >  4:477157/1‑186 (MQ=255)
                                                                                                            aaaGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                                                                                                                                                       <  2:477162/186‑1 (MQ=255)
                                                                                                             aaGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGTCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGAGCCTGGCGTGCCCCAGTACCGGTTGTCGTGGTTGGTAAt                                                                                                                                                             >  1:1393/1‑275 (MQ=255)
                                                                                                                 gCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACCCGCCACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGCCCTGGTTAGTGGCGCGATCCCTCTTTGCTTAAACCTAAAACCTAAGCGCGCCTGGCACGCCCCCGTACCGGTTGCCGCTGTTGGTAAtcc                                                                                                                                                           >  4:449316/1‑273 (MQ=255)
                                                                                                                           tCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTAt                                                                                                                                                                                                                                                                              <  8:8747/149‑1 (MQ=255)
                                                                                                                                    tcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCAc                                                                                                                                                                                                                                                                     <  7:168113/149‑1 (MQ=255)
                                                                                                                                    tcgAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCAc                                                                                                                                                                                                                                                                     >  5:178445/1‑149 (MQ=255)
                                                                                                                                      gAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCActc                                                                                                                                                                                                                                                                   <  6:259045/149‑1 (MQ=255)
                                                                                                                                      gAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTCTGTGGCGGCTATTGCTGCCActc                                                                                                                                                                                                                                                                   <  7:13172/149‑1 (MQ=255)
                                                                                                                                      gAGGATGGTGTCATTGTGGAACGCGGTACACATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCActc                                                                                                                                                                                                                                                                   <  7:319901/149‑1 (MQ=255)
                                                                                                                                        ggATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaag                                                                                                                                                                                                                                                                 <  7:55716/149‑5 (MQ=255)
                                                                                                                                        ggATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaag                                                                                                                                                                                                                                                                 <  5:172768/149‑5 (MQ=255)
                                                                                                                                           tGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTg                                                                                                                                                                                                                                                              <  5:179194/149‑1 (MQ=255)
                                                                                                                                               gtCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTg                                                                                                                                                                                                                                                          >  7:46316/1‑149 (MQ=255)
                                                                                                                                                tCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGt                                                                                                                                                                                                                                                         <  8:153407/149‑1 (MQ=255)
                                                                                                                                                 cATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTa                                                                                                                                                                                                                                                        <  7:69767/149‑1 (MQ=255)
                                                                                                                                                  aTTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTAt                                                                                                                                                                                                                                                       <  7:88134/149‑1 (MQ=255)
                                                                                                                                                   ttGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATg                                                                                                                                                                                                                                                      >  8:52750/1‑149 (MQ=255)
                                                                                                                                                   ttGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATg                                                                                                                                                                                                                                                      >  8:297417/1‑149 (MQ=255)
                                                                                                                                                      tgGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaag                                                                                                                                                                                                                                                                 >  8:130520/1‑131 (MQ=255)
                                                                                                                                                      tgGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATAGCCGCCAcaaag                                                                                                                                                                                                                                                                 <  7:130520/135‑5 (MQ=255)
                                                                                                                                                       ggAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCt                                                                                                                                                                                                                                                  >  6:86917/1‑149 (MQ=255)
                                                                                                                                                        gAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTg                                                                                                                                                                                                                                                 >  5:33212/1‑149 (MQ=255)
                                                                                                                                                        gAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAATGGAAAAACCCCGGTCGTTGTctgg                                                                                                                   >  2:373457/1‑275 (MQ=255)
                                                                                                                                                        gAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTctgg                                                                                                                   >  2:358733/1‑275 (MQ=255)
                                                                                                                                                           cgcgGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCa                                                                                                                                                                                                                                                                      >  1:148090/1‑125 (MQ=255)
                                                                                                                                                            gcgGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCg                                                                                                                                                                                            <  2:173118/198‑1 (MQ=255)
                                                                                                                                                               gTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCActct                                                                                                                                                                                                                                                                  <  3:148090/125‑1 (MQ=255)
                                                                                                                                                                      tAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGAGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGAt                                                                                                                                                                                                                                   >  7:269059/1‑149 (MQ=255)
                                                                                                                                                                              tGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCttt                                                                                                                                                                                                                           <  6:178445/149‑1 (MQ=255)
                                                                                                                                                                               gCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTg                                                                                                                                                                                                                          <  8:46316/149‑1 (MQ=255)
                                                                                                                                                                                cTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGc                                                                                                                                                                                                                         <  6:33212/149‑1 (MQ=255)
                                                                                                                                                                                 ttGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                                                                                                                                                        <  5:420675/149‑1 (MQ=255)
                                                                                                                                                                                 ttGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCt                                                                                                                                                                                                                        <  7:297417/149‑1 (MQ=255)
                                                                                                                                                                                      cACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGTGCTGGGGACTCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTAtt                                                                                     <  2:431183/275‑1 (MQ=255)
                                                                                                                                                                                      cACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATaaa                                                                                                                                                                                                                   >  7:419994/1‑149 (MQ=255)
                                                                                                                                                                                         cgcgGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATcc                                                                                                                                                                                                                                 >  4:78738/1‑132 (MQ=255)
                                                                                                                                                                                         cgcgGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTa                                                                                                                                                                                                                <  5:216283/149‑1 (MQ=255)
                                                                                                                                                                                           cggcgTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGt                                                                                                                                                                                                                                                            >  2:514095/1‑103 (MQ=255)
                                                                                                                                                                                             gcgTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCt                                                                                                                                                                                                                             <  2:78739/132‑1 (MQ=255)
                                                                                                                                                                                               gTTTACGCGCAACTTCACAAAATGCTGTTTGGCCAATGATCAAAGAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTa                                                                                                                                                                                                                                                        <  4:514090/103‑1 (MQ=255)
                                                                                                                                                                                                tttACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGa                                                                                                                                                                                                         >  7:358584/1‑149 (MQ=255)
                                                                                                                                                                                                tttACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCCCCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGa                                                                                                           >  4:365536/1‑243 (MQ=255)
                                                                                                                                                                                                tttACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGa                                                                                                           <  2:365541/243‑1 (MQ=255)
                                                                                                                                                                                                  tACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAACAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCTCGATCCGTCTTTGCTATAAACTCAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCACCGGAAAAACCCCGTTCGTTGTCTGGCTGGTGGAACAGTTGCACCAGCCCGGTATTCGCGTGGGGGTc                                                                         >  2:10147/1‑275 (MQ=255)
                                                                                                                                                                                                     gcgcAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAgcgc                                                                                                                                                                                                    >  5:151165/1‑149 (MQ=255)
                                                                                                                                                                                                              cacaAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGaa                                                                                                                                                                                                        <  8:282406/136‑1 (MQ=255)
                                                                                                                                                                                                              cacaAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGaa                                                                                                                                                                                                        >  7:282406/1‑136 (MQ=255)
                                                                                                                                                                                                                 aaaaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtat                                                                                                                                                                                                                      >  8:112468/1‑119 (MQ=255)
                                                                                                                                                                                                                 aaaaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCtat                                                                                                                                                                                                                      <  7:112468/119‑1 (MQ=255)
                                                                                                                                                                                                                 aaaaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGatat                                                          <  2:375491/275‑1 (MQ=255)
                                                                                                                                                                                                                   aaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATagccgc                                                                                                                                                                                                                                                                        <  8:413707/67‑1 (MQ=255)
                                                                                                                                                                                                                   aaTGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCAGCAATagccgc                                                                                                                                                                                                                                                                        >  7:413707/1‑67 (MQ=255)
                                                                                                                                                                                                                      gCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGt                                                                                                                                                                                   >  5:387776/1‑149 (MQ=255)
                                                                                                                                                                                                                         gTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTAcc                                                                                                                                                                                >  7:85048/1‑149 (MQ=255)
                                                                                                                                                                                                                         gTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTAcc                                                                                                                                                                                >  6:32281/1‑149 (MQ=255)
                                                                                                                                                                                                                            ttgCCAATGTGTGAAAAAATCGGGTCTGGTGCATCCCCTTTTTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACGGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCGGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGc                                               <  4:358728/273‑1 (MQ=255)
                                                                                                                                                                                                                             ggCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCCGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGGATCCCGGGGATATGGTGGTAAGGCt                                              >  2:256001/1‑275 (MQ=255)
                                                                                                                                                                                                                               ccgaTGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTCCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGa                                            <  2:190876/272‑1 (MQ=255)
                                                                                                                                                                                                                                   tGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACtctc                                                                                                                                                                                                                                                                 >  2:257985/1‑58 (MQ=255)
                                                                                                                                                                                                                                     aTCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggt                                                                                                                                                                    >  5:397052/1‑149 (MQ=255)
                                                                                                                                                                                                                                       cGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTgg                                                                                                                                                                                                                                                             <  4:257981/58‑1 (MQ=255)
                                                                                                                                                                                                                                        gAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGGATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttgg                                                                                                                                                                 >  8:366429/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aaaaaaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttggt                                                                                                                                                                >  7:259296/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aaaaaaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGAGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGtggttggt                                                                                                                                                                <  6:189752/149‑1 (MQ=255)
                                                                                                                                                                                                                                           aaaaTCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTaa                                                                                                                                                              <  5:317124/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 tggtctggtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCAc                                                                                                                                                        <  5:32281/149‑1 (MQ=255)
                                                                                                                                                                                                                                                      tggtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAg                                                                                                                                                   >  6:193934/1‑149 (MQ=255)
                                                                                                                                                                                                                                                      tggtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAg                                                                                                                                                   >  6:457752/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        gtGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAggc                                                                                                                                                 >  8:458131/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             tccgcttgtaggcggcTATTCGTGCCACGCTCCTGGTTGTAAGGCCTGGTAAGTGGCGCGATCCGTCTTGTCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATAc              <  2:227906/265‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCaa                                                                                                                                            <  5:457752/149‑1 (MQ=255)
                                                                                                                                                                                                                                                             tCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCaa                                                                                                                                            >  8:367629/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                ccTTTGTGGCGTCTATTGCTTCCTCTCTCCTGGTTGTATGGCCGGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAAGCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACcac           <  2:235248/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGa                                                                                                                                        >  5:161380/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  tttGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCTGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAAATAccctaa         >  2:244254/1‑271 (MQ=255)
                                                                                                                                                                                                                                                                  tttGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGCTATGGTGGTAAGGCTGAATCTTAACCGCTGTTATTGTCGGCAGTTACCacaa         >  2:20079/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                     gtgGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGaaaaa                                                                                                                                    >  7:12743/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                      tgGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAAc                                                                                                                                   >  5:109445/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        gcggcTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAAccc                                                                                                                                 >  6:254243/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        gcggcTATCGGTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGGGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTAGCCGCTGTTATTGTCGGCAGATACCACAACAGca    <  2:217855/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cggcTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAAcccc                                                                                                                                <  7:367629/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         cggcTATTGCTGCCACTCTCCGGGTTGTATGGCCTGGTGATGGGCGCGACCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGGAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGcaca  <  2:13177/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                 
TACCTCGGCTCTGGATACCGAATCCGAACGTGCGATTCAGGCGGCACTGGATGAGTTGCAGAAAAACCGTACCTCTCTGGTGATTGCCCACCGCTTGTCTACCATTGAAAAGGCAGACGAAATCGTGGTCGTCGAGGATGGTGTCATTGTGGAACGCGGTACGCATAACGATTTGCTTGAGCACCGCGGCGTTTACGCGCAACTTCACAAAATGCAGTTTGGCCAATGATCGAAAAAATCTGGTCTGGTGAATCCCCTTTGTGGCGGCTATTGCTGCCACTCTCCTGGTTGTATGGCCTGGTGAGTGGCGCGATCCGTCTTTGCTATAAACTAAAACTGAAGCGCGCCTGGCGTGCCCCCGTACCGGTTGTCGTGGTTGGTAATCTCACCGCAGGCGGCAACGGAAAAACCCCGGTCGTTGTCTGGCTGGTGGAACAGTTGCAACAGCGCGGTATTCGCGTGGGGGTCGTATCGCGGGGATATGGTGGTAAGGCTGAATCTTATCCGCTGTTATTGTCGGCAGATACCACAACAGCACA  >  NC_000913/968141‑968679

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: