Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,019,203 G→T 33.1% E27* (GAG→TAG)  guaD → guanine deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,019,2030GT33.1% 432.0 / 241.9 254E27* (GAG→TAG) guaDguanine deaminase
Reads supporting (aligned to +/- strand):  ref base G (96/74);  new base T (41/43);  total (137/117)
Fisher's exact test for biased strand distribution p-value = 2.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 257 total aligned reads displayed.

CCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCG  >  CP009273/3019059‑3019352
                                                                                                                                                |                                                                                                                                                     
ccGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGAtt                                                                                                                                                  <  2:752948/150‑1 (MQ=255)
ccGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGAtt                                                                                                                                                  <  2:4520454/150‑1 (MQ=255)
 cGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTg                                                                                                                                                 <  1:3360846/150‑1 (MQ=255)
     aGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATCACCCGGAATAGATTGCctc                                                                                                                                             >  2:479735/1‑150 (MQ=255)
     aGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCctc                                                                                                                                             >  2:479600/1‑150 (MQ=255)
     aGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCctc                                                                                                                                             >  2:193633/1‑150 (MQ=255)
      ggAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCctct                                                                                                                                            >  1:2847094/1‑150 (MQ=255)
         aaaCGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcg                                                                                                                                         >  1:4772155/1‑150 (MQ=255)
         aaaCGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTctgcg                                                                                                                                         <  2:671594/150‑1 (MQ=255)
          aaCGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTctgcgc                                                                                                                                        >  2:3370889/1‑150 (MQ=255)
                 tGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGtt                                                                                                                                 <  2:3590736/150‑1 (MQ=255)
                    cTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTAt                                                                                                                              >  1:1278862/1‑150 (MQ=255)
                     tGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCACTGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTAtt                                                                                                                             <  2:260259/150‑1 (MQ=255)
                     tGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTAtt                                                                                                                             <  2:3627982/150‑1 (MQ=255)
                     tGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTAtt                                                                                                                             <  2:4145683/150‑1 (MQ=255)
                       gTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGa                                                                                                                           >  2:237763/1‑150 (MQ=255)
                        tATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAg                                                                                                                          <  2:1311278/150‑1 (MQ=255)
                        tATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAg                                                                                                                          <  2:4687008/150‑1 (MQ=255)
                           tACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGAt                                                                                                                       <  2:1014697/150‑1 (MQ=255)
                           tACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGAt                                                                                                                       <  2:1000790/150‑1 (MQ=255)
                            aCCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATg                                                                                                                      <  2:4767479/150‑1 (MQ=255)
                             ccTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATgg                                                                                                                     <  2:1159203/150‑1 (MQ=255)
                               tCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGtt                                                                                                                   <  1:4455063/150‑1 (MQ=255)
                                   cGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatt                                                                                                               <  1:1024697/150‑1 (MQ=255)
                                   cGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatt                                                                                                               >  2:530889/1‑150 (MQ=255)
                                   cGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatt                                                                                                               >  2:2158367/1‑150 (MQ=255)
                                    ggAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTttatta                                                                                                              <  1:2187292/150‑1 (MQ=255)
                                     gAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTAc                                                                                                             <  2:2494124/150‑1 (MQ=255)
                                      aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACt                                                                                                            <  2:546029/150‑1 (MQ=255)
                                      aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACt                                                                                                            >  1:1720228/1‑150 (MQ=255)
                                      aaGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACt                                                                                                            >  1:2221408/1‑150 (MQ=255)
                                       agagATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTc                                                                                                           >  1:4362516/1‑150 (MQ=255)
                                       agagATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTc                                                                                                           >  1:4696384/1‑150 (MQ=255)
                                        gagaTGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCa                                                                                                          <  2:4409256/150‑1 (MQ=255)
                                         agaTGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCAt                                                                                                         <  1:234157/150‑1 (MQ=255)
                                             ggATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTaaa                                                                                                     >  1:1866747/1‑150 (MQ=255)
                                                    cAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGaa                                                                                              <  1:375202/150‑1 (MQ=255)
                                                     aGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGaaa                                                                                             >  1:3870593/1‑150 (MQ=255)
                                                       tAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAg                                                                                           <  2:979992/150‑1 (MQ=255)
                                                         aaGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTg                                                                                         >  1:1160092/1‑150 (MQ=255)
                                                         aaGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTg                                                                                         >  1:2131003/1‑150 (MQ=255)
                                                             aGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAAt                                                                                     >  2:3492144/1‑150 (MQ=255)
                                                             aGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAAt                                                                                     <  2:1710696/150‑1 (MQ=255)
                                                              gTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTGATTGAGGATGGTGTATTACTCATTAAACAGGGAAAAGTGGAATg                                                                                    >  2:4111875/1‑150 (MQ=255)
                                                                 aaTGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGtt                                                                                 >  1:4121599/1‑150 (MQ=255)
                                                                  atgatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGttt                                                                                >  1:596428/1‑150 (MQ=255)
                                                                   tgatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTg                                                                               >  2:909811/1‑150 (MQ=255)
                                                                    gatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTgg                                                                              <  2:3710581/150‑1 (MQ=255)
                                                                    gatgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTgg                                                                              <  2:4188953/150‑1 (MQ=255)
                                                                     atgTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGc                                                                             >  2:283308/1‑150 (MQ=255)
                                                                       gTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGa                                                                           >  2:1576462/1‑150 (MQ=255)
                                                                          aGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATg                                                                        >  1:3031563/1‑150 (MQ=255)
                                                                           ggAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATgg                                                                       >  1:1087206/1‑150 (MQ=255)
                                                                            gagaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATggg                                                                      >  2:2030567/1‑150 (MQ=255)
                                                                             agaACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGa                                                                     >  2:4349689/1‑150 (MQ=255)
                                                                                acacacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCCGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGaaaa                                                                  >  1:2662183/1‑150 (MQ=255)
                                                                                 cacacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAAc                                                                 <  1:1015809/150‑1 (MQ=255)
                                                                                  acacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACg                                                                <  2:2985758/150‑1 (MQ=255)
                                                                                  acacGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACg                                                                <  2:3980372/150‑1 (MQ=255)
                                                                                    acGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGa                                                              <  1:3268825/150‑1 (MQ=255)
                                                                                      gTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGaaa                                                            >  2:2544203/1‑150 (MQ=255)
                                                                                       ttAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAg                                                           <  1:193633/150‑1 (MQ=255)
                                                                                           aaGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATc                                                       >  1:2742258/1‑150 (MQ=255)
                                                                                            aGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCa                                                      <  1:788170/150‑1 (MQ=255)
                                                                                              tggTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGACAGCATCaaa                                                    >  2:2039219/2‑150 (MQ=255)
                                                                                                   cGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCt                                               >  2:3302736/1‑150 (MQ=255)
                                                                                                   cGAGGCAGTTTTATTGATGTCACCAGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCt                                               <  2:4559233/150‑1 (MQ=255)
                                                                                                    gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                              >  1:1589380/1‑150 (MQ=255)
                                                                                                    gAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTg                                              >  1:1867054/1‑150 (MQ=255)
                                                                                                     aGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGa                                             <  1:4594258/150‑1 (MQ=255)
                                                                                                      ggCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGac                                            >  2:2207486/1‑150 (MQ=255)
                                                                                                        cAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGacac                                          <  1:3370889/150‑1 (MQ=255)
                                                                                                         aGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACAcc                                         <  2:3031563/150‑1 (MQ=255)
                                                                                                             ttATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTc                                     <  2:1261820/150‑1 (MQ=255)
                                                                                                             ttATTGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTc                                     <  2:1260897/150‑1 (MQ=255)
                                                                                                             ttATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTc                                     <  2:3004531/150‑1 (MQ=255)
                                                                                                              tATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTcg                                    <  1:3311111/150‑1 (MQ=255)
                                                                                                                 tGATGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGt                                 >  2:3106108/1‑150 (MQ=255)
                                                                                                                 tGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAACGCATCAAATTCCTGACACCATTCGCGt                                 >  2:1097698/1‑150 (MQ=255)
                                                                                                                    tGTCACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTgcg                              >  2:4542631/1‑150 (MQ=255)
                                                                                                                       cACCCGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGa                           >  1:3472279/1‑150 (MQ=255)
                                                                                                                       cACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGa                           <  2:3462388/150‑1 (MQ=255)
                                                                                                                       cACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGa                           <  2:1518413/150‑1 (MQ=255)
                                                                                                                         cccGTACGATCGATAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACt                         <  1:1543279/150‑1 (MQ=255)
                                                                                                                         cccGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACt                         >  2:2536426/1‑150 (MQ=255)
                                                                                                                          ccGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGACAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTa                        >  2:2470557/1‑150 (MQ=255)
                                                                                                                              aCGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgc                    <  1:2030567/150‑1 (MQ=255)
                                                                                                                              aCGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgc                    <  1:4435471/150‑1 (MQ=255)
                                                                                                                              aCGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgc                    >  2:2904472/1‑150 (MQ=255)
                                                                                                                               cgatcgatAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATcgcg                   <  2:4362516/150‑1 (MQ=255)
                                                                                                                                    gatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaaa              >  1:4159713/1‑150 (MQ=255)
                                                                                                                                    gatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaaa              >  1:4469818/1‑150 (MQ=255)
                                                                                                                                    gatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaaa              >  1:2467494/1‑150 (MQ=255)
                                                                                                                                    gatAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCaaa              <  2:860078/150‑1 (MQ=255)
                                                                                                                                     atAACCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAAc             >  1:3087826/1‑150 (MQ=255)
                                                                                                                                       aaCCCGGAATAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTg           >  1:977381/1‑150 (MQ=255)
                                                                                                                                       aaCCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTg           >  1:1360448/1‑150 (MQ=255)
                                                                                                                                         cccGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGAt         <  1:2153706/150‑1 (MQ=255)
                                                                                                                                           cGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAg       <  2:2063850/150‑1 (MQ=255)
                                                                                                                                           cGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAg       <  1:2468907/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
CCGACAGGAAAACGTTCTGCCTGGTATTACCTCAGCGGAAGAGATGGATTAACAGTAAAGGAGTCAATGATGTCAGGAGAACACACGTTAAAAGCGGTACGAGGCAGTTTTATTGATGTCACCCGTACGATCGATAACCCGGAAGAGATTGCCTCTGCGCTGCGGTTTATTGAGGATGGTTTATTACTCATTAAACAGGGAAAAGTGGAATGGTTTGGCGAATGGGAAAACGGAAAGCATCAAATTCCTGACACCATTCGCGTGCGCGACTATCGCGGCAAACTGATAGTACCG  >  CP009273/3019059‑3019352

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