Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,837,997 C→T 5.4% A225V (GCC→GTC)  ade → cryptic adenine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,837,9970CT 5.4% 676.3 / 16.3 186A225V (GCC→GTC) adecryptic adenine deaminase
Reads supporting (aligned to +/- strand):  ref base C (83/93);  new base T (2/8);  total (85/101)
Fisher's exact test for biased strand distribution p-value = 1.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 187 total aligned reads displayed.

ACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAAT  >  CP009273/3837850‑3838143
                                                                                                                                                   |                                                                                                                                                  
    tCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCtat                                                                                                                                                <  1:1637472/148‑1 (MQ=255)
     ccGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCtata                                                                                                                                               <  1:2293346/148‑1 (MQ=255)
      cGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTatat                                                                                                                                              <  1:1212841/148‑1 (MQ=255)
      cGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTatat                                                                                                                                              <  2:2573458/148‑1 (MQ=255)
       gCAGGTTACCGTCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCATAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCGGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATAtt                                                                                                                                             <  2:2566421/148‑1 (MQ=255)
       gCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATAtt                                                                                                                                             <  2:308073/148‑1 (MQ=255)
        cAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTa                                                                                                                                            <  2:3050402/148‑1 (MQ=255)
           gTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTg                                                                                                                                         <  2:2442366/148‑1 (MQ=255)
           gTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTg                                                                                                                                         <  2:1228952/148‑1 (MQ=255)
           gTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTg                                                                                                                                         >  1:1892151/1‑148 (MQ=255)
           gTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTg                                                                                                                                         >  1:3084623/1‑148 (MQ=255)
            ttACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGc                                                                                                                                        >  1:300673/1‑148 (MQ=255)
             tACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCg                                                                                                                                       <  1:1698243/148‑1 (MQ=255)
               ccGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGTCTATATTACTGCggg                                                                                                                                     >  2:1406512/1‑148 (MQ=255)
                cGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGt                                                                                                                                    <  1:561656/148‑1 (MQ=255)
                 ggCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTa                                                                                                                                   <  2:2782941/148‑1 (MQ=255)
                  gCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAt                                                                                                                                  >  2:2472518/1‑148 (MQ=255)
                  gCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAt                                                                                                                                  <  1:2808856/148‑1 (MQ=255)
                     ttGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGa                                                                                                                               >  2:2529654/1‑148 (MQ=255)
                      tGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGaa                                                                                                                              <  2:1960308/148‑1 (MQ=255)
                         aGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAAc                                                                                                                           <  2:1078516/148‑1 (MQ=255)
                              tgatgGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCa                                                                                                                      <  2:2258134/148‑1 (MQ=255)
                                      tACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGTCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                                                                                                              <  1:36388/148‑1 (MQ=255)
                                      tACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGTCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                                                                                                              <  1:1040396/148‑1 (MQ=255)
                                      tACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                                                                                                              <  2:264982/148‑1 (MQ=255)
                                      tACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                                                                                                              >  2:1671709/1‑148 (MQ=255)
                                      tACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTCAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGtt                                                                                                              <  1:2869071/148‑1 (MQ=255)
                                         ccTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATc                                                                                                           <  1:792900/148‑1 (MQ=255)
                                           tGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAg                                                                                                         >  2:1220864/1‑148 (MQ=255)
                                             gCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCt                                                                                                       >  1:2287466/1‑148 (MQ=255)
                                               gTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTgg                                                                                                     <  1:1737173/148‑1 (MQ=255)
                                               gTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTgg                                                                                                     <  1:3242295/148‑1 (MQ=255)
                                                tAATTAGCGGGCATAATGCTCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAATAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGga                                                                                                    <  1:843195/148‑1 (MQ=255)
                                                tAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGga                                                                                                    >  1:3078301/1‑148 (MQ=255)
                                                tAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGga                                                                                                    <  2:1474687/148‑1 (MQ=255)
                                                   ttAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAtt                                                                                                                                 <  1:174105/116‑1 (MQ=255)
                                                   ttAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAtt                                                                                                                                 <  1:176446/116‑1 (MQ=255)
                                                   ttAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTAtt                                                                                                                                 >  2:176446/1‑116 (MQ=255)
                                                   ttAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaaga                                                                                                 >  1:537101/1‑148 (MQ=255)
                                                    tAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaagaa                                                                                                >  2:2256289/1‑148 (MQ=255)
                                                    tAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGgaagaa                                                                                                <  1:2384829/148‑1 (MQ=255)
                                                      gCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAgg                                                                                              <  1:177824/148‑1 (MQ=255)
                                                      gCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAgg                                                                                              <  2:1841166/148‑1 (MQ=255)
                                                       cGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGa                                                                                             >  1:1984425/1‑148 (MQ=255)
                                                       cGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGa                                                                                             >  1:2828025/1‑148 (MQ=255)
                                                         ggCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACg                                                                                           <  2:2537680/148‑1 (MQ=255)
                                                           cAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGcc                                                                                         >  1:218768/1‑148 (MQ=255)
                                                                tGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGaaa                                                                                    >  2:2004694/1‑148 (MQ=255)
                                                                 gcgcTGCTCTATAAACTGGATGCATTTCGCCACCTGACGATGTACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGTCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAAt                                                                                   <  1:1792553/148‑1 (MQ=255)
                                                                     tgctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGTCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACa                                                                               <  1:3143641/148‑1 (MQ=255)
                                                                     tgctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACa                                                                               <  2:469177/148‑1 (MQ=255)
                                                                       ctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAAc                                                                             >  2:3236349/1‑148 (MQ=255)
                                                                       ctCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAAc                                                                             >  2:1190993/1‑148 (MQ=255)
                                                                        tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACt                                                                            <  1:248647/148‑1 (MQ=255)
                                                                        tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACt                                                                            <  1:1917571/148‑1 (MQ=255)
                                                                        tCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACt                                                                            <  2:1861245/148‑1 (MQ=255)
                                                                            tAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGc                                                                        >  2:3216763/1‑148 (MQ=255)
                                                                            tAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGc                                                                        >  1:1241989/1‑148 (MQ=255)
                                                                            tAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCACGGGTTTGGGTGGTAAATAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGc                                                                        >  2:883908/1‑148 (MQ=255)
                                                                             aaaCTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCa                                                                       >  1:992470/1‑148 (MQ=255)
                                                                              aaCTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCAt                                                                      >  2:2958486/1‑148 (MQ=255)
                                                                               aCTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATg                                                                     >  1:3144862/1‑148 (MQ=255)
                                                                                  ggATGCATGTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCg                                                                  <  2:2632168/148‑1 (MQ=255)
                                                                                   gATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGt                                                                 <  2:3271951/148‑1 (MQ=255)
                                                                                         tttCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtgatga                                                           >  1:734303/1‑148 (MQ=255)
                                                                                         tttCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATACTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTtgatga                                                           >  1:2552392/1‑148 (MQ=255)
                                                                                                  cTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAg                                                  <  2:2186576/148‑1 (MQ=255)
                                                                                                  cTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAg                                                  >  1:2899796/1‑148 (MQ=255)
                                                                                                     aCGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTCACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                                               <  1:909589/148‑1 (MQ=255)
                                                                                                     aCGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGt                                               <  2:240450/148‑1 (MQ=255)
                                                                                                          ggACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCt                                          >  2:2256572/1‑148 (MQ=255)
                                                                                                           gACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTg                                         >  2:3063074/1‑148 (MQ=255)
                                                                                                            aCGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGc                                        >  1:2091914/1‑148 (MQ=255)
                                                                                                            aCGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGc                                        >  1:151089/1‑148 (MQ=255)
                                                                                                                tCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGc                                    >  1:1130321/1‑148 (MQ=255)
                                                                                                                tCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGc                                    >  2:1838006/1‑148 (MQ=255)
                                                                                                                 cACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCa                                   >  2:488999/1‑148 (MQ=255)
                                                                                                                 cACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCAACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTAGTTGATGATCCGCGAAGGGTCCGCTGCCCGCa                                   >  2:2226564/1‑148 (MQ=255)
                                                                                                                  aCTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCaa                                  <  2:597883/148‑1 (MQ=255)
                                                                                                                    tGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAAtc                                <  2:1610832/148‑1 (MQ=255)
                                                                                                                         gggTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAAc                           >  2:205408/1‑148 (MQ=255)
                                                                                                                          ggTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAAcg                          <  1:1872526/148‑1 (MQ=255)
                                                                                                                             ttGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACgcgc                       >  2:409812/1‑148 (MQ=255)
                                                                                                                                ggtggtAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTgg                    >  1:2215540/1‑148 (MQ=255)
                                                                                                                                      aaaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGt              >  1:2668954/1‑148 (MQ=255)
                                                                                                                                      aaaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGt              >  1:1015734/1‑148 (MQ=255)
                                                                                                                                       aaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGtt             >  1:3119492/1‑148 (MQ=255)
                                                                                                                                       aaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGtt             >  2:1081172/1‑148 (MQ=255)
                                                                                                                                       aaGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGtt             >  2:2121209/1‑148 (MQ=255)
                                                                                                                                        aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg            <  2:3291453/148‑1 (MQ=255)
                                                                                                                                        aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg            >  1:435847/1‑148 (MQ=255)
                                                                                                                                        aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg            <  2:2822784/148‑1 (MQ=255)
                                                                                                                                        aGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTg            >  2:1306167/1‑148 (MQ=255)
                                                                                                                                         gAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGa           <  2:2998077/148‑1 (MQ=255)
                                                                                                                                          aaCTTAACGTCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAt          <  1:1466282/148‑1 (MQ=255)
                                                                                                                                          aaCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGAt          <  2:533993/148‑1 (MQ=255)
                                                                                                                                           aCTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATc         >  1:3108237/1‑148 (MQ=255)
                                                                                                                                             ttAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCaa       >  1:1296818/1‑148 (MQ=255)
                                                                                                                                                 cGCCTATATTACTGCGGGTATTGAAAACTTCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGaa   <  2:3133244/148‑1 (MQ=255)
                                                                                                                                                 cGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGaa   <  1:2256303/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACCATCCGCAGGTTACCGGCCTTGCAGAAATGATGGACTACCCTGGCGTAATTAGCGGGCAGAATGCGCTGCTCGATAAACTGGATGCATTTCGCCACCTGACGCTGGACGGTCACTGCCCGGGTTTGGGTGGTAAAGAACTTAACGCCTATATTACTGCGGGTATTGAAAACTGCCACGAAAGTTATCAGCTGGAAGAAGGACGCCGGAAATTACAACTCGGCATGTCGTTGATGATCCGCGAAGGGTCCGCTGCCCGCAATCTCAACGCGCTGGCACCGTTGATCAACGAAT  >  CP009273/3837850‑3838143

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: