Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 274,735 A→T 90.5% pseudogene (138/348 nt) afuB ← pseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273274,7350AT90.5% 1064.1 / 78.1 325pseudogene (138/348 nt)afuBpseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (20/11);  new base T (169/125);  total (189/136)
Fisher's exact test for biased strand distribution p-value = 5.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

Only 100 of 329 total aligned reads displayed.

ATGAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTC  >  CP009273/274586‑274884
                                                                                                                                                     |                                                                                                                                                     
   aaaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATg                                                                                                                                                    >  2:3214465/1‑150 (MQ=255)
    aaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATgg                                                                                                                                                   >  2:4173779/1‑150 (MQ=255)
      ataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcg                                                                                                                                                 >  2:1700970/1‑150 (MQ=255)
        aatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgcg                                                                                                                                               >  2:858271/1‑150 (MQ=255)
        aatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgcg                                                                                                                                               >  2:2111108/1‑150 (MQ=255)
          tCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGca                                                                                                                                             >  1:2095657/1‑150 (MQ=255)
          tCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGca                                                                                                                                             >  1:3516689/1‑150 (MQ=255)
            gCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGcaca                                                                                                                                           <  1:1421721/150‑1 (MQ=255)
            gCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGcaca                                                                                                                                           >  2:61614/1‑150 (MQ=255)
            gCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGcaca                                                                                                                                           <  2:2811353/150‑1 (MQ=255)
               aGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAg                                                                                                                                        >  2:1928091/1‑150 (MQ=255)
               aGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAg                                                                                                                                        <  1:1706504/150‑1 (MQ=255)
                 cATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCt                                                                                                                                      <  1:395397/150‑1 (MQ=255)
                  aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTg                                                                                                                                     <  2:1573854/150‑1 (MQ=255)
                  aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTg                                                                                                                                     >  2:4308844/1‑150 (MQ=255)
                   tCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGt                                                                                                                                    >  2:2018794/1‑150 (MQ=255)
                   tCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGt                                                                                                                                    >  2:1538431/1‑150 (MQ=255)
                    caccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTa                                                                                                                                   >  1:703505/1‑150 (MQ=255)
                    caccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTa                                                                                                                                   <  2:3282799/150‑1 (MQ=255)
                    caccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTa                                                                                                                                   >  1:3557610/1‑150 (MQ=255)
                     accacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTaa                                                                                                                                  <  2:3010575/150‑1 (MQ=255)
                     accacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTaa                                                                                                                                  >  1:3725864/1‑150 (MQ=255)
                       cacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAAt                                                                                                                                >  2:2594431/1‑150 (MQ=255)
                       cacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAAt                                                                                                                                >  1:859628/1‑150 (MQ=255)
                          gATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAg                                                                                                                             >  1:1262097/1‑150 (MQ=255)
                            tCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCCCAAAGCTGTAAATCAgcg                                                                                                                           >  1:3219175/1‑150 (MQ=255)
                            tCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcg                                                                                                                           >  1:757028/1‑150 (MQ=255)
                            tCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcg                                                                                                                           <  2:68189/150‑1 (MQ=255)
                             cAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcgc                                                                                                                          >  2:2539411/1‑150 (MQ=255)
                             cAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcgc                                                                                                                          >  2:1042783/1‑150 (MQ=255)
                             cAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcgc                                                                                                                          >  2:4144922/1‑150 (MQ=255)
                                aaTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCga                                                                                                                       >  1:3419802/1‑150 (MQ=255)
                                aaTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCga                                                                                                                       >  1:922203/1‑150 (MQ=255)
                                 aTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCgag                                                                                                                      <  1:2827544/150‑1 (MQ=255)
                                      aCCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGAt                                                                                                                 >  1:1169172/1‑150 (MQ=255)
                                          tAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGcc                                                                                                             >  1:1945817/1‑150 (MQ=255)
                                          tAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGcc                                                                                                             <  2:1425094/150‑1 (MQ=255)
                                            ggCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCgg                                                                                                           >  1:240591/1‑150 (MQ=255)
                                              cAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGgc                                                                                                         >  2:3137820/1‑150 (MQ=255)
                                                aTCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGgcgc                                                                                                       >  2:3093741/1‑150 (MQ=255)
                                                  cGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcag                                                                                                     >  2:3944891/1‑150 (MQ=255)
                                                   gCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagc                                                                                                    >  1:3895026/1‑150 (MQ=255)
                                                    cTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagca                                                                                                   >  1:4192347/1‑150 (MQ=255)
                                                     tACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGcagcag                                                                                                  >  1:4501614/1‑150 (MQ=255)
                                                      acacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCagcagc                                                                                                 <  1:4129554/150‑1 (MQ=255)
                                                      acacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCagcagc                                                                                                 >  1:3195730/1‑150 (MQ=255)
                                                        acCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCgg                                                                                               <  1:3106461/150‑1 (MQ=255)
                                                          cGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGAGGGTATCAGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGTCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgca                                                                                             >  1:2057147/1‑150 (MQ=255)
                                                            tATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcagc                                                                                           >  1:595393/1‑150 (MQ=255)
                                                                cGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGa                                                                                       <  1:2631516/150‑1 (MQ=255)
                                                                cGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAAGGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGa                                                                                       <  2:3631542/150‑1 (MQ=255)
                                                                  ccATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtg                                                                                     >  2:1723300/1‑150 (MQ=255)
                                                                      cTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCg                                                                                 >  1:1849002/1‑150 (MQ=255)
                                                                        tCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTa                                                                               <  1:4276543/150‑1 (MQ=255)
                                                                        tCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTa                                                                               >  2:3866647/1‑150 (MQ=255)
                                                                          cACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAACGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAAt                                                                             <  2:2095657/150‑1 (MQ=255)
                                                                          cACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAAt                                                                             <  2:4030434/150‑1 (MQ=255)
                                                                          cACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAAt                                                                             >  2:2887164/1‑150 (MQ=255)
                                                                            cgcgGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATgg                                                                           <  2:2380551/150‑1 (MQ=255)
                                                                              cgGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTa                                                                         >  2:4652940/1‑150 (MQ=255)
                                                                              cgGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTa                                                                         >  2:1620514/1‑150 (MQ=255)
                                                                               ggTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTAc                                                                        <  2:1167499/150‑1 (MQ=255)
                                                                               ggTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTAc                                                                        >  1:1079404/1‑150 (MQ=255)
                                                                                    aGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAgg                                                                   >  2:4391639/1‑150 (MQ=255)
                                                                                       aTGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGaa                                                                <  1:1415247/150‑1 (MQ=255)
                                                                                          tAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCg                                                             >  2:2586947/1‑150 (MQ=255)
                                                                                           aaGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCgg                                                            >  2:263863/1‑150 (MQ=255)
                                                                                            aGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgc                                                           >  1:2355932/1‑150 (MQ=255)
                                                                                                gTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgcgcgc                                                       >  1:2906997/1‑150 (MQ=255)
                                                                                                 tCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCa                                                      >  2:1266110/1‑150 (MQ=255)
                                                                                                    cTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGc                                                   >  1:1484811/1‑150 (MQ=255)
                                                                                                       cgcgTGTATCCGGCGTGACGAGGAACACAATGGCGCTGCCGGTGGTTATGGCGCGCACAAAGATGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtga                                                <  1:1667177/150‑1 (MQ=255)
                                                                                                       cgcgGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtga                                                <  1:4071997/150‑1 (MQ=255)
                                                                                                         cgGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagt                                              <  2:2503678/150‑1 (MQ=255)
                                                                                                          gggTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagtg                                             >  2:1755071/1‑150 (MQ=255)
                                                                                                              aTCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGc                                         >  2:4460208/1‑150 (MQ=255)
                                                                                                                ccGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCtt                                       <  1:2511930/150‑1 (MQ=255)
                                                                                                                 cGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTc                                      >  2:4630495/1‑150 (MQ=255)
                                                                                                                  ggCGTGATGAGGAACACATTGGCGATGACGGTGGTTATGGCGCGCACTAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGTATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCa                                     <  1:2159710/150‑1 (MQ=255)
                                                                                                                  ggCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCa                                     <  2:4790933/150‑1 (MQ=255)
                                                                                                                    cGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATc                                   >  1:4771210/1‑150 (MQ=255)
                                                                                                                    cGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATc                                   <  1:4682950/150‑1 (MQ=255)
                                                                                                                      tGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGa                                 <  2:1849309/150‑1 (MQ=255)
                                                                                                                        aCGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGc                               >  1:3591903/1‑150 (MQ=255)
                                                                                                                         cGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCg                              >  1:1009031/1‑150 (MQ=255)
                                                                                                                         cGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCg                              >  2:2672960/1‑150 (MQ=255)
                                                                                                                         cGAGGAACACAATGGCGCTGACGGGGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCg                              >  1:1833256/1‑150 (MQ=255)
                                                                                                                          gAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCTCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGa                             >  1:3973507/1‑150 (MQ=255)
                                                                                                                          gAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGa                             >  2:3829422/1‑150 (MQ=255)
                                                                                                                          gAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGa                             <  1:1928119/150‑1 (MQ=255)
                                                                                                                             gAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGAttt                          <  2:1945116/150‑1 (MQ=255)
                                                                                                                              aaCACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTa                         >  1:4570828/1‑150 (MQ=255)
                                                                                                                              aaCACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTa                         <  1:472231/150‑1 (MQ=255)
                                                                                                                              aaCACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTa                         >  2:2270142/1‑150 (MQ=255)
                                                                                                                                cacaATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATc                       <  2:3416779/150‑1 (MQ=255)
                                                                                                                                   aaTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAt                    >  1:2291105/1‑150 (MQ=255)
                                                                                                                                     tGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAttt                  >  1:3280004/1‑150 (MQ=255)
                                                                                                                                      ggCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTg                 <  2:2229928/150‑1 (MQ=255)
                                                                                                                                                gtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAAt       <  2:3510406/150‑1 (MQ=255)
                                                                                                                                                     aaTGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTc  <  2:4501807/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTC  >  CP009273/274586‑274884

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