Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,553,563 T→G 17.1% intergenic (+138/‑53) intZ → / → yffL CPZ‑55 prophage; putative phage integrase/CPZ‑55 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,553,5630TG17.1% 638.4 / 148.7 276intergenic (+138/‑53)intZ/yffLCPZ‑55 prophage; putative phage integrase/CPZ‑55 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (127/101);  new base G (33/14);  total (161/115)
Fisher's exact test for biased strand distribution p-value = 7.50e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

Only 100 of 316 total aligned reads displayed.

GAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAG  >  CP009273/2553414‑2553712
                                                                                                                                                     |                                                                                                                                                     
gAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattatt                                                                                                                                                       <  1:2638583/150‑1 (MQ=255)
  aCCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattcc                                                                                                                                                     <  2:1303559/150‑1 (MQ=255)
    cGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattccat                                                                                                                                                   >  2:5330339/1‑150 (MQ=255)
     ggAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattccatt                                                                                                                                                  >  2:1027304/1‑150 (MQ=255)
     ggAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattccatt                                                                                                                                                  >  1:2421510/1‑150 (MQ=255)
        atgaatAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAattattccattatt                                                                                                                                               >  2:2387165/1‑150 (MQ=255)
           aataatTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAAttattccattattcca                                                                                                                                            >  2:680085/1‑150 (MQ=255)
           aataatTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTGTTTCGCTAAttattccattattcca                                                                                                                                            >  1:4629116/1‑150 (MQ=255)
            ataatTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCaa                                                                                                                                           >  2:3576103/1‑150 (MQ=255)
            ataatTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCaa                                                                                                                                           >  1:3038386/1‑150 (MQ=255)
               atTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTa                                                                                                                                        >  2:1968510/1‑150 (MQ=255)
               atTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTAGTCCATTATTCCAATTa                                                                                                                                        <  1:5669881/150‑1 (MQ=255)
                 tAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAg                                                                                                                                      <  1:2982572/150‑1 (MQ=255)
                     cGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGa                                                                                                                                  >  2:5239998/1‑150 (MQ=255)
                       gCTTTTTTGCAAATCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAAt                                                                                                                                <  1:2356900/150‑1 (MQ=255)
                       gCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAAt                                                                                                                                <  1:2331684/150‑1 (MQ=255)
                         ttttttGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATaa                                                                                                                              >  1:3662973/1‑150 (MQ=255)
                           ttttGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAat                                                                                                                            >  2:4829474/1‑150 (MQ=255)
                           ttgtGCAACGCGCAATTTTGCATTTTGTATAAGTCTTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAat                                                                                                                            <  1:3814490/147‑1 (MQ=255)
                            tttGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAatt                                                                                                                           <  1:1739629/150‑1 (MQ=255)
                             ttGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAatta                                                                                                                          >  2:1457686/1‑150 (MQ=255)
                                 aaaGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCa                                                                                                                      >  1:2303649/1‑150 (MQ=255)
                                     cgcAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAatat                                                                                                                  >  2:2115026/1‑150 (MQ=255)
                                       cAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCa                                                                                                                >  1:4478915/1‑150 (MQ=255)
                                        aaTTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCaa                                                                                                               >  1:4892702/1‑150 (MQ=255)
                                         aTTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCAATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAg                                                                                                              >  2:3408847/1‑150 (MQ=255)
                                          ttttGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTAAATATCAAGa                                                                                                             >  2:3912881/1‑150 (MQ=255)
                                           tttGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGaa                                                                                                            <  2:5027938/150‑1 (MQ=255)
                                            ttGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAAc                                                                                                           <  2:1977484/150‑1 (MQ=255)
                                            ttGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAAc                                                                                                           <  2:4562966/150‑1 (MQ=255)
                                             tGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACt                                                                                                          >  2:4079144/1‑150 (MQ=255)
                                               cATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTa                                                                                                        >  1:1890765/1‑150 (MQ=255)
                                                     gTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATcaa                                                                                                  >  2:3168315/1‑148 (MQ=255)
                                                      tataAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCAt                                                                                                 >  1:1030379/1‑150 (MQ=255)
                                                      tataAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCAt                                                                                                 >  1:734154/1‑150 (MQ=255)
                                                       ataAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATg                                                                                                >  1:4037933/1‑150 (MQ=255)
                                                        taAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGt                                                                                               >  1:1336187/1‑150 (MQ=255)
                                                          aGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGttt                                                                                             >  2:1888193/1‑150 (MQ=255)
                                                             aTTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCa                                                                                          >  1:4962532/1‑150 (MQ=255)
                                                                tgttttgtGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACg                                                                                       >  2:2803638/1‑150 (MQ=255)
                                                                 gttttgtGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGTTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAATAACTTAAATCCATGTTTGCAACGa                                                                                      >  1:1936182/1‑150 (MQ=255)
                                                                  ttttgtGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGaa                                                                                     >  2:1756645/1‑150 (MQ=255)
                                                                  ttttgtGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGaa                                                                                     >  1:5408432/1‑150 (MQ=255)
                                                                   tttgtGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGaaa                                                                                    <  2:3449676/150‑1 (MQ=255)
                                                                     tgtgtgGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGa                                                                                  <  1:1468363/150‑1 (MQ=255)
                                                                     tgtgtgGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGa                                                                                  >  1:3295611/1‑150 (MQ=255)
                                                                      gtgtgGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGTTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGAt                                                                                 >  1:2711810/1‑150 (MQ=255)
                                                                      gtgtgGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGAt                                                                                 >  1:1287772/1‑150 (MQ=255)
                                                                        gtgGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAAGATCAAGAACTTAAATCCATGTTTGCAACGAAAGATcc                                                                               <  2:2040713/150‑1 (MQ=255)
                                                                         tgGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCa                                                                              <  2:2336346/150‑1 (MQ=255)
                                                                          ggTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCag                                                                             <  1:255387/150‑1 (MQ=255)
                                                                             ttAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGt                                                                          <  1:1577890/150‑1 (MQ=255)
                                                                             ttAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGt                                                                          >  1:3296556/1‑150 (MQ=255)
                                                                              tAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGtt                                                                         >  1:2646910/1‑150 (MQ=255)
                                                                               aaCTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTc                                                                        <  2:2191428/150‑1 (MQ=255)
                                                                               aaCTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAACTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTc                                                                        <  2:1694321/150‑1 (MQ=255)
                                                                               aaCTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTc                                                                        >  2:828359/1‑150 (MQ=255)
                                                                                 cTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGa                                                                      >  1:903709/1‑150 (MQ=255)
                                                                                  tATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCTGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGaa                                                                     <  2:5632149/150‑1 (MQ=255)
                                                                                   aTCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGaaa                                                                    <  2:3940430/150‑1 (MQ=255)
                                                                                   aTCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGaaa                                                                    >  2:4347630/1‑150 (MQ=255)
                                                                                     cACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAAt                                                                  <  1:2391523/150‑1 (MQ=255)
                                                                                      aCCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATa                                                                 <  2:4407336/150‑1 (MQ=255)
                                                                                       ccTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAg                                                                <  1:3576103/150‑1 (MQ=255)
                                                                                        cTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAgg                                                               <  1:344698/150‑1 (MQ=255)
                                                                                         tGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGa                                                              >  2:1062766/1‑150 (MQ=255)
                                                                                          gATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGTCAATCGGAt                                                             >  1:4168927/1‑150 (MQ=255)
                                                                                           aTAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATc                                                            >  2:3710254/1‑150 (MQ=255)
                                                                                                 aaCCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAAcaca                                                      >  1:5068989/1‑150 (MQ=255)
                                                                                                   ccTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACaaa                                                    >  1:514957/1‑150 (MQ=255)
                                                                                                      gTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGTTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATaa                                                 <  1:277312/150‑1 (MQ=255)
                                                                                                      gTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATaa                                                 >  2:1811864/1‑150 (MQ=255)
                                                                                                        aTTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATaaaa                                               >  2:378289/1‑150 (MQ=255)
                                                                                                         ttCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAAGAGGATCAACACAAATAAAAg                                              >  1:2952040/1‑150 (MQ=255)
                                                                                                         ttCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATACATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAg                                              >  1:81292/1‑150 (MQ=255)
                                                                                                          tctaTTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGc                                             <  2:4296555/147‑1 (MQ=255)
                                                                                                          tCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGc                                             <  1:5330823/150‑1 (MQ=255)
                                                                                                            aaTTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGccc                                           >  2:2614252/1‑150 (MQ=255)
                                                                                                            aaTTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGccc                                           >  1:1761320/1‑150 (MQ=255)
                                                                                                            aaTTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGccc                                           >  1:3453664/1‑150 (MQ=255)
                                                                                                             aTTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAATAGcccc                                          >  2:4935233/1‑150 (MQ=255)
                                                                                                              tttGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGccccc                                         <  2:5365388/150‑1 (MQ=255)
                                                                                                                    ggATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtg                                   <  1:386856/150‑1 (MQ=255)
                                                                                                                     gATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgg                                  <  2:2499562/150‑2 (MQ=255)
                                                                                                                      aTTCGTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgct                                 <  2:4686742/150‑1 (MQ=255)
                                                                                                                        tCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAAtgctgc                               <  1:2090419/150‑1 (MQ=255)
                                                                                                                          gTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAACATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                             >  2:2019536/1‑150 (MQ=255)
                                                                                                                          gTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                             >  1:4809800/1‑150 (MQ=255)
                                                                                                                          gTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                             >  2:2266618/1‑150 (MQ=255)
                                                                                                                          gTCCGGCTTTTCTTTTTCGCTAATTATGCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCaa                             >  2:3677639/1‑150 (MQ=255)
                                                                                                                            ccGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACa                           >  1:3949374/1‑150 (MQ=255)
                                                                                                                                    tcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                   >  1:1126634/1‑150 (MQ=255)
                                                                                                                                    tcttttTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCgggg                   <  1:5531587/150‑1 (MQ=255)
                                                                                                                                      tttttCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGc                 >  1:523784/1‑150 (MQ=255)
                                                                                                                                       ttttCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCg                <  1:3793122/150‑1 (MQ=255)
                                                                                                                                          tCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa             <  1:1442585/150‑1 (MQ=255)
                                                                                                                                          tCGCTAATTATTCCATTATTCCAACTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTa             >  1:5190663/1‑150 (MQ=255)
                                                                                                                                                attattccattattccAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGaa       >  1:649651/1‑150 (MQ=255)
                                                                                                                                                  tattccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAAc     <  1:2179512/150‑1 (MQ=255)
                                                                                                                                                   attccattattccaATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACa    <  2:4895083/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACCGGAATGAATAATTAAGCGGCTTTTTTGCAAAGCGCAATTTTGCATTTTGTATAAGTATTTGTTTTGTGTGGTTTAACTATCACCTGATAAGAAACCTGTATTCAATTTGCCGGATTCGTCCGGCTTTTCTTTTTCGCTAATTATTCCATTATTCCAATTAAGTGGAATAAATTATTCAATATCAAGAACTTAAATCCATGTTTGCAACGAAAGATCCAGAGTTCGAAAATAGGATCAACACAAATAAAAGCCCCCGGAATGCTGCAACATGCCGGGGGCGGTATGAGAAACAAG  >  CP009273/2553414‑2553712

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: