Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 260,917 T→G 5.5% M1M (ATG→CTG) † yafX ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273260,9170TG 5.5% 656.9 / 2.8 209M1M (ATG→CTG) †yafXCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (103/88);  new base G (3/15);  total (106/103)
Fisher's exact test for biased strand distribution p-value = 2.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

ATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTT  >  CP009273/260769‑261065
                                                                                                                                                      |                                                                                                                                                    
   tGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATggg                                                                                                                                                   >  2:2559506/1‑149 (MQ=255)
    gCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGt                                                                                                                                                  <  2:1407427/149‑1 (MQ=255)
    gCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGt                                                                                                                                                  <  1:2510086/149‑1 (MQ=255)
        cTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGt                                                                                                                                              <  2:2056975/149‑1 (MQ=255)
           cAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGtctc                                                                                                                                           <  1:4494876/149‑1 (MQ=255)
            aGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTcc                                                                                                                                          <  2:3702647/149‑1 (MQ=255)
            aGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTcc                                                                                                                                          <  2:3584987/149‑1 (MQ=255)
            aGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTcc                                                                                                                                          >  2:2931723/1‑149 (MQ=255)
             ggAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCg                                                                                                                                         >  2:3432202/1‑149 (MQ=255)
                  aGCCACTGCGGGATAGCGCCACAGTACTGTGGCCAGACATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAGATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCGGTGTTGTCATGGGTCTGTCTCCGTCAAt                                                                                                                                    >  1:4044076/1‑149 (MQ=255)
                   gCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATa                                                                                                                                   >  1:1661432/1‑149 (MQ=255)
                    ccACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaa                                                                                                                                  <  1:1881823/149‑1 (MQ=255)
                     cACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaa                                                                                                                                 <  2:892916/149‑1 (MQ=255)
                     cACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaa                                                                                                                                 <  2:2413213/149‑1 (MQ=255)
                       cTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATaaaaa                                                                                                                               >  1:101886/1‑149 (MQ=255)
                       cTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCAGGGGTATGTCTCCGTCAATaaaaa                                                                                                                               >  1:3878941/1‑149 (MQ=255)
                        tGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGGCCGGTCAAAGTCGGGTGGGTTAGCGGGTGCGAGGTCGTGCTGCGTCGGTGTTGTCATGGGTCTGTCTCCGTCAATCAAAAc                                                                                                                              >  2:2601528/1‑149 (MQ=255)
                           gggATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGcc                                                                                                                           >  1:2816956/1‑149 (MQ=255)
                           gggATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGcc                                                                                                                           <  2:2331155/149‑1 (MQ=255)
                            ggATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCa                                                                                                                          <  1:605333/149‑1 (MQ=255)
                            ggATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCa                                                                                                                          >  2:3683478/1‑149 (MQ=255)
                             gATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAg                                                                                                                         >  1:235698/1‑149 (MQ=255)
                             gATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAg                                                                                                                         >  1:1986833/1‑149 (MQ=255)
                                     cAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGGGTGACGGTCTTTTAAAATTCGGTTTGGTTAGGGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTAATGGTTTTGTCTCCGTCAATAAAAACGCCAGCGGCGAtg                                                                                                                 <  2:2812396/149‑1 (MQ=21)
                                     cAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtg                                                                                                                 >  2:372525/1‑149 (MQ=255)
                                      aaaGTTCGGTGGCCATAAAGCTATACGCGGTTCTTGGGGGACCGGGGTTAGGTTCTTTTCAATTGGGGTTGGGTTAGGGGGTGCTTGGTCGTGCTGCTTTTTTTTTGTCATGGTTTTGTCTCCTTAAATAAAAACGCCAGCGGCGAtgg                                                                                                                <  2:3122017/149‑1 (MQ=255)
                                       aaGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtggc                                                                                                               >  2:858207/1‑149 (MQ=255)
                                       aaGTACTGGGGCCAGAAACCGATACGCTGTTCAGCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGGGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtggc                                                                                                               <  2:4450940/149‑1 (MQ=255)
                                         gTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtggctg                                                                                                             >  2:1974140/1‑149 (MQ=255)
                                          tACTGTGGCCATAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGAtggctgg                                                                                                            <  2:1068735/149‑1 (MQ=255)
                                             tggggCCAGAAATCGATACGCTGTTCATGGGGGACCGGCGTGACGGTCAGTTCAAATTGGGGTTGGTTAGGGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGt                                                                                                         <  1:3811998/146‑1 (MQ=21)
                                                   cAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCAGGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATgggg                                                                                                   >  2:2311510/1‑149 (MQ=255)
                                                     gAAATCTATACGGTGTTCATGGGGGCCCGGCGTGCCGGTCAGTTCAATTTGGGGTTGGTTTGCGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCTTGGGTTGGTCTCCGTCAATAAACACGCCAGCGGCGATGGCTGGCGTATGGGGat                                                                                                 <  1:2201484/149‑1 (MQ=21)
                                                            aTACGCTGTTCTTGGGGGACGGGCGTGACGTTCTTTTAAAATTGGGGTTGTTTAGGGGGTGCGTGGTGGTGCTGCTTCTGTGTTGTCTTGGTTTTGTCTCCTAAAAAAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgt                                                                                          <  1:658527/149‑1 (MQ=16)
                                                             tACGGTGTTCATGGGGGACCGGGGTTAGGTTATTTTAAATTTGGGTTTGTTTTGGGGGTGGGAGGTGTTGCTGGGTCTGTGTTGTCTTGGTTTTTTATCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgtg                                                                                         <  1:361943/149‑1 (MQ=12)
                                                             tACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgtg                                                                                         >  2:1915134/1‑149 (MQ=255)
                                                             tACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGACAATAAAAAAGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgtg                                                                                         >  1:2718080/1‑149 (MQ=255)
                                                              acggggTTCATCGGGGACCGGCGTGACGGTCAGTTAAATTGCGGTTGGGTTTGGGGGTGCGCGGTCGGGCTGCTTCGTTGTTGTAATGGGTAGGTATCCTTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgtgt                                                                                        <  2:3868928/144‑1 (MQ=21)
                                                              aCGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgtgt                                                                                        >  2:3026055/1‑149 (MQ=255)
                                                              aCGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgtgt                                                                                        >  2:3947447/1‑149 (MQ=255)
                                                               cGCTGTTTAGCGGGCACGGGGTTGAGGGTCTTTTCAATTTGGGGTTTTTTCGGGGGTGGGTGGTGTTGGTGGTTGTTTGTTTTCTGGGTTTTGTGCCCTTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAgtgtg                                                                                       <  1:4430608/149‑1 (MQ=255)
                                                                 cTGTTCCTCGGGGACGGGCGTGACGGTCAGTTCAATTTGGGGTTGGTTCGGGGGTGCGTGGTGTGGCTGGGTCTGTGTTGTCTTGGGTATGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGtt                                                                                     <  2:101886/149‑1 (MQ=21)
                                                                 cTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGtt                                                                                     <  2:2555215/149‑1 (MQ=255)
                                                                  tGTTCAGGGGGGACGGGGGTGAGGTTCTGTTAAATTTGGGGTTGTTTAGGGGGTGCGTGGTGGTGCTGGGTGTGTGTTGTATTGGTTATTTCTCGTTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTc                                                                                    <  1:2919930/149‑1 (MQ=16)
                                                                   gttattCGGGGACCGGCGTGACTGTCAGTTAAAATTGGGTTTGTTTAGCGGGTGCTAGGTCTTGCTGCGTCTGTGTTGTAATTGGTATGTATCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCg                                                                                   <  2:1692695/144‑1 (MQ=21)
                                                                     tCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCggg                                                                                 >  2:2694290/1‑149 (MQ=255)
                                                                      cttCGGGGACCGGCGTGACGGTCAGTTCAAATGCGGGTTGGTTCGGGGGTGCGTGGTGGTGCTGGGTGGGTGTTGTCATGGGTATGTCTCCGTCAATAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCgggg                                                                                <  2:490916/147‑1 (MQ=25)
                                                                      caggggggACCGGCGTGACGTTCAGTTCAATTTGGGTTTGGTTCGGGGGTGGGTGGTGGTGCTGGGTCTGTGTTGTATTGGTTTTGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCgggg                                                                                <  1:2126571/145‑1 (MQ=21)
                                                                        tCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAg                                                                              >  1:4310236/1‑149 (MQ=255)
                                                                          ggggACCGGCGTGACGGTCTTTTCAAAGTGGGGTTGGTTTGCGGGTGCGAGGTCGTGCTGCGTCTGTGGGGTCTTGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGt                                                                            <  2:3918057/149‑1 (MQ=25)
                                                                             gACCGGCGTGAGGGTCAGTTCAATTTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGaa                                                                         <  2:3294866/149‑1 (MQ=255)
                                                                               ccGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATg                                                                       >  2:1173311/1‑149 (MQ=255)
                                                                                cGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGc                                                                      >  1:1701757/1‑149 (MQ=255)
                                                                                   cGTGAGGGTCTGTTCAAATGGGGGTTGGTTCGGGGGTGCGTGGTGGTGCTGCGTCTGTGTTGTCTTGGTTTTGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCggg                                                                   <  1:4060873/149‑1 (MQ=21)
                                                                                     tGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTa                                                                 >  2:2187612/1‑149 (MQ=255)
                                                                                           tCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGACTGCGGGTAAATGCt                                                           >  2:1555708/1‑149 (MQ=255)
                                                                                            cAGTTCAATTTCGGGTTGGTTCGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCtt                                                          <  2:3535044/149‑1 (MQ=255)
                                                                                            cAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCtt                                                          >  2:4205560/1‑149 (MQ=255)
                                                                                                   atggggggTTGGTTGGGGGGTGGGTGGTGGTGCTGCGTGTGTTTTGTCTTGGTTTGGCCCCCTAAAAAAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATc                                                   <  1:3886521/144‑1 (MQ=21)
                                                                                                    aTTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCa                                                  >  2:2675802/1‑149 (MQ=255)
                                                                                                     ttCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAg                                                 >  2:2281388/1‑149 (MQ=255)
                                                                                                     ttCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAg                                                 >  2:37796/1‑149 (MQ=255)
                                                                                                      gggggTTGGTTGGGGGGTGGGTGGTGGGGGGGGGTGGGTGTGGTCTGGGTTTTGTCTCCGAAAAAAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGc                                                <  1:3468978/147‑1 (MQ=19)
                                                                                                                  tAGCGGGTGCGTGGTCGTGCTGCGTCTGTGTTGTCTTGGGTATGTCTCCGAAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCgt                                      <  1:2083002/149‑1 (MQ=255)
                                                                                                                   tgCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCgtg                                     <  1:3304231/148‑1 (MQ=255)
                                                                                                                        gTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGt                                >  2:1632477/1‑149 (MQ=255)
                                                                                                                         tGCGTGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAACGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTa                               <  2:3364162/149‑1 (MQ=255)
                                                                                                                         tGCGAGTTGGTGCTGCTTGGGTGTTTTCTGGGGTATGTCTCCGACAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTa                               <  1:3271393/149‑1 (MQ=37)
                                                                                                                         tGCGAGGTCGTGCTGCGTCTGTGTTGTCTTGGGTCTGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTa                               <  1:3683478/149‑1 (MQ=255)
                                                                                                                         tGCGAGGTCGTGCTGCGTCTGTGTTGTCCTGGGAATGTCTCCGACAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTa                               <  2:786982/149‑1 (MQ=255)
                                                                                                                         tGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTa                               >  2:793210/1‑149 (MQ=255)
                                                                                                                         tGCGAGGTCGTGCGGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTa                               <  2:4008558/149‑1 (MQ=255)
                                                                                                                         tGCGAGGTCGGGCTGCGTCTGTGTTGTCAGGGGTAGGTCACCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTa                               <  1:971667/149‑1 (MQ=255)
                                                                                                                             cggTCGTGCTGCGTCTGTGTTGGCTGGGGTATGCCCCCGTCAACAAACACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGa                           <  1:112826/148‑1 (MQ=38)
                                                                                                                             aGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGa                           >  2:3980881/1‑149 (MQ=255)
                                                                                                                              ggtggtGCTGCGTGGGTGTTTTCTGGGGTTTGTCTCCTTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGAt                          <  2:3446671/149‑1 (MQ=37)
                                                                                                                              ggTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGAt                          >  2:3046939/1‑149 (MQ=255)
                                                                                                                               gTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATg                         >  1:4389140/1‑149 (MQ=255)
                                                                                                                                tCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGc                        >  1:3132642/1‑149 (MQ=255)
                                                                                                                                tCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTAAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGc                        <  1:1417586/149‑1 (MQ=255)
                                                                                                                                  gggctggGTCTGTGTGGACAGGGGTATGTCGCCTAAAATAACAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCt                      <  1:78329/147‑1 (MQ=37)
                                                                                                                                   tgctgcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTg                     >  2:2359450/1‑149 (MQ=255)
                                                                                                                                   tgctgcGTCTGTGTTGACATGGGTATGTCTCCGTCAAAAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTg                     <  1:740268/149‑1 (MQ=255)
                                                                                                                                    gctgcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGa                    >  1:148534/1‑149 (MQ=255)
                                                                                                                                      tgcGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAgg                  <  2:2081451/149‑1 (MQ=255)
                                                                                                                                      tgcGGCTGTGTTGTCATGGCTATTTCTCCGTCAATAAAAACGACAGCGGCGATGGCTGGCGTATGGGGATAGAAAGTGTGTTCTGGGAGGTGACTGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTAAGGATGCCTGAgg                  >  2:153576/1‑149 (MQ=38)
                                                                                                                                         gTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGc               >  2:1638945/1‑149 (MQ=255)
                                                                                                                                         gTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGc               <  2:4344924/149‑1 (MQ=255)
                                                                                                                                         gTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGc               >  2:4543343/1‑149 (MQ=255)
                                                                                                                                          tCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCg              >  1:3064609/1‑149 (MQ=255)
                                                                                                                                          tCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCg              >  2:2488031/1‑149 (MQ=255)
                                                                                                                                              tgttgtCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGata          >  1:81962/1‑149 (MQ=255)
                                                                                                                                              tgttgtCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGata          >  2:3627806/1‑149 (MQ=255)
                                                                                                                                              tgttgtCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGata          >  1:3716921/1‑149 (MQ=255)
                                                                                                                                                  gtCATGGGTATGGCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGc      <  2:3687710/149‑1 (MQ=255)
                                                                                                                                                   tCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCg     >  2:3380428/1‑149 (MQ=255)
                                                                                                                                                     aTGGGTATGTCTCCGTCAATAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGtt  >  1:3553565/1‑149 (MQ=255)
                                                                                                                                                      tGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGtt  <  2:1595382/149‑1 (MQ=255)
                                                                                                                                                      tGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGtt  <  2:3763601/149‑1 (MQ=255)
                                                                                                                                                      tGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGtt  <  2:1860089/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
ATATGCGGCTCCAGGAGTAGCCACTGCGGGATAGCGCCAAAGTACTGTGGCCAGAAATCGATACGCTGTTCATCGGGGACCGGCGTGACGGTCAGTTCAAATTCGGGTTGGTTAGCGGGTGCGAGGTCGTGCTGCGTCTGTGTTGTCATGGGTATGTCTCCGTCAATAAAAACGCCAGCGGCGATGGCTGGCGTATGGGGATATAAAGTGTGTTCGGGGAGGTGAATGCGGGTAAATGCTTCGCGATCAGCGGGTGGCCGTGTCTGTACGGATGCCTGAGGTGCGGATATAGCGGTT  >  CP009273/260769‑261065

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