Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 700,407 C→A 8.2% N336K (AAC→AAA nagE → fused N‑acetyl glucosamine specific PTS enzyme: IIC, IIB, and IIA components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273700,4070CA 8.2% 791.2 / 28.8 235N336K (AAC→AAAnagEfused N‑acetyl glucosamine specific PTS enzyme: IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base C (85/125);  new base A (6/13);  total (91/144)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 235 total aligned reads displayed.

CTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACT  >  CP009273/700262‑700554
                                                                                                                                                 |                                                                                                                                                   
cTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTc                                                                                                                                                  <  1:2532004/149‑1 (MQ=255)
 tCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCt                                                                                                                                                 >  1:1833141/1‑149 (MQ=255)
 tCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCt                                                                                                                                                 <  1:1089410/149‑1 (MQ=255)
  ccGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTg                                                                                                                                                <  1:2313708/149‑1 (MQ=255)
       gctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgc                                                                                                                                           <  1:524090/149‑1 (MQ=255)
       gctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgc                                                                                                                                           <  2:3253670/149‑1 (MQ=255)
        ctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgct                                                                                                                                          <  1:2006314/149‑1 (MQ=255)
          gTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgctgc                                                                                                                                        <  1:470071/149‑1 (MQ=255)
           tACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATgctgct                                                                                                                                       <  2:3898319/149‑1 (MQ=255)
           tACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATgctgct                                                                                                                                       >  1:3513796/1‑149 (MQ=255)
            aCCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGctgctg                                                                                                                                      <  1:2413268/149‑1 (MQ=255)
              ctcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGt                                                                                                                                    <  2:2829374/149‑1 (MQ=255)
               tcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTg                                                                                                                                   <  1:1735227/149‑1 (MQ=255)
               tcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTg                                                                                                                                   <  1:3279770/149‑1 (MQ=255)
               tcctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTg                                                                                                                                   <  1:2518752/149‑1 (MQ=255)
                 ctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGAt                                                                                                                                 >  2:3954617/1‑149 (MQ=255)
                 ctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGAt                                                                                                                                 <  2:771235/149‑1 (MQ=255)
                  tGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATg                                                                                                                                <  2:741200/149‑1 (MQ=255)
                   gcacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATgg                                                                                                                               >  1:1894926/1‑149 (MQ=255)
                   gcacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTGTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATgg                                                                                                                               <  2:2012624/149‑1 (MQ=255)
                     acgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGc                                                                                                                             >  2:3796417/1‑149 (MQ=255)
                       gcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGt                                                                                                                           <  1:1030759/149‑1 (MQ=255)
                       gcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGt                                                                                                                           >  1:2053844/1‑149 (MQ=255)
                         acTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTa                                                                                                                         >  1:2417352/1‑149 (MQ=255)
                         acTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTa                                                                                                                         >  1:346943/1‑149 (MQ=255)
                           tgctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtc                                                                                                                       >  2:931095/1‑149 (MQ=255)
                            gctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtct                                                                                                                      >  1:3839343/1‑149 (MQ=255)
                             ctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtctt                                                                                                                     <  1:2699106/149‑1 (MQ=255)
                              tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttc                                                                                                                    <  1:2541684/149‑1 (MQ=255)
                               gACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttct                                                                                                                   <  1:3941144/149‑1 (MQ=255)
                               gACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttct                                                                                                                   >  2:3055583/1‑149 (MQ=255)
                                 ccGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCttcttc                                                                                                                 >  2:2977500/1‑149 (MQ=255)
                                   ggTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGc                                                                                                               <  1:2165619/149‑1 (MQ=255)
                                    gTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCt                                                                                                              >  2:2523276/1‑149 (MQ=255)
                                    gTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTGTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCt                                                                                                              <  2:3684162/149‑1 (MQ=255)
                                        cAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCt                                                                                                          <  2:1720055/149‑1 (MQ=255)
                                         aGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTa                                                                                                         <  1:2989647/149‑1 (MQ=255)
                                          gCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTAc                                                                                                        >  1:1532620/1‑149 (MQ=255)
                                            cTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACtt                                                                                                      >  1:553634/1‑149 (MQ=255)
                                              gTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCg                                                                                                    >  2:3564013/1‑149 (MQ=255)
                                               tttGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCgt                                                                                                   <  1:365278/149‑1 (MQ=255)
                                                ttGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCgtg                                                                                                  <  1:3540374/149‑1 (MQ=255)
                                                ttGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCGGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCgtg                                                                                                  <  2:1553689/149‑1 (MQ=255)
                                                    ggCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGgtgt                                                                                              <  1:905834/149‑1 (MQ=255)
                                                    ggCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGgtgt                                                                                              >  2:592008/1‑149 (MQ=255)
                                                     gCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGtt                                                                                             >  1:911945/1‑149 (MQ=255)
                                                        aCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          >  2:1829325/1‑149 (MQ=255)
                                                        aCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          >  2:1829294/1‑149 (MQ=255)
                                                          gctgctGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGtt                                                                                        <  1:3325475/149‑1 (MQ=255)
                                                             gctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCgg                                                                                     >  1:760640/1‑149 (MQ=255)
                                                               tgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGtt                                                                                   >  2:1622033/1‑149 (MQ=255)
                                                               tgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGGGGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCgggt                                                                                   >  1:2838905/1‑147 (MQ=255)
                                                                gggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                  <  1:3655086/149‑1 (MQ=255)
                                                                gggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                  <  2:2248734/149‑1 (MQ=255)
                                                                gggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                  >  1:172185/1‑149 (MQ=255)
                                                                gggTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTa                                                                                  <  1:1587436/149‑1 (MQ=255)
                                                                   tATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATcc                                                                               >  2:1471039/1‑149 (MQ=255)
                                                                    aTCCACGCGGGCTTCTCGTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCg                                                                              <  2:716800/149‑1 (MQ=255)
                                                                         cgcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGt                                                                         <  2:3513796/149‑1 (MQ=255)
                                                                            gggCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCa                                                                      <  2:3191846/149‑1 (MQ=255)
                                                                              gCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAc                                                                    <  2:129553/149‑1 (MQ=255)
                                                                              gCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAAc                                                                    <  2:128926/149‑1 (MQ=255)
                                                                                ttctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                  <  1:2598265/149‑1 (MQ=255)
                                                                                ttctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                  <  1:530217/149‑1 (MQ=255)
                                                                                   tctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaa                                                               >  1:601422/1‑149 (MQ=255)
                                                                                     tttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                             <  2:2616411/149‑1 (MQ=255)
                                                                                       tctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACg                                                           >  1:2898597/1‑149 (MQ=255)
                                                                                        ctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGc                                                          <  2:2548989/149‑1 (MQ=255)
                                                                                           tGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCgg                                                       >  2:497324/1‑149 (MQ=255)
                                                                                              gggggCTATCGACTACGCGTTTATTTATAACGTGCCGGCCGCCAGCCAGAACGTCTGGTTGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCTGTTGGGTTATCCGCTTTTTCAACCTGAAAACGCCGGGTc                                                    <  2:411879/149‑1 (MQ=255)
                                                                                              gggggCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTc                                                    <  2:601414/149‑1 (MQ=255)
                                                                                              gggggCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTc                                                    >  1:3048802/1‑149 (MQ=255)
                                                                                                  gCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGa                                                >  2:268749/1‑149 (MQ=255)
                                                                                                     aTCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTGTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGa                                             <  2:98505/149‑1 (MQ=255)
                                                                                                      tCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGAt                                            >  2:272730/1‑149 (MQ=255)
                                                                                                       cGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATa                                           >  1:1578991/1‑149 (MQ=255)
                                                                                                         aCTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCTGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaaa                                         <  2:526741/149‑1 (MQ=255)
                                                                                                           tACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaaga                                       <  1:4105561/149‑1 (MQ=255)
                                                                                                           tACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaaga                                       <  1:4104751/149‑1 (MQ=255)
                                                                                                            aCGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCTTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTACAACCTGAAAACGCCGGGTCGTGAAGATAaagaa                                      <  1:639595/149‑1 (MQ=255)
                                                                                                            aCGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaa                                      <  2:112433/149‑1 (MQ=255)
                                                                                                                gTTGATGTATAACCGGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGTTGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACg                                  <  1:3425193/149‑1 (MQ=255)
                                                                                                                 ttGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACga                                 <  2:2467693/149‑1 (MQ=255)
                                                                                                                 ttGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACga                                 <  2:4023070/149‑1 (MQ=255)
                                                                                                                   gATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgaga                               <  2:3975160/149‑1 (MQ=255)
                                                                                                                   gATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgaga                               >  2:2237050/1‑149 (MQ=255)
                                                                                                                          aaCCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTAc                        >  1:1666436/1‑149 (MQ=255)
                                                                                                                             cTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTga                     >  2:1418066/1‑149 (MQ=255)
                                                                                                                             cTGCCGGCCGCCAGCCAGAACGTCTGGAGGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGACGATAAAGAAGACGAGATCGTTACTga                     >  2:3289306/1‑149 (MQ=255)
                                                                                                                                ccggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaaga                  <  2:1894926/149‑1 (MQ=255)
                                                                                                                                 cggccgCCAGCCATAATGTTTGGTTGCTGCTGGTGATGGGCTTTATCTTCTTCGCTATCTACTTCGTGGTTTTCATTTCGGTTATCCGAATGTTCAACCTTAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaagaa                 <  1:3579991/149‑1 (MQ=255)
                                                                                                                                 cggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTgaagaa                 <  2:2511802/149‑1 (MQ=255)
                                                                                                                                  ggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGaagaag                >  1:2574626/1‑149 (MQ=255)
                                                                                                                                        cagccagAAGGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACa          <  1:3263951/149‑1 (MQ=255)
                                                                                                                                        cagcaatAACGTGTGGATGCTGGTGGTGTTGGGGGTTATCTTGTTCGATATCTACTTCGTGGTTTTCAGTTGGGTTATCCGATTGTTCAACCTCAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACa          <  1:1120208/149‑1 (MQ=255)
                                                                                                                                        cagcaagAACGTCTGGATGCTGCGGGTGTTGGGCTTTATGTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGAATGTACAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACa          <  1:1563210/149‑1 (MQ=255)
                                                                                                                                          gccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGc        >  2:277064/1‑149 (MQ=255)
                                                                                                                                          gccagAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGc        >  1:851731/1‑149 (MQ=255)
                                                                                                                                               aaCGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGTTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAacac   <  2:3616749/149‑1 (MQ=255)
                                                                                                                                                aCGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACt  >  2:629054/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACT  >  CP009273/700262‑700554

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