Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 410,300 C→A 23.0% A304S (GCG→TCG)  sbcC ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273410,3000CA23.0% 644.2 / 156.1 256A304S (GCG→TCG) sbcCexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (102/95);  new base A (28/31);  total (130/126)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 266 total aligned reads displayed.

CTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTTCTTGCTGCAATT  >  CP009273/410151‑410449
                                                                                                                                                     |                                                                                                                                                     
   tAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCt                                                                                                                                                    <  2:1067712/150‑1 (MQ=255)
    aaCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                                   >  1:1000697/1‑150 (MQ=255)
    aaCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtg                                                                                                                                                   >  1:639424/1‑150 (MQ=255)
     aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCtgt                                                                                                                                                  >  2:2096994/1‑150 (MQ=255)
     aCCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCtgt                                                                                                                                                  <  1:3394021/150‑1 (MQ=255)
      ccAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtgtg                                                                                                                                                 >  2:2314323/1‑150 (MQ=255)
      ccAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCAGCTGACTGCTTCGCCGCGAGGTGGCGAATGCTCGCGCGAAGCGCCATAGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCtgtg                                                                                                                                                 >  2:3055394/1‑150 (MQ=255)
       cAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCTCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTgtgt                                                                                                                                                >  1:2150266/1‑150 (MQ=255)
       cAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTgtgt                                                                                                                                                >  2:1632401/1‑150 (MQ=255)
        aGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGtt                                                                                                                                               >  2:4263347/1‑150 (MQ=255)
        aGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGtt                                                                                                                                               >  1:1026203/1‑150 (MQ=255)
         ggTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTc                                                                                                                                              <  2:1615638/150‑1 (MQ=255)
         ggTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTc                                                                                                                                              <  2:3736680/150‑1 (MQ=255)
           ttttCAGGATTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCTAATGCTCGCGCGACGCGCATTTGTGCTCTGTAAGCGAGTATTTAATTCTTCAATCTGCTGGCGAATTTGCGCCAGCGCCGCGCTGTGTTCTg                                                                                                                                            <  1:1516850/148‑1 (MQ=255)
            aTTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGc                                                                                                                                           >  2:2528073/1‑150 (MQ=255)
                tggCTTTGTTGCTGCTGGTTTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCTAATGCTCGCGCGACGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATTTGCGCCAGCGCCGCGCTGTGTTCTGCGATg                                                                                                                                       <  1:2835399/149‑1 (MQ=255)
                 ggCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGc                                                                                                                                      >  1:1647582/1‑150 (MQ=255)
                        ttgCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGACGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCa                                                                                                                               <  2:2824304/150‑1 (MQ=255)
                        ttgCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCa                                                                                                                               <  2:559457/150‑1 (MQ=255)
                            tgctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAgtgt                                                                                                                           >  1:296339/1‑150 (MQ=255)
                             gctgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGtgtg                                                                                                                          >  1:3304638/1‑150 (MQ=255)
                               tgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGa                                                                                                                        >  2:2583533/1‑150 (MQ=255)
                               tgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGa                                                                                                                        <  1:4422953/150‑1 (MQ=255)
                               tgctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGa                                                                                                                        >  1:669358/1‑150 (MQ=255)
                                gctgTAATTCTGCGGCGTGCTTGGCCGCGTGGGGGCGAATGCTCGCGCGAAGCGCATTTGGGCTCTGTAAGCGTGTTTTTACTTCTTCAATCTGGTGTCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGAc                                                                                                                       <  1:3503839/150‑1 (MQ=255)
                                 ctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACg                                                                                                                      <  2:100553/150‑1 (MQ=255)
                                 ctgTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACg                                                                                                                      >  2:2390193/1‑150 (MQ=255)
                                       ttCTGCTGACTGCGTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGAtt                                                                                                                <  2:2952733/150‑1 (MQ=255)
                                        tCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGAttt                                                                                                               <  2:3369285/150‑1 (MQ=255)
                                         ctgctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTc                                                                                                              >  1:2562869/1‑150 (MQ=255)
                                         ctgctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTCCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTc                                                                                                              <  1:3504718/150‑1 (MQ=255)
                                           gctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGt                                                                                                            <  2:2096480/150‑1 (MQ=255)
                                            ctgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTg                                                                                                           <  1:3397993/150‑1 (MQ=255)
                                             tgACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGc                                                                                                          >  1:2235377/1‑150 (MQ=255)
                                              gACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGcc                                                                                                         <  2:628320/150‑1 (MQ=255)
                                                  gCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCGGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGtt                                                                                                     <  2:3473951/150‑1 (MQ=255)
                                                     tCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTtgtg                                                                                                  >  2:2472525/1‑150 (MQ=255)
                                                        ccgcGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCa                                                                                               <  2:526423/150‑1 (MQ=255)
                                                         cgcgTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAg                                                                                              >  2:2608924/1‑150 (MQ=255)
                                                           cgTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGAc                                                                                            >  2:1814590/1‑150 (MQ=255)
                                                               gtggCGAATGCTCGCGCGAAGCGCAATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAg                                                                                        <  2:1708057/150‑1 (MQ=255)
                                                                tggCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAgc                                                                                       >  1:3814930/1‑150 (MQ=255)
                                                                 ggCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAgcg                                                                                      <  2:2456319/150‑1 (MQ=255)
                                                                   cGAATGCTCGCGCGACGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCGGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGcgcc                                                                                    <  1:1229750/150‑1 (MQ=255)
                                                                      aTGCTCGCGCGAAGCGCATTGGTGCTCTGTAAGCGTGTATTTAATTCTTCAATCTGCTGGCGAATTTGCGCCAGCGCCGCGCTGTGTTCTGCGTTGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCgccgcc                                                                                 <  1:292461/150‑1 (MQ=255)
                                                                       tGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCa                                                                                >  1:4430095/1‑150 (MQ=255)
                                                                             cgcgAAGCGCCATTGTGCTCGGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAg                                                                          <  1:2496312/150‑1 (MQ=255)
                                                                              gcgtcgcgcCATTGTGCTCTGTAAGCTTTTTTTTAATACTAAAATCGGCTGGCGAATAGGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAgg                                                                         <  1:1911318/145‑1 (MQ=255)
                                                                               cgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGt                                                                        >  2:4202970/1‑150 (MQ=255)
                                                                               cgAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCGGTGTTCTGCGATGCGTTCCCAGTGTGGACGACGATTTCGTGCCGGTTGTGCCAGACTCAGCGCCGCCAGTTGAGGt                                                                        >  2:1064370/1‑150 (MQ=255)
                                                                                   gcgcCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTgcg                                                                    >  2:2045068/1‑150 (MQ=255)
                                                                                     gcCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCt                                                                  >  1:4284484/1‑150 (MQ=255)
                                                                                      ccATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCtt                                                                 >  2:3295795/1‑150 (MQ=255)
                                                                                      ccATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCtt                                                                 >  1:3213745/1‑150 (MQ=255)
                                                                                      ccATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCtt                                                                 <  2:4429359/150‑1 (MQ=255)
                                                                                          tgtgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttc                                                             <  2:2150266/150‑1 (MQ=255)
                                                                                           gtgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttct                                                            >  1:3045776/1‑150 (MQ=255)
                                                                                           gtgCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttct                                                            <  1:2876359/150‑1 (MQ=255)
                                                                                              ctctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttct                                                         >  1:3562298/1‑150 (MQ=255)
                                                                                               tctGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTttcttctt                                                        <  2:1181449/150‑1 (MQ=255)
                                                                                                 tGTAAGCGATTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCg                                                      <  1:3990508/150‑1 (MQ=255)
                                                                                                  gTAGGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCgg                                                     >  1:711600/1‑150 (MQ=255)
                                                                                                    aaGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCt                                                   <  2:3863853/150‑1 (MQ=255)
                                                                                                    aaGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCt                                                   >  1:2405857/1‑150 (MQ=255)
                                                                                                     aGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTa                                                  <  2:3602264/150‑1 (MQ=255)
                                                                                                      gCGAGTATTTCCTTCTTCAATCGGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGGTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTaa                                                 <  1:4027616/150‑1 (MQ=255)
                                                                                                        gAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACg                                               <  1:2939239/150‑1 (MQ=255)
                                                                                                              ttACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGtt                                         <  2:4054435/150‑1 (MQ=255)
                                                                                                               tACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTg                                        >  1:1579067/1‑150 (MQ=255)
                                                                                                               tACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTg                                        >  2:1022349/1‑150 (MQ=255)
                                                                                                                 cttcttCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCa                                      >  2:1211275/1‑150 (MQ=255)
                                                                                                                   tcttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAg                                    <  2:1547968/150‑1 (MQ=255)
                                                                                                                   tcttcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAg                                    <  2:3829391/150‑1 (MQ=255)
                                                                                                                      tcAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGcc                                 <  1:1011146/150‑1 (MQ=255)
                                                                                                                       cAATCTGCTGGCGAATTTGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCct                                >  2:937343/1‑150 (MQ=255)
                                                                                                                       cAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCct                                >  1:635521/1‑150 (MQ=255)
                                                                                                                         aTCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgc                              <  1:3702888/150‑1 (MQ=255)
                                                                                                                         aTCTGCGGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgc                              <  1:4370296/150‑1 (MQ=255)
                                                                                                                          tCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGGTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGACAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgct                             <  2:2860796/150‑1 (MQ=255)
                                                                                                                           ctgctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgctg                            >  2:1446459/1‑150 (MQ=255)
                                                                                                                           cggctgGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCctgctg                            <  1:1161783/148‑1 (MQ=255)
                                                                                                                             gctgGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAc                          >  1:1897813/1‑150 (MQ=255)
                                                                                                                              ctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcg                         <  1:1487354/150‑1 (MQ=255)
                                                                                                                              ctgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcg                         <  1:3770260/150‑1 (MQ=255)
                                                                                                                              ctgGCGAATATGCGCCAGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcg                         >  1:3535113/1‑150 (MQ=255)
                                                                                                                               tgGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcgg                        >  1:83985/1‑150 (MQ=255)
                                                                                                                                ggCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcggc                       <  2:133729/150‑1 (MQ=255)
                                                                                                                                ggCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGAcggc                       <  1:2390193/150‑1 (MQ=255)
                                                                                                                                   gAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACggcggc                    >  2:857467/1‑150 (MQ=255)
                                                                                                                                   gAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACggcggc                    >  2:859320/1‑150 (MQ=255)
                                                                                                                                    aaTATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggct                   <  2:3791826/150‑1 (MQ=255)
                                                                                                                                     atatGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctg                  <  1:177640/150‑1 (MQ=255)
                                                                                                                                       atGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCggctggc                >  1:2068054/1‑150 (MQ=255)
                                                                                                                                          cgcCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGcttc             >  2:2719617/1‑150 (MQ=255)
                                                                                                                                              aGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTTCTtgc         >  1:4130253/1‑150 (MQ=255)
                                                                                                                                              aGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGcttctttc         >  1:3130105/1‑148 (MQ=255)
                                                                                                                                              aGCGCCGAGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTTCTtgc         >  2:2693979/1‑150 (MQ=255)
                                                                                                                                               gcgcCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTTCTtgct        <  1:2771142/150‑1 (MQ=255)
                                                                                                                                                 gccgcGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTTCTtgctgc      >  1:426685/1‑150 (MQ=255)
                                                                                                                                                   cgcgCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTTCTTGCTGCaa    >  1:473296/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTGTAACCAGGTATTCAGGCTTTGTTGCTGCTGCTGTAATTCTGCTGACTGCTTCGCCGCGTGGTGGCGAATGCTCGCGCGAAGCGCCATTGTGCTCTGTAAGCGAGTATTTACTTCTTCAATCTGCTGGCGAATATGCGCCAGCGCCGCGCTGTGTTCTGCGATGCGTTCCCAGTGTGGACGAAGATTTCGTGCCGGTTGTGCCAGACTAAGCGCCGCCAGTTGAGGTTGCGCTTTTTCTTCTTCGGCTAACGCCTGTTGCAAGGCCTGCTGACGGCGGCTGGCTTCTTGCTGCAATT  >  CP009273/410151‑410449

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: