Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,192,525 C→T 5.0% S155N (AGC→AAC)  ymfD ← e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,192,5250CT 5.0% 1147.7 / 31.4 320S155N (AGC→AAC) ymfDe14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (155/148);  new base T (10/6);  total (166/154)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 327 total aligned reads displayed.

CGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCG  >  CP009273/1192377‑1192674
                                                                                                                                                    |                                                                                                                                                     
cGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCt                                                                                                                                                      >  2:467973/1‑150 (MQ=255)
 gACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCtt                                                                                                                                                     >  2:2075449/1‑150 (MQ=255)
   cTTTAATATTTCAAGGCCATGTGAAGAACATATTTATTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCtttt                                                                                                                                                   <  2:1939463/150‑1 (MQ=255)
    tttAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTta                                                                                                                                                  <  2:3542960/150‑1 (MQ=255)
     atAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTtat                                                                                                                                                 <  2:389666/149‑1 (MQ=255)
           tttCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGAt                                                                                                                                           >  1:4894364/1‑150 (MQ=255)
             tCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGAttt                                                                                                                                         <  1:459651/150‑1 (MQ=255)
              cAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTa                                                                                                                                        <  2:2385890/150‑1 (MQ=255)
               aaGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTAc                                                                                                                                       >  1:2407351/1‑150 (MQ=255)
                aGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACa                                                                                                                                      >  2:5364110/1‑150 (MQ=255)
                 ggCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTTATTTACaa                                                                                                                                     >  1:4435943/1‑150 (MQ=255)
                 ggCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACaa                                                                                                                                     >  1:2484536/1‑150 (MQ=255)
                 ggCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTGTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTTAAGTATGAGCTTTTATATTGATTTACaa                                                                                                                                     >  1:1225891/1‑150 (MQ=255)
                  gCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                    >  2:2975797/1‑150 (MQ=255)
                  gCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                    >  1:4239308/1‑150 (MQ=255)
                  gCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                    <  2:4337670/150‑1 (MQ=255)
                  gCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAAt                                                                                                                                    <  1:4178682/150‑1 (MQ=255)
                     aTGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAag                                                                                                                                 <  1:2446530/150‑1 (MQ=255)
                     aTGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAag                                                                                                                                 <  1:5161887/150‑1 (MQ=255)
                        tgAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTGTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGt                                                                                                                              >  1:2816250/1‑150 (MQ=255)
                          aagaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGtct                                                                                                                            >  2:2178110/1‑150 (MQ=255)
                           agaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGtctc                                                                                                                           <  1:5445361/150‑1 (MQ=255)
                           agaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGtctc                                                                                                                           <  1:3231848/150‑1 (MQ=255)
                           agaaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGtctc                                                                                                                           >  1:3892581/1‑150 (MQ=255)
                             aaCATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCAc                                                                                                                         >  2:4491491/1‑150 (MQ=255)
                                     tCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAt                                                                                                                 <  1:2002679/150‑1 (MQ=255)
                                      cTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtt                                                                                                                >  2:1853647/1‑150 (MQ=255)
                                        tGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtttt                                                                                                              <  1:1020824/150‑1 (MQ=255)
                                        tGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtttt                                                                                                              <  2:3109841/150‑1 (MQ=255)
                                        tGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGAtttt                                                                                                              >  1:2589256/1‑150 (MQ=255)
                                           aCTGCGAGTTCATCGAGTAATCCATATTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGt                                                                                                           <  2:4948552/150‑1 (MQ=255)
                                            cTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTa                                                                                                          <  2:5205773/150‑1 (MQ=255)
                                            cTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTa                                                                                                          <  2:1332465/150‑1 (MQ=255)
                                            cTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTa                                                                                                          <  2:4064808/150‑1 (MQ=255)
                                                gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTatgaaatg                                                                                                      >  1:1334378/1‑143 (MQ=255)
                                                gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                      <  2:5392446/150‑1 (MQ=255)
                                                gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                      >  1:3835146/1‑150 (MQ=255)
                                                gAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtct                                                                                                      >  1:5423650/1‑150 (MQ=255)
                                                 aGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAAtctc                                                                                                     >  1:667515/1‑150 (MQ=255)
                                                   ttCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCtt                                                                                                   >  1:2144084/1‑150 (MQ=255)
                                                     cATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGa                                                                                                 <  2:2479964/150‑1 (MQ=255)
                                                       tCGTGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTATATTTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTTATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATg                                                                                               <  2:3805092/150‑1 (MQ=255)
                                                         gAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCt                                                                                             <  1:2646250/150‑1 (MQ=255)
                                                         gAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCt                                                                                             >  2:4784351/1‑150 (MQ=255)
                                                          aGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTa                                                                                            >  1:582074/1‑150 (MQ=255)
                                                            tAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAg                                                                                          >  1:2918475/1‑150 (MQ=255)
                                                              aTCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGaa                                                                                        >  1:411829/1‑150 (MQ=255)
                                                               tCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAAc                                                                                       <  2:3191419/150‑1 (MQ=255)
                                                                       aaGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGt                                                                               <  1:1444765/150‑1 (MQ=255)
                                                                         gAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCATGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCg                                                                             >  2:2725050/1‑150 (MQ=255)
                                                                          aaaCACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGa                                                                            >  2:2831235/1‑150 (MQ=255)
                                                                               cTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTTTCGATTGTa                                                                       >  1:5074119/1‑150 (MQ=255)
                                                                               cTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTa                                                                       <  2:3511038/150‑1 (MQ=255)
                                                                               cTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTa                                                                       >  1:467626/1‑150 (MQ=255)
                                                                                  tttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCTATTGTATcc                                                                    >  2:5171464/1‑150 (MQ=255)
                                                                                  tttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATcc                                                                    >  2:1498921/1‑150 (MQ=255)
                                                                                   ttGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATccc                                                                   >  1:818783/1‑150 (MQ=255)
                                                                                      aaTTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGAc                                                                >  2:1056663/1‑150 (MQ=255)
                                                                                      aaTTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGAc                                                                >  2:2234136/1‑150 (MQ=255)
                                                                                       aTTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACa                                                               <  2:3935773/150‑1 (MQ=255)
                                                                                       aTTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACa                                                               >  2:4272276/1‑150 (MQ=255)
                                                                                          tttGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAgg                                                            <  2:4605976/150‑1 (MQ=255)
                                                                                            tGTAAATGTATCCGTAAAGATTTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGaa                                                          <  1:2782806/150‑1 (MQ=255)
                                                                                            tGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGaa                                                          <  1:2086469/150‑1 (MQ=255)
                                                                                             gTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTAACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAAc                                                         <  1:2672253/150‑1 (MQ=255)
                                                                                               aaaTGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACgg                                                       >  1:1792452/1‑150 (MQ=255)
                                                                                                 aTGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCa                                                     >  2:3602899/1‑150 (MQ=255)
                                                                                                   gTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCaga                                                   <  1:1707248/150‑1 (MQ=255)
                                                                                                    tATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCagag                                                  >  1:2271936/1‑150 (MQ=255)
                                                                                                     aTCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATAGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGAATGTATCCCGACAAGGAACGGCagaga                                                 >  2:648942/1‑150 (MQ=255)
                                                                                                       ccGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAgagaga                                               >  1:2631212/1‑150 (MQ=255)
                                                                                                         gTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAAc                                             >  2:4845067/1‑150 (MQ=255)
                                                                                                         gTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGATCTATTTTGTGGATAGTATCCCGACCAGGAACGGCAgagagagc                                             >  2:3242824/1‑148 (MQ=255)
                                                                                                         gTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAAc                                             >  2:491430/1‑150 (MQ=255)
                                                                                                         gTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAAc                                             <  1:3107990/150‑1 (MQ=255)
                                                                                                          tAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACa                                            <  2:2568021/150‑1 (MQ=255)
                                                                                                            aaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGACCAtt                                          >  2:2279791/1‑150 (MQ=255)
                                                                                                            aaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAtt                                          >  1:675676/1‑150 (MQ=255)
                                                                                                            aaGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACAtt                                          >  2:5097414/1‑150 (MQ=255)
                                                                                                             aGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTa                                         <  2:1265387/150‑1 (MQ=255)
                                                                                                                   tttCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCaa                                   >  1:915042/1‑150 (MQ=255)
                                                                                                                   tttCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCaa                                   <  2:4597057/150‑1 (MQ=255)
                                                                                                                     tcgtcCTGGTTCGGATTGTTGGACGTCTGCGCGTGTATAGAGATTTCCACTCAGAGTCGCACCTCCTGATTGTAGTCCTCTCTTGATGCTAAGAACTATGTGGTCGATTGGATCCCGCCAAGGAACGGCAGAGAGAACCTTAGAGCACAg                                 <  2:3293533/147‑1 (MQ=255)
                                                                                                                      cttcCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAgg                                <  2:2142665/150‑1 (MQ=255)
                                                                                                                         ccTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATg                             >  2:1815133/1‑150 (MQ=255)
                                                                                                                         ccTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATg                             >  1:4309156/1‑150 (MQ=255)
                                                                                                                         ccTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATg                             <  1:330679/150‑1 (MQ=255)
                                                                                                                          cTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGa                            <  2:4239237/150‑1 (MQ=255)
                                                                                                                               tCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATTAAATCg                       >  1:1117983/1‑150 (MQ=255)
                                                                                                                                    ttttttGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAAcc                  >  2:1530465/1‑150 (MQ=255)
                                                                                                                                           aaGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatt           >  2:3055182/1‑150 (MQ=255)
                                                                                                                                           aaGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatt           <  1:3636023/150‑1 (MQ=255)
                                                                                                                                           aaGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatt           <  1:1670253/150‑1 (MQ=255)
                                                                                                                                            aGTAAGAGCATTTATATTTATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAatta          <  2:1168685/150‑1 (MQ=255)
                                                                                                                                              tATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAagtatt        >  1:3107426/1‑150 (MQ=255)
                                                                                                                                                 gAGTTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTg     <  2:4432571/150‑1 (MQ=255)
                                                                                                                                                   gCTTTTATTTTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTc   <  2:2885716/150‑1 (MQ=255)
                                                                                                                                                   gCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTc   >  1:3282668/1‑150 (MQ=255)
                                                                                                                                                   gCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTc   <  2:890834/150‑1 (MQ=255)
                                                                                                                                                    cTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCg  <  1:2364198/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGACTTTAATATTTCAAGGCCATGTGAAGAACATATTTCTTGCACTGCGAGTTCATCGAGTAATCCATAGTAAGAAACACTTTTTGAATTTTTGTAAATGTATCCGTAAAGATGTTTTCTTCCTGTTTCGTATTTTTTGAAGTATGAGCTTTTATATTGATTTACAATAAGAGTCTCACCTCCTGATTTTAGTAATCTCTTGATGCTAAGAACTATTTTGTCGATTGTATCCCGACAAGGAACGGCAGAGAGAACATTAGAGCAAAGGATGAAATCGTAACCCCCAATTATTTTGTCG  >  CP009273/1192377‑1192674

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