Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,522,558 G→T 5.2% A177A (GCG→GCT yrfG → GMP/IMP nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,522,5580GT 5.2% 1073.6 / 26.1 278A177A (GCG→GCTyrfGGMP/IMP nucleotidase
Reads supporting (aligned to +/- strand):  ref base G (112/151);  new base T (9/6);  total (121/157)
Fisher's exact test for biased strand distribution p-value = 2.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

GAGCATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACGCCA  >  CP009273/3522409‑3522707
                                                                                                                                                     |                                                                                                                                                     
  gCATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGAt                                                                                                                                                     <  2:3465289/150‑1 (MQ=255)
   cATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGAtt                                                                                                                                                    <  1:2798054/150‑1 (MQ=255)
      aCCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACTCATTTTGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCTAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATtcac                                                                                                                                                 <  2:1373127/150‑3 (MQ=255)
      aCCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATtctc                                                                                                                                                 >  2:1364956/1‑150 (MQ=255)
      aCCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATtctc                                                                                                                                                 <  1:4974214/150‑1 (MQ=255)
      aCCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATtctc                                                                                                                                                 <  1:4875404/150‑1 (MQ=255)
      aCCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATtctc                                                                                                                                                 <  1:2381113/150‑1 (MQ=255)
       ccGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCg                                                                                                                                                <  2:1158881/150‑1 (MQ=255)
        cGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGa                                                                                                                                               <  2:850154/150‑1 (MQ=255)
        cGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGa                                                                                                                                               >  1:2565565/1‑150 (MQ=255)
         ggTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGAt                                                                                                                                              >  2:2542518/1‑150 (MQ=255)
          gTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGAtg                                                                                                                                             >  2:2823027/1‑150 (MQ=255)
            cTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGAtgct                                                                                                                                           <  1:2065746/150‑1 (MQ=255)
                aCGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCcg                                                                                                                                       <  1:5131471/150‑1 (MQ=255)
                     cacCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCaa                                                                                                                                  >  1:1341802/1‑150 (MQ=255)
                        cTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAAttt                                                                                                                               <  2:639438/150‑1 (MQ=255)
                        cTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAAttt                                                                                                                               <  2:3265915/150‑1 (MQ=255)
                        cTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAAttt                                                                                                                               >  1:822830/1‑150 (MQ=255)
                         ttGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTg                                                                                                                              <  2:3973428/150‑1 (MQ=255)
                         ttGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTg                                                                                                                              <  1:3195655/150‑1 (MQ=255)
                         ttGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTg                                                                                                                              <  1:565173/150‑1 (MQ=255)
                            aTTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTa                                                                                                                           <  1:1169101/150‑1 (MQ=255)
                            aTTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTa                                                                                                                           >  1:883255/1‑150 (MQ=255)
                             tttATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTAt                                                                                                                          <  1:4820242/150‑1 (MQ=255)
                              ttattaCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTAtt                                                                                                                         <  1:4672478/150‑1 (MQ=255)
                               tattaCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTc                                                                                                                        <  1:5023578/150‑1 (MQ=255)
                                   aCTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTa                                                                                                                    >  1:631926/1‑150 (MQ=255)
                                        ccacccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGcc                                                                                                               <  1:3902403/150‑1 (MQ=255)
                                         cacccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCt                                                                                                              <  2:2142705/150‑1 (MQ=255)
                                          acccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGCTTCTCGATGCTGCCGCGCAAGTTGGTATTCGTTACTGCCTc                                                                                                             >  2:4504696/1‑150 (MQ=255)
                                          acccacACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTc                                                                                                             >  1:1080458/1‑150 (MQ=255)
                                              acacaTTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCg                                                                                                         <  1:1387043/150‑1 (MQ=255)
                                               cacaTTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGt                                                                                                        >  2:3903646/1‑150 (MQ=255)
                                               cacaTTTGGTTACCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGt                                                                                                        >  2:4424076/1‑150 (MQ=255)
                                                acaTTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTg                                                                                                       <  1:65723/150‑1 (MQ=255)
                                                  aTTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGAc                                                                                                     >  2:5217770/1‑150 (MQ=255)
                                                   tttGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACt                                                                                                    >  2:2601747/1‑150 (MQ=255)
                                                    ttGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTa                                                                                                   <  2:2593577/150‑1 (MQ=255)
                                                     tGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTaa                                                                                                  <  1:3962770/150‑1 (MQ=255)
                                                     tGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTaa                                                                                                  <  1:1798047/150‑1 (MQ=255)
                                                     tGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGATTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTaa                                                                                                  <  1:2108575/150‑1 (MQ=255)
                                                      ggTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAAt                                                                                                 <  2:1627659/150‑1 (MQ=255)
                                                      ggTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAAt                                                                                                 <  2:2393158/150‑1 (MQ=255)
                                                      ggTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAAt                                                                                                 <  2:4277860/150‑1 (MQ=255)
                                                       gTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGCTGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATc                                                                                                <  2:2612428/150‑1 (MQ=255)
                                                          aTCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTg                                                                                             <  1:47457/150‑1 (MQ=255)
                                                          aTCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTg                                                                                             <  1:191902/150‑1 (MQ=255)
                                                           tCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGa                                                                                            <  1:1560454/150‑1 (MQ=255)
                                                            ccGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGAt                                                                                           >  1:1020062/1‑150 (MQ=255)
                                                              gAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTc                                                                                         >  2:486831/1‑150 (MQ=255)
                                                               aaaGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTcc                                                                                        >  2:3304698/1‑150 (MQ=255)
                                                               aaaGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTcc                                                                                        <  1:2038307/150‑1 (MQ=255)
                                                                aaGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCAGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGGATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCg                                                                                       <  2:3125521/150‑1 (MQ=255)
                                                                  gagGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCggg                                                                                     >  2:5187352/1‑150 (MQ=255)
                                                                    ggATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGAt                                                                                   >  2:3822100/1‑150 (MQ=255)
                                                                       tCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGc                                                                                >  1:1882147/1‑150 (MQ=255)
                                                                        cAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGcc                                                                               >  2:943721/1‑150 (MQ=255)
                                                                        cAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGcc                                                                               <  2:2658669/150‑1 (MQ=255)
                                                                         aGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCg                                                                              <  1:3638206/150‑1 (MQ=255)
                                                                           cGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCgag                                                                            >  1:1398926/1‑150 (MQ=255)
                                                                           cGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCgag                                                                            >  1:4480314/1‑150 (MQ=255)
                                                                               tATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTCTTCGTTACTGCCTCGGCGTGACTCATCCTGATTCCGGGATTGCCGAGAAAc                                                                        >  2:1837310/1‑150 (MQ=255)
                                                                                     aTGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATc                                                                  >  2:736594/1‑150 (MQ=255)
                                                                                      tGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCa                                                                 <  2:2240004/150‑1 (MQ=255)
                                                                                      tGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCa                                                                 <  2:5151976/150‑1 (MQ=255)
                                                                                      tGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCa                                                                 <  2:2961089/150‑1 (MQ=255)
                                                                                           tGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGcc                                                            <  1:1271651/150‑1 (MQ=255)
                                                                                           tGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGcc                                                            <  1:4543724/150‑1 (MQ=255)
                                                                                            ggCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCa                                                           <  2:2074481/150‑1 (MQ=255)
                                                                                               cGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGTCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATcc                                                        <  2:4863139/150‑1 (MQ=255)
                                                                                               cGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATcc                                                        >  1:4307563/1‑150 (MQ=255)
                                                                                                      aCGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTg                                                 <  1:326003/150‑1 (MQ=255)
                                                                                                         ggTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAAt                                              >  1:3474805/1‑150 (MQ=255)
                                                                                                                aagctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAcc                                       <  2:3908297/150‑1 (MQ=255)
                                                                                                                 agctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAccg                                      >  2:1955847/1‑150 (MQ=255)
                                                                                                                   ctgaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTAccgcc                                    <  2:3559239/150‑1 (MQ=255)
                                                                                                                     gaaagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACcgccgc                                  >  1:2832934/1‑150 (MQ=255)
                                                                                                                       aagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCt                                <  2:4715470/150‑1 (MQ=255)
                                                                                                                       aagAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCt                                <  2:4213214/150‑1 (MQ=255)
                                                                                                                         gAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGa                              <  1:910612/150‑1 (MQ=255)
                                                                                                                          aaCGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGAt                             >  2:1897661/1‑150 (MQ=255)
                                                                                                                           aCGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATc                            >  1:3883433/1‑150 (MQ=255)
                                                                                                                             gCTGTTTATTGATGACAGCGAAGCTATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATccc                          <  1:3777755/150‑1 (MQ=255)
                                                                                                                                gTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTc                       <  2:4656866/150‑1 (MQ=255)
                                                                                                                                 tttATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCg                      <  2:1504601/150‑1 (MQ=255)
                                                                                                                                 tttATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCg                      >  1:3017378/1‑150 (MQ=255)
                                                                                                                                  ttATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGc                     >  2:4094922/1‑150 (MQ=255)
                                                                                                                                         tgaCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAAtgtg              >  1:1006245/1‑150 (MQ=255)
                                                                                                                                         tgaCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAAtgtg              <  2:4350186/150‑1 (MQ=255)
                                                                                                                                          gaCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGa             <  2:594942/150‑1 (MQ=255)
                                                                                                                                           aCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGaa            >  1:2640433/1‑150 (MQ=255)
                                                                                                                                           aCAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGaa            <  1:4139441/150‑1 (MQ=255)
                                                                                                                                            cAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAg           >  1:402991/1‑150 (MQ=255)
                                                                                                                                             agcgaagctaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTCATGTGAAgg          >  2:3629501/1‑150 (MQ=255)
                                                                                                                                               cgaagcgaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGgag        >  2:2355972/1‑150 (MQ=255)
                                                                                                                                                gaagcgaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGgaga       >  1:1895247/1‑150 (MQ=255)
                                                                                                                                                 aagcgaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTGACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGAc      <  2:1010409/150‑1 (MQ=255)
                                                                                                                                                  agcgaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACg     >  2:3717518/1‑150 (MQ=255)
                                                                                                                                                     gaTTCTCGATGCTGCCGCTCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACGCCa  <  2:2779835/150‑1 (MQ=255)
                                                                                                                                                     gaTTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACGCCa  >  2:1677714/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAGCATACCGGTCTGGACGCACACCTTGATTTATTACTTTCCACCCACACATTTGGTTATCCGAAAGAGGATCAGCGGTTATGGCATGCGGTGGCCGAAGCTACGGGTCTGAAAGCTGAAAGAACGCTGTTTATTGATGACAGCGAAGCGATTCTCGATGCTGCCGCGCAATTTGGTATTCGTTACTGCCTCGGCGTGACTAATCCTGATTCCGGGATTGCCGAGAAACAGTATCAACGCCATCCGTCACTGAATGACTACCGCCGCCTGATCCCCTCGCTAATGTGAAGGAGACGCCA  >  CP009273/3522409‑3522707

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: