Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 700,407 C→A 66.0% 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,4070CA66.0% 154.3 / 251.0 236N336K (AAC→AAAnagEfused N‑acetyl glucosamine specific PTS enzyme: IIC, IIB, and IIA components
Reads supporting (aligned to +/- strand):  ref base C (33/47);  new base A (66/87);  total (99/137)
Fisher's exact test for biased strand distribution p-value = 8.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

Only 100 of 238 total aligned reads displayed.

CTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACT  >  CP009273/700258‑700554
                                                                                                                                                     |                                                                                                                                                   
cTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAc                                                                                                                                                     <  1:3259367/150‑1 (MQ=255)
  tGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAACCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGt                                                                                                                                                   <  2:4920461/150‑1 (MQ=255)
    cTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCt                                                                                                                                                 <  2:5209272/150‑1 (MQ=255)
     tCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTg                                                                                                                                                >  2:3486196/1‑150 (MQ=255)
     tCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTg                                                                                                                                                >  2:3521783/1‑150 (MQ=255)
     tCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTg                                                                                                                                                >  2:3781747/1‑150 (MQ=255)
       cGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGa                                                                                                                                              >  1:377701/1‑150 (MQ=255)
         ctgctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGAtg                                                                                                                                            >  2:4023621/1‑150 (MQ=255)
          tgctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGAtgc                                                                                                                                           >  2:1625670/1‑150 (MQ=255)
            ctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGAtgctg                                                                                                                                         <  1:5039608/150‑1 (MQ=255)
            ctgTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGAtgctg                                                                                                                                         >  1:3163123/1‑150 (MQ=255)
              gTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATgctgct                                                                                                                                       <  1:720863/150‑1 (MQ=255)
                     ctGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATg                                                                                                                                <  1:4325721/150‑1 (MQ=255)
                      tGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATgg                                                                                                                               <  2:4438630/150‑1 (MQ=255)
                      tGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATgg                                                                                                                               >  1:4764688/1‑150 (MQ=255)
                       gcacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATggg                                                                                                                              <  1:2935973/150‑1 (MQ=255)
                       gcacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATggg                                                                                                                              <  2:3221658/150‑1 (MQ=255)
                       gcacgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATggg                                                                                                                              >  1:3181666/1‑150 (MQ=255)
                         acgcacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCg                                                                                                                            >  2:4886154/1‑150 (MQ=255)
                            cacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTa                                                                                                                         >  1:3014614/1‑150 (MQ=255)
                            cacTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTa                                                                                                                         <  2:3622583/150‑1 (MQ=255)
                              ctgctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTAtc                                                                                                                       >  1:1635859/1‑150 (MQ=255)
                                gctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTAtctt                                                                                                                     <  1:4176823/150‑1 (MQ=255)
                                gctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTAtctt                                                                                                                     <  2:2802186/150‑1 (MQ=255)
                                gctgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTAtctt                                                                                                                     <  1:1758153/150‑1 (MQ=255)
                                  tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttct                                                                                                                   <  1:969980/150‑1 (MQ=255)
                                  tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTAtcttct                                                                                                                   >  2:2527633/1‑150 (MQ=255)
                                  tgACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTAtcttct                                                                                                                   <  1:262582/150‑1 (MQ=255)
                                   gACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATcttctt                                                                                                                  <  1:330627/150‑1 (MQ=255)
                                    aCCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCttcttc                                                                                                                 <  2:338672/150‑1 (MQ=255)
                                     ccGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCg                                                                                                                <  1:480240/150‑1 (MQ=255)
                                       ggTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCt                                                                                                              >  1:548025/1‑150 (MQ=255)
                                          aTCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATc                                                                                                           >  2:2746200/1‑150 (MQ=255)
                                          aTCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTGTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATc                                                                                                           <  2:1569781/150‑1 (MQ=255)
                                               ccTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACtt                                                                                                      <  2:136399/150‑1 (MQ=255)
                                                 tGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCg                                                                                                    >  2:2211567/1‑150 (MQ=255)
                                                     tgtgGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATATTCTTCGCTATCTACTTCgtggt                                                                                                <  1:3999079/150‑1 (MQ=255)
                                                       tgGCAACGCTGCTGGGTATCTACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGTTGGGCGTTATCTTCTTCGCTATCTACTTCGTGgtgt                                                                                              <  1:3389751/150‑1 (MQ=255)
                                                       tgGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGgtgt                                                                                              >  1:3830905/1‑150 (MQ=255)
                                                        ggCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGtt                                                                                             <  1:3059151/150‑1 (MQ=255)
                                                        ggCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCCGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGtt                                                                                             >  1:3230809/1‑150 (MQ=255)
                                                          cAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCa                                                                                           <  1:5283164/150‑1 (MQ=255)
                                                          cAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCa                                                                                           <  1:1033105/150‑1 (MQ=255)
                                                           aaCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          >  2:2733928/1‑150 (MQ=255)
                                                           aaCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          <  2:4500145/150‑1 (MQ=255)
                                                           aaCGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAg                                                                                          >  1:1094651/1‑150 (MQ=255)
                                                            aCGCTGCTGGGTATCCACGCGGGCTTCTCGTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGt                                                                                         <  2:1924672/150‑1 (MQ=255)
                                                              gctgctGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTc                                                                                       <  2:3865813/150‑1 (MQ=255)
                                                               ctgctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCg                                                                                      <  1:3672234/150‑1 (MQ=255)
                                                                tgctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCgg                                                                                     >  2:1944057/1‑150 (MQ=255)
                                                                 gctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGt                                                                                    <  1:5005559/150‑1 (MQ=255)
                                                                  ctgGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGtt                                                                                   >  1:2231493/1‑150 (MQ=255)
                                                                        aTCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGc                                                                             >  2:333013/1‑150 (MQ=255)
                                                                        aTCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGc                                                                             >  1:360183/1‑150 (MQ=255)
                                                                          ccACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCAt                                                                           <  2:230134/150‑1 (MQ=255)
                                                                            aCGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGt                                                                         <  2:3983578/150‑1 (MQ=255)
                                                                             cgcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGtt                                                                        >  2:2028255/1‑150 (MQ=255)
                                                                             cgcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGtt                                                                        <  2:1112065/150‑1 (MQ=255)
                                                                              gcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCGTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTc                                                                       <  2:3756857/150‑1 (MQ=255)
                                                                              gcgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTc                                                                       <  1:1600305/150‑1 (MQ=255)
                                                                               cgGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCa                                                                      >  2:4755547/1‑150 (MQ=255)
                                                                                   cgtctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                  <  1:4846726/148‑1 (MQ=255)
                                                                                   cTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                  <  1:5097556/150‑1 (MQ=255)
                                                                                   cTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                  <  2:4952872/150‑1 (MQ=255)
                                                                                   cTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                  <  1:1764545/150‑1 (MQ=255)
                                                                                   cTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCt                                                                  >  1:2706333/1‑150 (MQ=255)
                                                                                     tctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGa                                                                <  1:1019245/150‑1 (MQ=255)
                                                                                     tctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGa                                                                <  2:291270/150‑1 (MQ=255)
                                                                                      ctctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaa                                                               >  2:5275416/1‑150 (MQ=255)
                                                                                        ctttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGaaaa                                                             >  2:2095086/1‑150 (MQ=255)
                                                                                         tttctctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAAc                                                            <  2:3116246/150‑1 (MQ=255)
                                                                                             tctGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCg                                                        >  2:2652401/1‑150 (MQ=255)
                                                                                                 cGGGGGCTATCTACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGTTGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTc                                                    <  2:741699/150‑1 (MQ=255)
                                                                                                        tATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGa                                             >  1:2663720/1‑150 (MQ=255)
                                                                                                        tATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGa                                             >  2:3693905/1‑150 (MQ=255)
                                                                                                          tCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGTTGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATa                                           <  1:2184749/150‑1 (MQ=255)
                                                                                                           ctaCTACGCGTTGATGTATAACCTGCCGGCCGCCCGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGAATGTACAACCTGAAAACGCCGGGTCGTGAAGATaa                                          <  1:2446247/148‑1 (MQ=255)
                                                                                                           cGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaa                                          >  2:1683415/1‑150 (MQ=255)
                                                                                                            gACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaaa                                         >  1:5079069/1‑150 (MQ=255)
                                                                                                            gACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATaaa                                         >  1:4439380/1‑150 (MQ=255)
                                                                                                                ccgcgTTGATGTATAACCTGCCGGCCGCCAGCCAGAATGTCTGGATGCTGCTGGTGTTGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                     <  1:4266334/149‑1 (MQ=255)
                                                                                                                aCGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                     <  1:3419937/150‑1 (MQ=255)
                                                                                                                aCGCGTTGATGTATAACCTGCCGACCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAaagaag                                     <  2:1515812/150‑1 (MQ=255)
                                                                                                                  gcgTTGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGATATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTACAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGAc                                   <  1:2122827/150‑1 (MQ=255)
                                                                                                                   cgTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACg                                  <  2:4292815/150‑1 (MQ=255)
                                                                                                                     ttGATGTATAACCTGCCGGCCGCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTACAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgag                                <  1:1305866/150‑1 (MQ=255)
                                                                                                                      tGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACgaga                               >  2:5312304/1‑150 (MQ=255)
                                                                                                                       gATTTATACCCTGCCGGCCGCCCGCATGAACTTGTGTTTGCTGCTGGTTTTGGGCTTTATCTTCTTCGCTATCTACTTGGTGTTTTTCTGTTCGGTTATCCGAATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGAt                              <  1:2527642/150‑1 (MQ=255)
                                                                                                                       gATGTATACCCTGCCGGCCGCCAGCAATAACGTCTGGTGGCTGGTGGTGATGGGCGTTATCTTCTTCGCTATATACTTCGTGGTGTTCAGTTGGGTTATCCGCATGTACAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGAt                              <  1:3421579/150‑1 (MQ=255)
                                                                                                                        aTGTATACCCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCGGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCCGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATc                             <  2:2992829/150‑1 (MQ=255)
                                                                                                                              aaCCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACt                       <  2:1279180/150‑1 (MQ=255)
                                                                                                                               aCCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGTTGGGCTTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGTTTACCCGATTGTACAACCTTAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTg                      <  1:486805/150‑1 (MQ=255)
                                                                                                                                     cggccgCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGaagaag                <  2:4873498/150‑1 (MQ=255)
                                                                                                                                      ggccgCCAGCCAGAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGc               >  2:3200611/1‑150 (MQ=255)
                                                                                                                                           ccagccagAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACa          <  2:4827381/150‑1 (MQ=255)
                                                                                                                                               ccagAAAGTCTGGTTGCTGCTGGTGTTGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCTGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCaa      <  1:3781915/150‑1 (MQ=255)
                                                                                                                                               ccagAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATACGCATGTTCAACCGGAAAACGCCGGGTCGGGATGTTAAAGAAGACGATATAGTTACTGAAGAATCCAAAAGCaa      >  1:5051829/1‑150 (MQ=255)
                                                                                                                                                cagAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAac     >  2:2586151/1‑150 (MQ=255)
                                                                                                                                                cagAAAGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAac     >  2:4149394/1‑150 (MQ=255)
                                                                                                                                                  gAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAacac   >  2:5178466/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CTTGCTCCGCTGCTGTACCTCCTGCACGCACTGCTGACCGGTATCAGCCTGTTTGTGGCAACGCTGCTGGGTATCCACGCGGGCTTCTCTTTCTCTGCGGGGGCTATCGACTACGCGTTGATGTATAACCTGCCGGCCGCCAGCCAGAACGTCTGGATGCTGCTGGTGATGGGCGTTATCTTCTTCGCTATCTACTTCGTGGTGTTCAGTTCGGTTATCCGCATGTTCAACCTGAAAACGCCGGGTCGTGAAGATAAAGAAGACGAGATCGTTACTGAAGAAGCCAACAGCAACACT  >  CP009273/700258‑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: