Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 274,735 A→T 89.6% 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,7350AT89.6% 1049.9 / 85.3 324pseudogene (138/348 nt)afuBpseudogene, CP4‑6 prophage; ABC transporter permease family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base A (23/11);  new base T (159/130);  total (182/142)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

ATGAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGA  >  CP009273/274586‑274881
                                                                                                                                                     |                                                                                                                                                  
aTGAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTa                                                                                                                                                    <  2:3797790/150‑1 (MQ=255)
aTGAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTa                                                                                                                                                    >  2:774683/1‑150 (MQ=255)
 tGAAAATAATCGCCAGCATCACCACGATCAGAATGGAAACGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGAGGTTGAGTATGTAAGCCGTAGCTACGAGGGTATCCGGCGTTACGAGGAACACAATGGCGCTGAAGGTGGTaa                                                                                                                                                   >  1:2376769/1‑150 (MQ=255)
     aataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgc                                                                                                                                               <  1:1737387/150‑1 (MQ=255)
     aataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgc                                                                                                                                               <  2:2585101/150‑1 (MQ=255)
     aataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGgc                                                                                                                                               <  2:206884/150‑1 (MQ=255)
      ataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcg                                                                                                                                              >  1:513658/1‑150 (MQ=255)
      ataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcg                                                                                                                                              >  1:4273137/1‑150 (MQ=255)
      ataatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcg                                                                                                                                              >  1:84158/1‑150 (MQ=255)
       taatCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGgcgc                                                                                                                                             <  1:4944494/150‑1 (MQ=255)
            gCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGcaca                                                                                                                                        >  2:1398741/1‑150 (MQ=255)
              cAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACaaa                                                                                                                                      >  2:5113000/1‑150 (MQ=255)
               aGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAg                                                                                                                                     >  2:284122/1‑150 (MQ=255)
                        acGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATc                                                                                                                            >  2:3657424/1‑150 (MQ=255)
                          gATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAg                                                                                                                          >  2:2646763/1‑150 (MQ=255)
                          gATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAg                                                                                                                          >  1:388/1‑150 (MQ=255)
                           aTCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                         <  1:989158/150‑1 (MQ=255)
                           aTCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                         >  2:2377019/1‑150 (MQ=255)
                           aTCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                         <  2:2999285/150‑1 (MQ=255)
                           aTCAGAATGGAACCGTAGGCAATAGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgc                                                                                                                         >  1:689834/1‑150 (MQ=255)
                            tCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcg                                                                                                                        >  2:1969655/1‑150 (MQ=255)
                             cAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACACTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAgcgc                                                                                                                       >  1:1455232/1‑150 (MQ=255)
                                aaTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCga                                                                                                                    <  2:3182603/150‑1 (MQ=255)
                                aaTGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCga                                                                                                                    >  1:4444568/1‑150 (MQ=255)
                                   ggAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGagag                                                                                                                 >  1:2893462/1‑150 (MQ=255)
                                   ggAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGagag                                                                                                                 >  1:71813/1‑150 (MQ=255)
                                    gAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAgg                                                                                                                >  1:3700951/1‑150 (MQ=255)
                                     aaCCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGa                                                                                                               >  1:3542496/1‑150 (MQ=255)
                                     aaCCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGa                                                                                                               >  2:3045483/1‑150 (MQ=255)
                                       ccGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATc                                                                                                             <  2:2281574/150‑1 (MQ=255)
                                        cGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCg                                                                                                            >  1:3065115/1‑150 (MQ=255)
                                          tAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGcc                                                                                                          <  1:1222347/150‑1 (MQ=255)
                                          tAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTAGCCGGCGTGACGAGGAACACACTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGcc                                                                                                          <  2:565784/150‑1 (MQ=255)
                                            ggCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCgg                                                                                                        <  1:112455/150‑1 (MQ=255)
                                            ggCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCgg                                                                                                        <  1:4210193/150‑1 (MQ=255)
                                             gCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCggg                                                                                                       >  2:4324932/1‑150 (MQ=255)
                                                 tCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgca                                                                                                   <  1:557083/150‑1 (MQ=255)
                                                  cGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCgcag                                                                                                  <  2:4278576/150‑1 (MQ=255)
                                                     tACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGcagcag                                                                                               >  2:1519298/1‑150 (MQ=255)
                                                     tACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGcagcag                                                                                               >  2:4517799/1‑150 (MQ=255)
                                                      acacCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCagcagc                                                                                              >  2:2389588/1‑150 (MQ=255)
                                                         ccGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgc                                                                                           <  1:4018218/150‑1 (MQ=255)
                                                          cGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgca                                                                                          >  1:1286134/1‑150 (MQ=255)
                                                          cGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgca                                                                                          <  1:3158978/150‑1 (MQ=255)
                                                             aTTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGgcagca                                                                                       >  2:70063/1‑150 (MQ=255)
                                                               tCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAgg                                                                                     >  1:1874998/1‑150 (MQ=255)
                                                               tCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAgg                                                                                     >  1:2006371/1‑150 (MQ=255)
                                                               tCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAgg                                                                                     <  2:4427079/150‑1 (MQ=255)
                                                                  ccATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtg                                                                                  >  2:3593617/1‑150 (MQ=255)
                                                                  ccATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtg                                                                                  >  1:2562338/1‑150 (MQ=255)
                                                                  ccATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGAAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtg                                                                                  <  1:2681467/150‑1 (MQ=255)
                                                                   cATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtgt                                                                                 <  1:3345180/150‑1 (MQ=255)
                                                                    aTCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGAtgtg                                                                                >  1:1555016/1‑150 (MQ=255)
                                                                     tCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGc                                                                               >  2:2128508/1‑150 (MQ=255)
                                                                        tCCACGCGGGTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTa                                                                            >  2:4627147/1‑150 (MQ=255)
                                                                                gTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACg                                                                    >  1:1076128/1‑150 (MQ=255)
                                                                                 ttGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGc                                                                   >  2:239931/1‑150 (MQ=255)
                                                                                 ttGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGc                                                                   >  2:4188138/1‑150 (MQ=255)
                                                                                  tGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCa                                                                  >  1:4452677/1‑150 (MQ=255)
                                                                                   gAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGGTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCCGCGGCAGCCGGATGTGCGTAATGGTACGCAg                                                                 >  2:796824/1‑150 (MQ=255)
                                                                                     ggATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAggg                                                               <  1:2563875/150‑1 (MQ=255)
                                                                                       aTGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGaa                                                             >  1:3312824/1‑150 (MQ=255)
                                                                                        tGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAAc                                                            >  2:3250148/1‑150 (MQ=255)
                                                                                        tGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAAc                                                            >  2:3893889/1‑150 (MQ=255)
                                                                                          tAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCg                                                          <  2:165667/150‑1 (MQ=255)
                                                                                           aaGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCgg                                                         <  2:1475815/150‑1 (MQ=255)
                                                                                           aaGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCgg                                                         <  1:3348818/150‑1 (MQ=255)
                                                                                            aGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgc                                                        <  2:1153088/150‑1 (MQ=255)
                                                                                               cGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGgcgcg                                                     >  2:2947770/1‑150 (MQ=255)
                                                                                                  cGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAg                                                  >  2:1828090/1‑150 (MQ=255)
                                                                                                  cGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAg                                                  >  1:743093/1‑150 (MQ=255)
                                                                                                   gCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAgg                                                 <  2:4131122/150‑1 (MQ=255)
                                                                                                          gggTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagtg                                          >  2:2976504/1‑150 (MQ=255)
                                                                                                           ggTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagtga                                         <  1:4761886/150‑1 (MQ=255)
                                                                                                            gTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCtgagtgag                                        >  1:2301813/1‑150 (MQ=255)
                                                                                                               tCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCt                                     <  1:1029999/150‑1 (MQ=255)
                                                                                                               tCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCt                                     >  1:2759290/1‑150 (MQ=255)
                                                                                                                ccGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCtt                                    >  1:4116291/1‑150 (MQ=255)
                                                                                                                    cGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATc                                >  2:472484/1‑150 (MQ=255)
                                                                                                                     gTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCg                               >  1:1582689/1‑150 (MQ=255)
                                                                                                                     gTGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCg                               <  1:1655350/150‑1 (MQ=255)
                                                                                                                      tGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGa                              >  2:3366350/1‑150 (MQ=255)
                                                                                                                      tGACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGa                              >  2:4838595/1‑150 (MQ=255)
                                                                                                                       gACGAGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAg                             <  2:5044323/150‑1 (MQ=255)
                                                                                                                       gACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAg                             >  1:2233029/1‑150 (MQ=255)
                                                                                                                           aGGAACACAATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGAt                         >  1:4201065/1‑150 (MQ=255)
                                                                                                                             gAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGAttt                       >  2:4961632/1‑150 (MQ=255)
                                                                                                                                cacaCTGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGAAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATc                    <  1:1922574/150‑1 (MQ=255)
                                                                                                                                cacaATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATc                    >  1:3660837/1‑150 (MQ=255)
                                                                                                                                 acaATGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCg                   <  1:1304865/150‑1 (MQ=255)
                                                                                                                                     tGGCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGAttt               <  2:3450878/150‑1 (MQ=255)
                                                                                                                                      ggCGCTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTg              >  2:1881842/1‑150 (MQ=255)
                                                                                                                                         gcTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGcc           <  1:1539536/150‑1 (MQ=255)
                                                                                                                                          cTGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGccc          >  1:2513256/1‑150 (MQ=255)
                                                                                                                                           tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg         <  1:4697158/150‑1 (MQ=255)
                                                                                                                                           tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg         >  1:4495109/1‑150 (MQ=255)
                                                                                                                                           tGACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCg         <  1:4798272/150‑1 (MQ=255)
                                                                                                                                            gACGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCgg        >  1:874979/1‑150 (MQ=255)
                                                                                                                                              cGGTGGTTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGAAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTa      >  1:5022680/1‑150 (MQ=255)
                                                                                                                                               ggtggtTATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTaa     >  2:4057401/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGAAAATAATCGCCAGCATCACCACGATCAGAATGGAACCGTAGGCAATCGCTACACCGTATTCGCCATCTTCCACGCGGTTGAGGATGTAAGCCGTCGCTACGCGGGTATCCGGCGTGACGAGGAACACAATGGCGCTGACGGTGGTAATGGCGCGCACAAAGCTGTAAATCAGCGCCGAGAGGATCGCCGGGCGCAGCAGCGGCAGCAGGATGTGCGTAATGGTACGCAGGGAACCGGCGCGCAGGCTGAGTGAGGCTTCATCGAGCGATTTATCGATTTGGCCCGGTAATGA  >  CP009273/274586‑274881

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: