Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,615,947 G→T 100% E73D (GAG→GAT)  radA → DNA repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,615,9470GT97.2% 553.8 / 5.1 178E73D (GAG→GAT) radADNA repair protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base T (81/92);  total (85/93)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 181 total aligned reads displayed.

AGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAG  >  CP009273/4615804‑4616094
                                                                                                                                               |                                                                                                                                                   
agtgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGc                                                                                                                                                 <  2:2311017/148‑1 (MQ=255)
  tggctttCATGCCTGGAACACCATCACCGAGGTGCGTTTTGCTGCGTCGCCAATGGTGGCGGGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcg                                                                                                                                               <  2:299684/142‑1 (MQ=255)
  tgcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcg                                                                                                                                               <  2:2096684/148‑1 (MQ=255)
  ggcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcg                                                                                                                                               <  2:2096386/147‑1 (MQ=255)
   gcCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCcgc                                                                                                                                              >  1:1732698/1‑148 (MQ=255)
     cTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCCGCGGCTAGGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGt                                                                                                                                            <  1:995658/148‑1 (MQ=255)
     cTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGt                                                                                                                                            <  1:895637/148‑1 (MQ=255)
       gTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGttt                                                                                                                                          <  1:3165050/148‑1 (MQ=255)
        tCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGtttt                                                                                                                                         <  2:1108136/148‑1 (MQ=255)
          aggCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTc                                                                                                                                       <  2:182153/146‑1 (MQ=255)
           tGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTcc                                                                                                                                      <  2:1673556/148‑1 (MQ=255)
            gCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCa                                                                                                                                     <  2:2375548/148‑1 (MQ=255)
            gCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCa                                                                                                                                     <  2:2562287/148‑1 (MQ=255)
                 gAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGCGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGa                                                                                                                                <  1:3300741/148‑1 (MQ=255)
                 gAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCNGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGa                                                                                                                                >  1:1806206/1‑148 (MQ=255)
                 gAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGa                                                                                                                                >  1:1478265/1‑148 (MQ=255)
                 gAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGa                                                                                                                                >  2:1307310/1‑148 (MQ=255)
                  aaCACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGAt                                                                                                                               >  1:933876/1‑148 (MQ=255)
                      ccATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                           <  1:2816367/148‑1 (MQ=255)
                      ccATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                           <  1:548483/148‑1 (MQ=255)
                      ccATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                           <  1:3526848/148‑1 (MQ=255)
                      ccATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTaa                                                                                                                           >  2:3094071/1‑148 (MQ=255)
                        aTCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAag                                                                                                                         >  2:1218486/1‑148 (MQ=255)
                        aTCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAag                                                                                                                         <  1:1501997/148‑1 (MQ=255)
                        aTCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAag                                                                                                                         <  2:3308216/148‑1 (MQ=255)
                          ccccGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAagag                                                                                                                       <  1:2733541/146‑1 (MQ=255)
                          cACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAagag                                                                                                                       >  1:2029095/1‑148 (MQ=255)
                           aCCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                      <  2:1913306/148‑1 (MQ=255)
                           aCCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGt                                                                                                                      <  2:3425717/148‑1 (MQ=255)
                            ccGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCGGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGtt                                                                                                                     <  2:870905/148‑1 (MQ=255)
                            ccGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGtt                                                                                                                     <  2:2412360/148‑1 (MQ=255)
                              gAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCg                                                                                                                   <  1:2499927/148‑1 (MQ=255)
                              gAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCg                                                                                                                   <  1:2967458/148‑1 (MQ=255)
                              gAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCg                                                                                                                   <  2:3366516/148‑1 (MQ=255)
                                ggTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCTATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGAc                                                                                                                 >  1:153226/1‑148 (MQ=255)
                                      tCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTa                                                                                                           >  1:197874/1‑148 (MQ=255)
                                       cTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTAc                                                                                                          >  2:1459331/1‑148 (MQ=255)
                                           ctgcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAgg                                                                                                      >  2:1316860/1‑148 (MQ=255)
                                             gcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcg                                                                                                    <  2:1347140/148‑1 (MQ=255)
                                             gcGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcg                                                                                                    <  2:3054623/148‑1 (MQ=255)
                                              cGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcgg                                                                                                   >  1:564957/1‑148 (MQ=255)
                                              cGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcgg                                                                                                   >  1:2911963/1‑148 (MQ=255)
                                              cGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcgg                                                                                                   >  1:3206756/1‑148 (MQ=255)
                                                tCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcg                                                                                                 >  2:337133/1‑148 (MQ=255)
                                                 cGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAggcggcgg                                                                                                >  2:775457/1‑148 (MQ=255)
                                                     aaTGGTGGCGCGTACCGAGCGTCTCCGCGGCTAGGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtg                                                                                            <  1:438773/148‑1 (MQ=255)
                                                     aaTGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtg                                                                                            >  1:3043939/1‑148 (MQ=255)
                                                     aaTGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtg                                                                                            >  1:3043883/1‑148 (MQ=255)
                                                      aTGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCgtgg                                                                                           >  2:529757/1‑148 (MQ=255)
                                                         gtggCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGc                                                                                        <  2:3372518/148‑1 (MQ=255)
                                                            gcgcGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAg                                                                                     >  2:2795694/1‑148 (MQ=255)
                                                               cgTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGaa                                                                                  <  2:3226186/148‑1 (MQ=255)
                                                               cgTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGaa                                                                                  >  1:1006862/1‑148 (MQ=255)
                                                                gTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAg                                                                                 <  1:1307310/148‑1 (MQ=255)
                                                                  aaCGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTg                                                                               <  2:1358628/148‑1 (MQ=255)
                                                                    cGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGcc                                                                             <  2:1653917/148‑1 (MQ=255)
                                                                       gCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCAtt                                                                          <  2:1841037/148‑1 (MQ=255)
                                                                        cGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTc                                                                         >  2:507751/1‑148 (MQ=255)
                                                                        cGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTc                                                                         >  2:364743/1‑148 (MQ=255)
                                                                          tctcCGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTg                                                                       <  1:1545038/148‑1 (MQ=255)
                                                                           ctcAGCTGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                      >  2:1405732/1‑148 (MQ=255)
                                                                           ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                      >  2:1405711/1‑148 (MQ=255)
                                                                           ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGa                                                                      >  1:1190511/1‑148 (MQ=255)
                                                                           ctcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGACTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCTGCGGCGTGGTGACAGGAAGTGCCATTCTGa                                                                      >  2:2722665/1‑148 (MQ=255)
                                                                            tcAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGAt                                                                     >  1:1181255/1‑148 (MQ=255)
                                                                              aGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTg                                                                   >  2:480248/1‑148 (MQ=255)
                                                                               gCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTACAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTgg                                                                  >  2:2809751/1‑148 (MQ=255)
                                                                               gCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTgg                                                                  <  1:2791905/148‑1 (MQ=255)
                                                                                cGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGc                                                                 <  1:2307645/148‑1 (MQ=255)
                                                                                     aTGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCAGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTaa                                                            <  2:146722/148‑1 (MQ=255)
                                                                                       gCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc                                                          >  1:2702338/1‑148 (MQ=255)
                                                                                       gCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc                                                          <  2:2641648/148‑1 (MQ=255)
                                                                                       gCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAAcc                                                          >  1:2702262/1‑148 (MQ=255)
                                                                                            tAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGt                                                     >  1:3344119/1‑148 (MQ=255)
                                                                                               cgcCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCg                                                  >  1:1123068/1‑148 (MQ=255)
                                                                                               cgcCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCg                                                  >  2:2865612/1‑148 (MQ=255)
                                                                                                    gggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaa                                             >  2:2417674/1‑148 (MQ=255)
                                                                                                    gggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaa                                             >  1:2647336/1‑148 (MQ=255)
                                                                                                     ggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaaa                                            >  2:2111646/1‑148 (MQ=255)
                                                                                                     ggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaaa                                            <  1:2951823/148‑1 (MQ=255)
                                                                                                     ggTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGaaa                                            <  1:3243451/148‑1 (MQ=255)
                                                                                                        ggCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATcc                                         <  2:1806206/148‑1 (MQ=255)
                                                                                                        ggCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATcc                                         <  1:1459331/148‑1 (MQ=255)
                                                                                                         gCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCa                                        <  1:337133/148‑1 (MQ=255)
                                                                                                            aaaGTCCAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgc                                     >  1:368256/1‑148 (MQ=255)
                                                                                                                tCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACgctgct                                 <  1:241636/148‑1 (MQ=255)
                                                                                                                 ccAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTa                                <  2:2728419/148‑1 (MQ=255)
                                                                                                                 ccAGAAACTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTa                                <  2:3505726/148‑1 (MQ=255)
                                                                                                                     aaaCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGc                            >  1:2200555/1‑148 (MQ=255)
                                                                                                                       aCTCTCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCaa                          <  2:173732/148‑1 (MQ=255)
                                                                                                                        ctctCCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCaaa                         <  2:2172680/148‑1 (MQ=255)
                                                                                                                          ctcCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACg                       >  1:50872/1‑148 (MQ=255)
                                                                                                                          ctcCGATATCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACg                       <  2:1028663/148‑1 (MQ=255)
                                                                                                                               atatCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCtgtg                  >  2:1062474/1‑148 (MQ=255)
                                                                                                                                  tCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaa               >  2:1514327/1‑148 (MQ=255)
                                                                                                                                  tCAGCCTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCaa               <  1:196531/148‑1 (MQ=255)
                                                                                                                                      ccTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTg           >  2:2046407/1‑148 (MQ=255)
                                                                                                                                      ccTTGAAGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTg           >  1:3153003/1‑148 (MQ=255)
                                                                                                                                           aaGATCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCca      >  2:801653/1‑148 (MQ=255)
                                                                                                                                            agatCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCcag     >  1:2217710/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
AGTGCCTGTCATGCCTGGAACACCATCACCGAGGTGCGTCTTGCTGCGTCGCCAATGGTGGCGCGTAACGAGCGTCTCAGCGGCTATGCCGGTAGCGCCGGGGTGGCAAAAGTCCAGAAACTCTCCGATATCAGCCTTGAAGAGCTGCCGCGTTTTTCCACCGGATTTAAAGAGTTCGACCGCGTACTAGGCGGCGGCGTGGTGCCAGGAAGTGCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAG  >  CP009273/4615804‑4616094

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