Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,753,064 G→A 6.4% D241N (GAC→AAC)  yfjI → CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,753,0640GA 6.4% 1008.6 / 37.5 282D241N (GAC→AAC) yfjICP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (136/128);  new base A (8/10);  total (144/138)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 283 total aligned reads displayed.

CACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAA  >  CP009273/2752916‑2753211
                                                                                                                                                    |                                                                                                                                                   
cacacTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCg                                                                                                                                                     >  2:910304/1‑149 (MQ=255)
  cacTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGAc                                                                                                                                                   >  1:2145945/1‑149 (MQ=255)
  cacTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGAc                                                                                                                                                   <  1:2899760/149‑1 (MQ=255)
    cTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATAGACAg                                                                                                                                                 >  1:208822/1‑149 (MQ=255)
      tAACGAGCTGCCTTTTATCAATAAGATGTGGTATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGa                                                                                                                                               <  1:4922364/149‑1 (MQ=255)
       aaCGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCAACAGGaa                                                                                                                                              >  1:4304988/1‑149 (MQ=255)
        aCGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCCGCCGAATGTTTTTAAGGGTTATATCGACAGGaaa                                                                                                                                             >  1:545702/1‑149 (MQ=255)
        aCGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGaaa                                                                                                                                             <  1:3273815/149‑1 (MQ=255)
        aCGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCAACAGGaaa                                                                                                                                             >  2:112073/1‑149 (MQ=255)
          gAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATAGACAGGAAAgg                                                                                                                                           >  1:2397078/1‑149 (MQ=255)
           aGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGga                                                                                                                                          <  1:3784645/149‑1 (MQ=255)
            gCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGgag                                                                                                                                         <  2:4684212/149‑1 (MQ=255)
              tGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCGAATGTTTTTAAGGGTTATATAGACAGGAAAGGAGat                                                                                                                                       >  1:4364983/1‑149 (MQ=255)
               gCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGata                                                                                                                                      <  1:4990314/149‑1 (MQ=255)
                 cttttATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATg                                                                                                                                    <  1:3389923/147‑1 (MQ=255)
                      aTCAATATGTTGTGGTATGTTTAAATATTTACGGTGGAAAGTAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCaaa                                                                                                                               <  1:3046610/149‑1 (MQ=255)
                      aTCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCaaa                                                                                                                               >  2:1187057/1‑149 (MQ=255)
                      aTCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCaaa                                                                                                                               <  2:2655979/149‑1 (MQ=255)
                       tCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAg                                                                                                                              >  1:1480107/1‑149 (MQ=255)
                        cAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAgg                                                                                                                             <  2:535125/149‑1 (MQ=255)
                         aaTAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAggg                                                                                                                            <  2:1242727/149‑1 (MQ=255)
                         aaTAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAggg                                                                                                                            <  1:217041/149‑1 (MQ=255)
                          aTAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAgggg                                                                                                                           <  2:2585263/149‑1 (MQ=255)
                            aaGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTATTTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGggat                                                                                                                         <  1:5129208/149‑1 (MQ=255)
                                tgtgGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGggattgga                                                                                                                     <  2:1895264/149‑1 (MQ=255)
                                tgtgGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCAACAGGAAAGGAGATATGGCAAAGGggattgga                                                                                                                     <  2:1298617/149‑1 (MQ=255)
                                  tgGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGggattggatt                                                                                                                   <  1:2877500/149‑1 (MQ=255)
                                  tgGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGggattggatt                                                                                                                   <  1:4641240/149‑1 (MQ=255)
                                        ggTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGc                                                                                                             >  2:4853669/1‑149 (MQ=255)
                                         gTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCa                                                                                                            <  1:625312/149‑1 (MQ=255)
                                             aaTATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCAACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGt                                                                                                        <  2:1786807/149‑1 (MQ=255)
                                              atatTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTg                                                                                                       <  1:3877267/149‑1 (MQ=255)
                                                atTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGcc                                                                                                     >  2:1620619/1‑149 (MQ=255)
                                                   tACGGTGGAAAGGAAAAACGAGCCCGGGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCa                                                                                                  <  1:824708/149‑1 (MQ=255)
                                                      ggTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCAtgt                                                                                               <  1:1129299/149‑1 (MQ=255)
                                                         ggaaaggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGcc                                                                                            >  1:1613822/1‑149 (MQ=255)
                                                         ggaaaggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGcc                                                                                            <  2:2439030/149‑1 (MQ=255)
                                                             aggaaaAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGcc                                                                                        <  2:365849/149‑1 (MQ=255)
                                                                 aaaaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCt                                                                                    <  1:861615/149‑1 (MQ=255)
                                                                  aaaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTTCACGGTGCCTCATGTGCCAGCCTGCtt                                                                                   >  1:3948633/1‑149 (MQ=255)
                                                                  aaaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCAACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCtt                                                                                   >  1:126678/1‑149 (MQ=255)
                                                                   aaCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTTCACGGTGCCTCATGTGCCAGCCTGCTTc                                                                                  >  2:2094115/1‑149 (MQ=255)
                                                                    aCGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCa                                                                                 >  2:309628/1‑149 (MQ=255)
                                                                      gAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAac                                                                               <  2:4252593/149‑1 (MQ=255)
                                                                        gCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAacac                                                                             >  1:3593479/1‑149 (MQ=255)
                                                                           cGAGAAATTTATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTACTGTTTTTAAGGTTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAg                                                                          <  1:4864752/149‑1 (MQ=255)
                                                                             agaAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGt                                                                        <  1:48395/149‑1 (MQ=255)
                                                                               aaaTTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTaa                                                                      >  2:4498625/1‑149 (MQ=255)
                                                                               aaaTTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTaa                                                                      <  1:2282163/149‑1 (MQ=255)
                                                                                aattaattAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAAc                                                                     <  2:975849/149‑1 (MQ=255)
                                                                                  ttaattaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAg                                                                   >  1:564549/1‑149 (MQ=255)
                                                                                  ttaattaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAg                                                                   >  2:2905691/1‑149 (MQ=255)
                                                                                   taattaGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGa                                                                  <  2:2377720/149‑1 (MQ=255)
                                                                                     attaGAGATGCAAGAATAACACTGTCGCTGATGGTACAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGaaa                                                                >  2:1366674/1‑149 (MQ=255)
                                                                                         gagaTGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAAtt                                                            >  1:2984240/1‑149 (MQ=255)
                                                                                            aTGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCAACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTcc                                                         <  1:941406/149‑1 (MQ=255)
                                                                                             tGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCa                                                        >  2:3868109/1‑149 (MQ=255)
                                                                                             tGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCa                                                        >  2:1996830/1‑149 (MQ=255)
                                                                                              gCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCaa                                                       <  2:622991/149‑1 (MQ=255)
                                                                                               cAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAAc                                                      <  1:2613613/149‑1 (MQ=255)
                                                                                               cAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAAc                                                      <  1:1651173/149‑1 (MQ=255)
                                                                                                aaGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAAcc                                                     >  1:1672119/1‑149 (MQ=255)
                                                                                                aaGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAAcc                                                     >  2:3929594/1‑149 (MQ=255)
                                                                                                    aTAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAt                                                 >  1:2747042/1‑149 (MQ=255)
                                                                                                     tAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAtt                                                >  2:4075772/1‑149 (MQ=255)
                                                                                                     tAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAtt                                                >  2:826823/1‑149 (MQ=255)
                                                                                                     tAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAtt                                                >  2:3425692/1‑149 (MQ=255)
                                                                                                     tAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAttt                                               >  2:223056/1‑149 (MQ=255)
                                                                                                         acTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAAtttttt                                            <  2:2745151/149‑1 (MQ=255)
                                                                                                          cTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTc                                           <  2:3685812/149‑1 (MQ=255)
                                                                                                          cTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATAGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAATGTAACAGAAAAATTTCCAACACAAATTTTTc                                           >  2:4475566/1‑149 (MQ=255)
                                                                                                           tGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCa                                          >  2:4408588/1‑149 (MQ=255)
                                                                                                              cGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCACTTTTTTCAAAt                                       >  2:4358082/1‑149 (MQ=255)
                                                                                                              cGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAAt                                       <  2:3503242/149‑1 (MQ=255)
                                                                                                              cGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAAt                                       <  2:2183275/149‑1 (MQ=255)
                                                                                                                cTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGa                                     >  2:4431572/1‑149 (MQ=255)
                                                                                                                 tGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGaa                                    >  1:978565/1‑149 (MQ=255)
                                                                                                                  gATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAAc                                   >  1:352422/1‑149 (MQ=255)
                                                                                                                   aTGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACa                                  <  1:2005724/149‑1 (MQ=255)
                                                                                                                    tGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAt                                 <  1:1055011/149‑1 (MQ=255)
                                                                                                                     ggTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAtt                                <  2:1477411/149‑1 (MQ=255)
                                                                                                                     ggTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAtt                                <  2:5015878/149‑1 (MQ=255)
                                                                                                                         cAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGcc                            >  1:4117890/1‑149 (MQ=255)
                                                                                                                             cTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTa                        >  1:256379/1‑149 (MQ=255)
                                                                                                                              tAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTAt                       >  2:3830998/1‑149 (MQ=255)
                                                                                                                              tAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTAt                       >  1:769216/1‑149 (MQ=255)
                                                                                                                               aaTGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTAtt                      >  1:4186510/1‑149 (MQ=255)
                                                                                                                               aaTGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTAtt                      >  1:2047487/1‑149 (MQ=255)
                                                                                                                               aaTGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTAtt                      >  2:1671637/1‑149 (MQ=255)
                                                                                                                                aTGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTAttt                     <  1:2259774/149‑1 (MQ=255)
                                                                                                                                 tGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTc                    <  1:4455775/149‑1 (MQ=255)
                                                                                                                                   tttttAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACAATTGCCGGTATTTCAc                  >  2:436840/1‑149 (MQ=255)
                                                                                                                                         aGGGTTATATCAACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGt            >  2:3469089/1‑149 (MQ=255)
                                                                                                                                             ttATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTa        >  2:1600384/1‑149 (MQ=255)
                                                                                                                                                tatCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATgg     >  1:2951706/1‑149 (MQ=255)
                                                                                                                                                 atCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGa    >  1:3168674/1‑149 (MQ=255)
                                                                                                                                                 atCGACAGGAAAGGAGATATGGCAAAGGGGATTGGAATTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGACCATTTGCCGGTATTTCACCAACGTCTTATGGa    >  1:4535706/1‑149 (MQ=255)
                                                                                                                                                  tCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaa   <  2:2175506/149‑1 (MQ=255)
                                                                                                                                                   cGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaaa  <  2:2471143/149‑1 (MQ=255)
                                                                                                                                                   cGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGaaa  <  2:315151/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACACTTAACGAGCTGCCGTTTATCAATAAGATGTGGGATGGTTCAATATTTACGGTGGAAAGGAAAAACGAGCCCGAGAAATTAATTAGAGATGCAAGAATAACACTGTCGCTGATGGTCCAGCCTAATGTTTTTAAGGGTTATATCGACAGGAAAGGAGATATGGCAAAGGGGATTGGATTTTTTGCACGGTGCCTCATGTGCCAGCCTGCTTCAACACAAGGTAACAGAAAAATTTCCAACCCAATTTTTTCAAATGAACATTTGCCGGTATTTCACCAACGTCTTATGGAAA  >  CP009273/2752916‑2753211

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: