Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 274,735 A→T 27.1% 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,7350AT27.1% 372.0 / 180.4 211pseudogene (138/348 nt)afuBpseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (87/67);  new base T (27/30);  total (114/97)
Fisher's exact test for biased strand distribution p-value = 2.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.08e-01

Only 100 of 216 total aligned reads displayed.

ATGAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACTC  >  CP009273/274586‑274884
                                                                                                                                                     |                                                                                                                                                     
  gAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAAt                                                                                                                                                     <  2:1541441/150‑1 (MQ=255)
  gAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAAt                                                                                                                                                     <  2:872219/150‑1 (MQ=255)
  gAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAAt                                                                                                                                                     >  2:1616338/1‑150 (MQ=255)
   aaaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATg                                                                                                                                                    >  2:2662945/1‑150 (MQ=255)
   aaaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATg                                                                                                                                                    >  1:4535873/1‑150 (MQ=255)
   aaaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATg                                                                                                                                                    >  1:1948505/1‑150 (MQ=255)
   aaaaTAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATg                                                                                                                                                    >  2:2778879/1‑150 (MQ=255)
     aatcatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGACCACAATGGCGCTGACGGTGGTACTGgc                                                                                                                                                  >  2:1611395/1‑150 (MQ=255)
     aataatCGCCAGCATCACCACGATCAGAATGGAAACGTAGGCAATCGCTACACCGTATTCGCCATCTTCCAAGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGgc                                                                                                                                                  >  2:552353/1‑150 (MQ=255)
       taatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGgcgc                                                                                                                                                <  1:1323473/150‑1 (MQ=255)
       taatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGgcgc                                                                                                                                                >  2:5010625/1‑150 (MQ=255)
       taatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGgcgc                                                                                                                                                >  2:3297444/1‑150 (MQ=255)
        aatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGgcgcg                                                                                                                                               <  2:4412741/150‑1 (MQ=255)
            gCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGcaca                                                                                                                                           <  1:3189691/150‑1 (MQ=255)
             ccAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACaa                                                                                                                                          >  2:4218382/1‑150 (MQ=255)
              cAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACaaa                                                                                                                                         >  2:2745186/1‑150 (MQ=255)
               aGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAg                                                                                                                                        >  2:3884740/1‑150 (MQ=255)
                 cATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCt                                                                                                                                      <  1:2144000/150‑1 (MQ=255)
                 cATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCt                                                                                                                                      >  2:4900402/1‑150 (MQ=255)
                  aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTg                                                                                                                                     <  2:4244164/150‑1 (MQ=255)
                  aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTg                                                                                                                                     <  2:4012048/150‑1 (MQ=255)
                  aTCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTg                                                                                                                                     >  1:4493570/1‑150 (MQ=255)
                   tCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGt                                                                                                                                    >  2:4868200/1‑150 (MQ=255)
                    caccacGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTa                                                                                                                                   >  1:560451/1‑150 (MQ=255)
                           aTCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                            <  1:5170680/150‑1 (MQ=255)
                            tCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAgcg                                                                                                                           <  1:3377749/150‑1 (MQ=255)
                             cAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAgcgc                                                                                                                          <  2:3321147/150‑1 (MQ=255)
                              aGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGcc                                                                                                                         >  1:3696599/1‑150 (MQ=255)
                                 aTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCgag                                                                                                                      >  2:2165040/1‑150 (MQ=255)
                                      aCCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGAt                                                                                                                 <  2:113998/150‑1 (MQ=255)
                                        cGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCg                                                                                                               <  2:1320293/150‑1 (MQ=255)
                                        cGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCg                                                                                                               <  2:4739416/150‑1 (MQ=255)
                                          tAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGcc                                                                                                             <  1:662958/150‑1 (MQ=255)
                                            ggCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCgg                                                                                                           >  1:2731084/1‑150 (MQ=255)
                                              cAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGgc                                                                                                         >  1:4032952/1‑150 (MQ=255)
                                               aaTCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCTCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGAGCGCCGGgcg                                                                                                        >  1:1189683/1‑150 (MQ=255)
                                               aaTCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGgcg                                                                                                        >  1:2411692/1‑150 (MQ=255)
                                                aTCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGgcgc                                                                                                       >  1:3133922/1‑150 (MQ=255)
                                                 tCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgca                                                                                                      <  2:3405518/150‑1 (MQ=255)
                                                    cTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagca                                                                                                   <  1:4156820/150‑1 (MQ=255)
                                                    cTACACCATATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcagca                                                                                                   >  1:2180399/1‑150 (MQ=255)
                                                       cacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCg                                                                                                <  2:381374/150‑1 (MQ=255)
                                                        acCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCgg                                                                                               >  2:3518924/1‑150 (MQ=255)
                                                          cGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgca                                                                                             >  1:5015291/1‑150 (MQ=255)
                                                          cGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgca                                                                                             <  1:840413/150‑1 (MQ=255)
                                                           gTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcag                                                                                            >  2:4544010/1‑150 (MQ=255)
                                                           gTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcag                                                                                            >  1:2784692/1‑150 (MQ=255)
                                                             aTTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcagca                                                                                          <  2:943926/150‑1 (MQ=255)
                                                             aTTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcagca                                                                                          >  1:4194294/1‑150 (MQ=255)
                                                              ttCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGcagcag                                                                                         >  2:2540643/1‑150 (MQ=255)
                                                               tCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAgg                                                                                        <  1:2504930/150‑1 (MQ=255)
                                                                  ncATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGCGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtg                                                                                     >  1:3132586/2‑150 (MQ=255)
                                                                  ccATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtg                                                                                     <  2:1686994/150‑1 (MQ=255)
                                                                  ccATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtg                                                                                     <  1:4461789/150‑1 (MQ=255)
                                                                   cATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtgt                                                                                    <  1:4894662/150‑1 (MQ=255)
                                                                    aTCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtgtg                                                                                   >  2:3472323/1‑150 (MQ=255)
                                                                      cTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCg                                                                                 >  2:1034464/1‑150 (MQ=255)
                                                                       ttCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGt                                                                                >  2:79555/1‑150 (MQ=255)
                                                                       ttCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGt                                                                                >  2:79411/1‑150 (MQ=255)
                                                                          cACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAAt                                                                             <  2:2251265/150‑1 (MQ=255)
                                                                               ggTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTAc                                                                        <  2:2098011/150‑1 (MQ=255)
                                                                                 ttGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGc                                                                      >  1:10749/1‑150 (MQ=255)
                                                                                   gAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAg                                                                    >  1:3041058/1‑150 (MQ=255)
                                                                                    aGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTCATGGTACGCAgg                                                                   >  2:4775690/1‑150 (MQ=255)
                                                                                        tGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAAc                                                               <  2:2383800/150‑1 (MQ=255)
                                                                                        tGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAAc                                                               <  1:1616338/150‑1 (MQ=255)
                                                                                          tAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCg                                                             <  2:1401830/150‑1 (MQ=255)
                                                                                          tAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCg                                                             <  2:3144388/150‑1 (MQ=255)
                                                                                            aGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgc                                                           >  2:2339184/1‑150 (MQ=255)
                                                                                                gTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGACCCGgcgcgc                                                       >  2:4109125/1‑150 (MQ=255)
                                                                                                 tCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTAGGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCa                                                      <  2:1800165/150‑1 (MQ=255)
                                                                                                 tCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCa                                                      <  1:1213168/150‑1 (MQ=255)
                                                                                                  cGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAg                                                     >  2:1900525/1‑150 (MQ=255)
                                                                                                    cTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGc                                                   >  1:3172997/1‑150 (MQ=255)
                                                                                                     tACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCt                                                  <  2:3219922/150‑1 (MQ=255)
                                                                                                      aCGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtg                                                 >  2:1494522/1‑150 (MQ=255)
                                                                                                      aCGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtg                                                 >  1:548320/1‑150 (MQ=255)
                                                                                                          gggTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagtg                                             <  1:594449/150‑1 (MQ=255)
                                                                                                           ggTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagtga                                            <  1:2935711/150‑1 (MQ=255)
                                                                                                           ggTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagtga                                            <  1:2433582/150‑1 (MQ=255)
                                                                                                              aTCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGc                                         <  1:4310503/150‑1 (MQ=255)
                                                                                                               tCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTCCGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCt                                        >  2:2478052/1‑150 (MQ=255)
                                                                                                                 cGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTc                                      <  2:636765/150‑1 (MQ=255)
                                                                                                                     gTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCg                                  >  2:4465881/1‑150 (MQ=255)
                                                                                                                       gACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAg                                >  2:1173075/1‑150 (MQ=255)
                                                                                                                               acacAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTAt                        >  1:3573045/1‑150 (MQ=255)
                                                                                                                                cacaATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATc                       >  2:1196744/1‑150 (MQ=255)
                                                                                                                                cacaATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATc                       <  1:2552944/150‑1 (MQ=255)
                                                                                                                                 acaATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCg                      <  1:4003348/150‑1 (MQ=255)
                                                                                                                                    aTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAtt                   <  2:2011688/150‑1 (MQ=255)
                                                                                                                                     tGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAttt                  <  2:5039463/150‑1 (MQ=255)
                                                                                                                                      ggCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTg                 >  2:3944999/1‑150 (MQ=255)
                                                                                                                                      ggCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTg                 <  2:3511757/150‑1 (MQ=255)
                                                                                                                                            gACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCgg           >  1:2205459/1‑150 (MQ=255)
                                                                                                                                              cGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTa         <  1:4195930/150‑1 (MQ=255)
                                                                                                                                                gtggtAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAAt       <  2:4386577/150‑1 (MQ=255)
                                                                                                                                                 tggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATg      >  1:1768029/1‑150 (MQ=255)
                                                                                                                                                   gtAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGAc    >  2:2185829/1‑150 (MQ=255)
                                                                                                                                                    tAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGACt   <  1:3276558/150‑1 (MQ=255)
                                                                                                                                                     aaTGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCGGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCAACGAGCGATTTATCGATTTGGCCCGGTAATGACTc  >  1:1760241/1‑150 (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: