Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,508,852 G→T 8.6% P141T (CCC→ACC)  dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,508,8520GT 8.6% 873.2 / 39.6 233P141T (CCC→ACC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base G (106/107);  new base T (13/7);  total (119/114)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 241 total aligned reads displayed.

GGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAG  >  CP009273/3508703‑3509000
                                                                                                                                                     |                                                                                                                                                    
ggTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAggg                                                                                                                                                      <  2:2217310/150‑1 (MQ=255)
    gCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGtttt                                                                                                                                                  >  2:2869976/1‑150 (MQ=255)
     cAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGttttt                                                                                                                                                 <  2:2113996/150‑1 (MQ=255)
        aCAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGt                                                                                                                                              <  1:3178282/150‑1 (MQ=255)
         cAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAACTAGGGTTTTTTGTa                                                                                                                                             >  2:4009903/1‑150 (MQ=255)
         cAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTa                                                                                                                                             <  2:1061376/150‑1 (MQ=255)
              gTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAACTAGGTTTTTTTGTAGCGGc                                                                                                                                        >  2:3849338/1‑150 (MQ=255)
              gTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGc                                                                                                                                        <  1:4081449/150‑1 (MQ=255)
               tGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGcc                                                                                                                                       <  1:3249809/150‑1 (MQ=255)
                gCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAACTAGGGATTTTTGTAGCGGCCg                                                                                                                                      >  1:3754858/1‑150 (MQ=255)
                gCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCg                                                                                                                                      <  2:556381/150‑1 (MQ=255)
                   tAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAAc                                                                                                                                   <  1:355077/150‑1 (MQ=255)
                   tAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAAc                                                                                                                                   <  1:2228807/150‑1 (MQ=255)
                    aaGGCGGATCGCAATAGACGACGGATGCATCATCAGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCTCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACg                                                                                                                                  <  2:4124207/150‑1 (MQ=255)
                     aGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACgg                                                                                                                                 >  1:422142/1‑150 (MQ=255)
                      ggCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGc                                                                                                                                >  2:3450630/1‑150 (MQ=255)
                      ggCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGc                                                                                                                                <  1:2839948/150‑1 (MQ=255)
                       gCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCa                                                                                                                               <  2:1761385/150‑1 (MQ=255)
                        cGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCAc                                                                                                                              <  2:712062/150‑1 (MQ=255)
                        cGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCAc                                                                                                                              >  1:1386779/1‑150 (MQ=255)
                             cGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTCTCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTaa                                                                                                                         >  1:1127547/1‑150 (MQ=255)
                              gCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTaaa                                                                                                                        >  2:2755013/1‑150 (MQ=255)
                                aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                      >  1:26319/1‑150 (MQ=255)
                                aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                      <  1:1572028/150‑1 (MQ=255)
                                aaTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACt                                                                                                                      <  2:1457361/150‑1 (MQ=255)
                                 aTAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGTCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTc                                                                                                                     >  2:2042560/1‑150 (MQ=255)
                                  tAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCa                                                                                                                    >  1:1498544/1‑150 (MQ=255)
                                    gacgacGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAcc                                                                                                                  >  1:3013045/1‑150 (MQ=255)
                                    gacgacGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAcc                                                                                                                  <  2:2337844/150‑1 (MQ=255)
                                     acgacgGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACcg                                                                                                                 <  2:1951614/150‑1 (MQ=255)
                                          ggATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAg                                                                                                            >  2:800895/1‑150 (MQ=255)
                                              gCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTg                                                                                                        <  2:436742/150‑1 (MQ=255)
                                                atcatcTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTa                                                                                                      <  2:3157787/150‑1 (MQ=255)
                                                  catcTGCGCGCGCCATGCTATCGGCGTAAGACTCACATTAGAAACAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAAc                                                                                                    <  2:667917/150‑1 (MQ=255)
                                                            cgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGcc                                                                                          <  1:2278212/150‑1 (MQ=255)
                                                            cgcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGcc                                                                                          <  1:3165535/150‑1 (MQ=255)
                                                             gcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCg                                                                                         >  1:1953536/1‑150 (MQ=255)
                                                             gcCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCg                                                                                         <  1:11859/150‑1 (MQ=255)
                                                               cATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGtt                                                                                       >  1:2374374/1‑150 (MQ=255)
                                                                   cTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTaa                                                                                   >  2:1235842/1‑150 (MQ=255)
                                                                    tATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAAc                                                                                  <  2:1371318/150‑1 (MQ=255)
                                                                     aTCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGACATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAAcc                                                                                 >  2:1269507/1‑150 (MQ=255)
                                                                     aTCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAAcc                                                                                 <  1:571943/150‑1 (MQ=255)
                                                                     aTCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAAcc                                                                                 <  1:246689/150‑1 (MQ=255)
                                                                       cGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGt                                                                               <  1:2986427/150‑1 (MQ=255)
                                                                          cGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGc                                                                            >  2:1416561/1‑150 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCgg                                                                          >  2:2004594/1‑150 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCgg                                                                          >  1:3353151/1‑150 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCgg                                                                          >  2:3271148/1‑150 (MQ=255)
                                                                            tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGAAACAACTCTGCTTCCGGGAAATAGTGTTTTTTGTAGCTGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCgg                                                                          >  2:3128737/1‑150 (MQ=255)
                                                                              aGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGtt                                                                        >  1:275806/1‑150 (MQ=255)
                                                                               gACTCACAATAGAAACAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTc                                                                       <  2:825399/150‑1 (MQ=255)
                                                                               gACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTc                                                                       >  2:1591726/1‑150 (MQ=255)
                                                                                aCTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCa                                                                      >  2:505260/1‑150 (MQ=255)
                                                                                   cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAat                                                                   <  2:3024525/150‑1 (MQ=255)
                                                                                   cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAat                                                                   <  2:1931949/150‑1 (MQ=255)
                                                                                     caATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAatat                                                                 <  2:3356765/150‑1 (MQ=255)
                                                                                     caATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAatat                                                                 >  1:3368249/1‑150 (MQ=255)
                                                                                      aaTAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAAtata                                                                >  1:3443903/1‑150 (MQ=255)
                                                                                        tAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaa                                                              >  2:2597232/1‑150 (MQ=255)
                                                                                        tAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAATCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaa                                                              >  1:2316194/1‑150 (MQ=255)
                                                                                          gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATaaaaa                                                            >  1:470157/1‑150 (MQ=255)
                                                                                           aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAAc                                                           >  2:3544178/1‑150 (MQ=255)
                                                                                            aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACa                                                          >  1:2520128/1‑150 (MQ=255)
                                                                                               aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTa                                                       >  1:1865057/1‑150 (MQ=255)
                                                                                                ggCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTAc                                                      <  2:3737528/150‑1 (MQ=255)
                                                                                                  ctttatGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGACGAACTGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACAGTAGCGGTTCTAATATAAAAACAGTACCg                                                    <  2:3209442/145‑1 (MQ=255)
                                                                                                  cATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCg                                                    <  2:237709/150‑1 (MQ=255)
                                                                                                         cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacg                                             <  2:337611/150‑1 (MQ=255)
                                                                                                         cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacg                                             <  1:890271/150‑1 (MQ=255)
                                                                                                         cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGGACCGCCCGacg                                             >  1:327453/1‑150 (MQ=255)
                                                                                                           cTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacgga                                           >  2:1821769/1‑150 (MQ=255)
                                                                                                           cTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacgga                                           >  2:1329298/1‑150 (MQ=255)
                                                                                                           cTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacgga                                           >  1:3073704/1‑150 (MQ=255)
                                                                                                               ttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGacggaacgg                                       <  2:123181/150‑1 (MQ=255)
                                                                                                                   gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATcc                                   <  2:3713363/150‑1 (MQ=255)
                                                                                                                   gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATcc                                   <  2:3713770/150‑1 (MQ=255)
                                                                                                                   gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATcc                                   <  2:3990056/150‑1 (MQ=255)
                                                                                                                    cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                  >  2:581512/1‑150 (MQ=255)
                                                                                                                     gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTg                                 <  1:230102/150‑1 (MQ=255)
                                                                                                                      aaGTGATACAACTCTGCTTCCGGGAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTgg                                >  1:3479741/1‑150 (MQ=255)
                                                                                                                      aaGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTgg                                >  1:1852630/1‑150 (MQ=255)
                                                                                                                      aaGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTgg                                <  2:3934208/150‑1 (MQ=255)
                                                                                                                      aaGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACAGCCCGACGGAACGGATCCTgg                                >  2:2681249/1‑150 (MQ=255)
                                                                                                                          gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCttt                            >  1:682278/1‑150 (MQ=255)
                                                                                                                            tACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttg                          >  2:1041107/1‑150 (MQ=255)
                                                                                                                               aaCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgttg                       >  1:1789685/1‑150 (MQ=255)
                                                                                                                                   ctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAct                   <  2:3530404/150‑1 (MQ=255)
                                                                                                                                       ttCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTc               <  1:1221011/150‑1 (MQ=255)
                                                                                                                                       ttCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTc               >  1:3198461/1‑150 (MQ=255)
                                                                                                                                        tCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcg              <  1:3959239/150‑1 (MQ=255)
                                                                                                                                         ccGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcgc             >  1:2059602/1‑150 (MQ=255)
                                                                                                                                           gggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCgg           <  2:3243335/150‑1 (MQ=255)
                                                                                                                                             gAAATAGGTTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGaa         >  2:1334683/1‑150 (MQ=255)
                                                                                                                                               aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACt       <  1:1595146/150‑1 (MQ=255)
                                                                                                                                                 tAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGa     <  1:2198963/150‑1 (MQ=255)
                                                                                                                                                   gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATa   <  1:852526/150‑1 (MQ=255)
                                                                                                                                                   gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATa   >  2:2081239/1‑150 (MQ=255)
                                                                                                                                                    ggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg  <  1:2065137/150‑1 (MQ=255)
                                                                                                                                                    ggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg  <  2:1179385/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTCGCAGACAGCGGTGCATAAGGCGGATCGCAATAGACGACGGATGCATCATCTGCGCGCGCCATGCTATCGGCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAG  >  CP009273/3508703‑3509000

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: ‑