Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 201,032 G→A 8.5% G21D (GGC→GAC)  lpxD → UDP‑3‑O‑(3‑hydroxymyristoyl)‑glucosamine N‑acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913201,0320GA 8.5% 574.0 / 33.9 186G21D (GGC→GAC) lpxDUDP‑3‑O‑(3‑hydroxymyristoyl)‑glucosamine N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (82/88);  new base A (8/8);  total (90/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.38e-03

Only 100 of 187 total aligned reads displayed.

TTGCTCAGAAAGCGCAGGCTTTTGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTACGCGCGCATGGCGCAAATTTTAGATACCACGCCGCAGCCCGCGCAGAACATTGCACCCAGTGCGGTGATCGACGCGACGGCGAAGC  >  NC_000913/200762‑201304
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
ttGCTCAGAAAGCGCAGGCTTTTGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTTATTTAGCCCAGCAGTTGGGTTCAGAACTACCCGGTGGTGGCg                                                                                                                                                                                                                                                                                >  4:533852/1‑273 (MQ=255)
ttGCTCAGAAAGCGCAGGCTTTTGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGat                                                                                                                                                                                                                                                                              <  2:216746/275‑1 (MQ=255)
  gCTCAGAAAGCGACGGCGTTTGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGata                                                                                                                                                                                                                                                                             <  4:72735/274‑1 (MQ=255)
      aGAAAGCGCAGGCTTTTGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGt                                                                                                                                                                                                                                                                         <  4:253442/274‑1 (MQ=255)
          agggaaGGCATTTGCGAAGGATCACGGACGAGGTCCCAACGGAGAACGGGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAAACCGTTGCGAACAGCCAGGATATCGATGTGGGTGTTGATGCAAACGCCGTTGCTTACAAAGGCAGCGATGTAAAAGACATCAATGCCGACGTACTGAAACAGGTTAAATAAGTAATGCAGTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtcat                                                                                                                                                                                                                                                                      <  4:470703/268‑1 (MQ=255)
           gcgcAGGCTTTTGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtcatca                                                                                                                                                                                                                                                                    <  2:541736/274‑1 (MQ=255)
            cgaaGGGTTTTGAGCAGGATCGCGCACGTAGTTCCAACGAAGAAAGCGGCAAACTGGTTACTCGTATCCAGACTGGTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCAcc                                                                                                                                                                                                                                                                  <  4:465380/272‑1 (MQ=255)
                  gttttGATCAGGATCGCGCACGTCGTTCAAACGAAGAACGCGGCAAACAGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGt                                                                                                                                                                                                                                                             <  2:416159/273‑1 (MQ=255)
                   ttttGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTg                                                                                                                                                                                                                                                           <  2:363540/275‑1 (MQ=255)
                             gATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTTCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGTTTTAGCGCATCAGTTGGATGCAGAACTCAACGGTGATGGCGATACCGTCATTCCAGGCGTTGCGGTCATgaa                                                                                                                                                                                                                                                 >  4:178521/1‑273 (MQ=255)
                                                                              tCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGAAGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGACTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCAc                                                                                                                                                                                                                                                                   >  1:73930/1‑208 (MQ=255)
                                                                               ccAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTg                                                                                                                                                                                               <  2:178522/275‑1 (MQ=255)
                                                                                   aCTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCCACAACAGGTCACATTACGTTCAGGGTTAACCCAAAATACCTTgccct                                                                                                                                                                                           >  4:480317/1‑271 (MQ=255)
                                                                                        tgtgAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAg                                                                                                                                                                                      >  4:521857/1‑275 (MQ=255)
                                                                                                              cAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGGAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAATAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGt                                                                                                                                                                <  2:92034/275‑1 (MQ=255)
                                                                                                                            ggTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAgg                                                                                                                                                  <  2:521858/275‑1 (MQ=255)
                                                                                                                                          aCGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTAGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGt                                                                                                                                                                                                                                                                         <  6:370443/142‑1 (MQ=255)
                                                                                                                                          aCGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTAGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGt                                                                                                                                                                                                                                                                         >  5:370443/1‑142 (MQ=255)
                                                                                                                                                  gCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTg                                                                                                                                                                                                                                                           <  5:901397/148‑1 (MQ=255)
                                                                                                                                                    ttACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGa                                                                                                                                                                                                                                                                               >  8:194779/1‑126 (MQ=255)
                                                                                                                                                       caacaGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGt                                                                                                                                                                                                                                                                         <  8:119685/129‑1 (MQ=255)
                                                                                                                                                       caacaGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGt                                                                                                                                                                                                                                                             >  8:742484/1‑141 (MQ=255)
                                                                                                                                                       caacaGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATACAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGt                                                                                                                                                                                                                                                             <  7:742484/141‑1 (MQ=255)
                                                                                                                                                          cagcagCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCAc                                                                                                                    <  2:160107/275‑1 (MQ=255)
                                                                                                                                                          cagcagCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCAc                                                                                                                    <  4:343590/275‑1 (MQ=255)
                                                                                                                                                                   tttAAAACACATCACGGGCGACTTACTGAAACAGGTAAAATAAGGAATGCCTGCAAGTCGACTGGCTAATTTAGGGCAGCAATCTGTGGCAGAACGACACGGTAAGGGCGAAAGCGTCATCACCGGAGTTGACTCCATGCAATCTGCACAAACCGGTCACATGACGTTCAGGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGTCGTCCGCGGGTGTCATGACCCAGGACGATCTTCCGTTCGCGAAAAGGGCCGCACTGGTAGTga                                                                                                           <  4:195654/273‑1 (MQ=255)
                                                                                                                                                                        aaGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCa                                                                                                                                                                                                                                                 >  5:651169/1‑136 (MQ=255)
                                                                                                                                                                        aaGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCg                                                                                                                                                                                                                                                                          <  8:1111042/111‑1 (MQ=255)
                                                                                                                                                                        aaGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCg                                                                                                                                                                                                                                                                          >  7:1111042/1‑111 (MQ=255)
                                                                                                                                                                        aaGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCGccatc                                                                                                                                                                                                                                                                     <  6:558722/116‑1 (MQ=255)
                                                                                                                                                                           aCATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACa                                                                                                                                                                                                                                    <  5:579222/146‑1 (MQ=255)
                                                                                                                                                                             aTCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTcac                                                                                                                                                                                                                              >  4:560782/1‑150 (MQ=255)
                                                                                                                                                                             aTCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGtt                                                                                                                                                               >  1:34748/1‑213 (MQ=255)
                                                                                                                                                                             aTCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTa                                                                                                 >  2:206753/1‑275 (MQ=255)
                                                                                                                                                                              tCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAg                                                                                                                                                                                                                                   >  6:657508/1‑144 (MQ=255)
                                                                                                                                                                                 cTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTa                                                                                                                                                                                                                          <  2:560783/150‑1 (MQ=255)
                                                                                                                                                                                 cTGCCGACGTACTANAACATGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTATCGCAGCAGTTGGATGCAGAACTACNCGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCa                                                                                                                                                           <  3:34748/213‑1 (MQ=255)
                                                                                                                                                                                      gACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTa                                                                                        >  4:388420/1‑275 (MQ=255)
                                                                                                                                                                                             tGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtcat                                                                                                                                                                                                                                                                      >  5:1095054/1‑94 (MQ=255)
                                                                                                                                                                                             tGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGtcat                                                                                                                                                                                                                                                                      <  6:1095054/94‑1 (MQ=255)
                                                                                                                                                                                             tGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGt                                                                                                                                                                                                                >  8:110496/1‑148 (MQ=255)
                                                                                                                                                                                               aaaCAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTa                                                                                                                                                                                                              <  7:110496/148‑1 (MQ=255)
                                                                                                                                                                                                aaCAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTaa                                                                                                                                                                                                             >  6:740363/1‑148 (MQ=255)
                                                                                                                                                                                                   aGGTTAAATAAGTAATGCCTTCAATTCGACTGGCGGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTACGCGCGCATGgcg                                                                           <  2:468046/275‑1 (MQ=255)
                                                                                                                                                                                                     gTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTACGCGCGTATGgcca                                                                          >  2:197488/1‑272 (MQ=255)
                                                                                                                                                                                                           tAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTg                                                                                                                                                                                                                                                           <  8:677876/91‑1 (MQ=255)
                                                                                                                                                                                                           tAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTg                                                                                                                                                                                                                                                           >  7:677876/1‑91 (MQ=255)
                                                                                                                                                                                                           tAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCGccatc                                                                                                                                                                                                                                                                     >  8:956914/1‑81 (MQ=255)
                                                                                                                                                                                                            aaGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCAt                                                                                                                                                                                                                   >  8:920459/1‑130 (MQ=255)
                                                                                                                                                                                                            aaGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCg                                                                                                                                                                                                 <  5:397097/148‑1 (MQ=255)
                                                                                                                                                                                                              gTAAGGGCTTCAATTCGACTGGCTGATTTAGCGAAGCAATTGGATGAAGTAATAAACGGTGAGGGCGATATCGTCAGCGCCGGCGTTGTGTCCATGCAATCTGAACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCAATTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAAGCCCTACCTGACTTACGCGCGGATGGCGCAAATTTTAGa                                                                <  4:195529/275‑1 (MQ=255)
                                                                                                                                                                                                               tAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTc                                                                                                                                                                      >  1:26961/1‑172 (MQ=255)
                                                                                                                                                                                                               tAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGAcc                                                                                                                                                      >  2:126353/1‑188 (MQ=255)
                                                                                                                                                                                                                aaTGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTcaca                                                                                                                                                                                                                             >  5:496935/1‑116 (MQ=255)
                                                                                                                                                                                                                 aTGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGc                                                                                                                                                                                                                                           >  7:856619/1‑101 (MQ=255)
                                                                                                                                                                                                                   gCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTg                                                                                                                                                                                               <  8:106925/143‑1 (MQ=255)
                                                                                                                                                                                                                   gCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTg                                                                                                                                                                                               >  7:106925/1‑143 (MQ=255)
                                                                                                                                                                                                                   gCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAgg                                                                                                                                                  <  4:126352/188‑1 (MQ=255)
                                                                                                                                                                                                                       tCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTg                                                                                                                                                                                                                                                           >  8:811991/1‑79 (MQ=255)
                                                                                                                                                                                                                       tCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTg                                                                                                                                                                                                                                                           <  7:811991/79‑1 (MQ=255)
                                                                                                                                                                                                                       tCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATg                                                                                                                                                                                                                                                   <  7:143518/87‑1 (MQ=255)
                                                                                                                                                                                                                       tCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTcaca                                                                                                                                                                                                                             >  5:659526/1‑109 (MQ=255)
                                                                                                                                                                                                                       tCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTaa                                                                                                                                                                                                             >  7:616411/1‑125 (MQ=255)
                                                                                                                                                                                                                       tCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTaa                                                                                                                                                                                                             <  8:616411/125‑1 (MQ=255)
                                                                                                                                                                                                                        cAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGAAGTGAAGAACCCCTACCTGACTTCCGCGCGCATGGCGCAAATTTTAGATACCAcgccg                                                      >  2:298224/1‑275 (MQ=255)
                                                                                                                                                                                                                         aaTTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTACCCCCGCCTGGCCCAAATTTTAGATACCAcgccc                                                      >  4:302888/1‑273 (MQ=255)
                                                                                                                                                                                                                         aaTTCGACTGGCTGATTTAGCGCAGCAGTTGGATCCAGAACTACACGGTGATGGCGATATCGCCATCACCGTCGTCGCGTCCATGCAATCTGCACAAACAGGTCACATTACATTCATGGTTAACTCAAAATACCGTGAGCATTTAGGTTTCTGCTAGGCGTTCGCGGTTGTCATGACCCAGTACGCACTTCCTTTCGCGATATGTGCCGCTATGCTAGTGAAGAATTTCTACCTTACTTCTGCCCTCATGGCGCCAATTTTaattatcacccccc                                                     >  4:136164/1‑262 (MQ=255)
                                                                                                                                                                                                                            tCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGGGATGACGATATCGccatc                                                                                                                                                                                                                                                                     <  7:946133/64‑1 (MQ=38)
                                                                                                                                                                                                                               aCTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAgg                                                                                                                                                                                                                                  <  7:310454/96‑1 (MQ=255)
                                                                                                                                                                                                                               aCTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCAt                                                                                                                                                                                          <  8:718353/136‑1 (MQ=255)
                                                                                                                                                                                                                                cTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCa                                                                                                                                                                                                                    >  5:784351/1‑109 (MQ=255)
                                                                                                                                                                                                                                 tGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTACGCGCGCATGGCGCAAATTTTAGATACCACGCCGCAGCCcgcg                                             <  2:133159/275‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACAt                                                                                                                                                                                                                            <  8:369570/99‑1 (MQ=255)
                                                                                                                                                                                                                                  ggCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCGCCATCACCgt                                                                                                                                                                                                                                                                <  6:154606/63‑2 (MQ=37)
                                                                                                                                                                                                                                  ggCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGACGATATCGCCATCACCgt                                                                                                                                                                                                                                                                >  5:154606/1‑62 (MQ=37)
                                                                                                                                                                                                                                    cTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTg                                                                                                                                                                                               >  7:377074/1‑126 (MQ=255)
                                                                                                                                                                                                                                    cTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCa                                                                                                                                                                            >  6:484055/1‑145 (MQ=255)
                                                                                                                                                                                                                                        tttAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTACGCGCGCATGGCGCAAATTTTAGATACAACGCCGCAGCCCGCGCAGaaaa                                      >  2:350267/1‑273 (MQ=255)
                                                                                                                                                                                                                                         ttAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCAt                                                                                                                                                                                                                                                    <  6:650942/68‑1 (MQ=255)
                                                                                                                                                                                                                                         ttAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCg                                                                                                                                                                                                 >  7:111333/1‑119 (MQ=255)
                                                                                                                                                                                                                                         ttAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCg                                                                                                                                                                                                 <  8:111333/119‑1 (MQ=255)
                                                                                                                                                                                                                                            gcgcAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTc                                                                                                                                                                                                                     >  6:566607/1‑96 (MQ=255)
                                                                                                                                                                                                                                               cagcagTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAATCAAGTCACATTACGTTCATGGTTAACCCAAAATa                                                                                                                                                                                                    <  7:419483/110‑1 (MQ=255)
                                                                                                                                                                                                                                               cagcagTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCaaa                                                                                                                                                                                                       >  7:900755/1‑107 (MQ=255)
                                                                                                                                                                                                                                               cagcagTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATa                                                                                                                                                                                                    >  8:419483/1‑110 (MQ=255)
                                                                                                                                                                                                                                                        gATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTcc                                                                                                                                                                     <  8:236578/132‑1 (MQ=255)
                                                                                                                                                                                                                                                          tGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATg                                                                                                                                                         >  6:483398/1‑142 (MQ=255)
                                                                                                                                                                                                                                                          tGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATg                                                                                                                                                         <  5:483398/142‑1 (MQ=255)
                                                                                                                                                                                                                                                           gCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCg                                                                                                                                                                        >  6:798660/1‑126 (MQ=255)
                                                                                                                                                                                                                                                           gCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCg                                                                                                                                                                        <  5:798660/126‑1 (MQ=255)
                                                                                                                                                                                                                                                              gAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCAt                                                                                                                                                                                                                   >  6:44501/1‑80 (MQ=255)
                                                                                                                                                                                                                                                              gAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAAc                                                                                                                                                                                                            <  7:994047/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCAt                                                                                                                                                                                                                   <  5:565900/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                aCTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCAt                                                                                                                                                                                                                   >  8:894515/1‑78 (MQ=255)
                                                                                                                                                                                                                                                                   acacGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGa                                                                                                                                              >  8:399449/1‑144 (MQ=255)
                                                                                                                                                                                                                                                                   acacGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGa                                                                                                                                              <  7:399449/144‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cacGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCAttt                                                                                                                                                                                        >  7:1067280/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                        gTGATGACGATATCGCCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGt                                                                                                                                                                >  7:226169/1‑121 (MQ=255)
                                                                                                                                                                                                                                                                         tGATGGCGTTATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACTTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTTTGCCAGGCGTCCGCGGTTGt                                                                                                                                                             <  8:7621/123‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGt                                                                                                                                                             >  7:7621/1‑123 (MQ=255)
                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                
TTGCTCAGAAAGCGCAGGCTTTTGAGCAGGATCGCGCACGTCGTTCCAACGAAGAACGCGGCAAACTGGTTACTCGTATCCAGACTGCTGTGAAATCCGTTGCCAACAGCCAGGATATCGATCTGGTTGTTGATGCAAACGCCGTTGCTTACAACAGCAGCGATGTAAAAGACATCACTGCCGACGTACTGAAACAGGTTAAATAAGTAATGCCTTCAATTCGACTGGCTGATTTAGCGCAGCAGTTGGATGCAGAACTACACGGTGATGGCGATATCGTCATCACCGGCGTTGCGTCCATGCAATCTGCACAAACAGGTCACATTACGTTCATGGTTAACCCAAAATACCGTGAGCATTTAGGCTTGTGCCAGGCGTCCGCGGTTGTCATGACCCAGGACGATCTTCCTTTCGCGAAAAGTGCCGCACTGGTAGTGAAGAATCCCTACCTGACTTACGCGCGCATGGCGCAAATTTTAGATACCACGCCGCAGCCCGCGCAGAACATTGCACCCAGTGCGGTGATCGACGCGACGGCGAAGC  >  NC_000913/200762‑201304

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