Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,578,959 C→A 6.8% intergenic (‑75/‑131) aegA ← / → narQ putative oxidoreductase, FeS binding subunit/NAD/FAD‑binding subunit/sensory histidine kinase in two‑component regulatory system with NarP (NarL)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,578,9590CA 6.8% 1058.6 / 30.9 283intergenic (‑75/‑131)aegA/narQputative oxidoreductase, FeS binding subunit/NAD/FAD‑binding subunit/sensory histidine kinase in two‑component regulatory system with NarP (NarL)
Reads supporting (aligned to +/- strand):  ref base C (135/128);  new base A (14/6);  total (149/134)
Fisher's exact test for biased strand distribution p-value = 1.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

GCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACC  >  CP009273/2578811‑2579106
                                                                                                                                                    |                                                                                                                                                   
gCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATc                                                                                                                                                     <  2:273486/149‑1 (MQ=255)
  cATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCgg                                                                                                                                                   >  2:507256/1‑149 (MQ=255)
        acaGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGt                                                                                                                                             <  2:2720953/149‑1 (MQ=255)
            gCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCAt                                                                                                                                         <  2:2460784/149‑1 (MQ=255)
            gCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCAt                                                                                                                                         <  2:1052152/149‑1 (MQ=255)
             cGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCAtg                                                                                                                                        <  1:3243087/149‑1 (MQ=255)
             cGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCAtg                                                                                                                                        <  1:2948127/149‑1 (MQ=255)
             cGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCAtg                                                                                                                                        <  1:2039899/149‑1 (MQ=255)
               aTTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCAtgtg                                                                                                                                      >  2:3928257/1‑149 (MQ=255)
                  tCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGccc                                                                                                                                   <  2:1395618/149‑1 (MQ=255)
                    acaAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAg                                                                                                                                 <  2:868608/149‑1 (MQ=255)
                     caAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGa                                                                                                                                >  1:1538360/1‑149 (MQ=255)
                        gCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATa                                                                                                                             <  2:2968577/149‑1 (MQ=255)
                        gCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGGTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATa                                                                                                                             >  2:3001177/1‑149 (MQ=255)
                         cATGACAACCCAGACATTGCTGACTGTGGGCCATAATAACACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATaa                                                                                                                            <  1:2475388/149‑1 (MQ=255)
                          aTGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAg                                                                                                                           <  1:2545030/149‑1 (MQ=255)
                           tGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGt                                                                                                                          <  2:2708330/149‑1 (MQ=255)
                           tGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGt                                                                                                                          <  2:252769/149‑1 (MQ=255)
                             aCAACCCAGACATTGCTGACTGTTGGCCATAATAACACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTaa                                                                                                                        <  2:3823558/149‑1 (MQ=255)
                             aCAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTaa                                                                                                                        <  1:79370/149‑1 (MQ=255)
                             aCAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTaa                                                                                                                        >  1:4282491/1‑149 (MQ=255)
                             aCAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTaa                                                                                                                        >  2:4348310/1‑149 (MQ=255)
                               aaCCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACt                                                                                                                      <  2:1535906/149‑1 (MQ=255)
                                aCCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTa                                                                                                                     <  1:2184936/149‑1 (MQ=255)
                                aCCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTa                                                                                                                     <  1:622864/149‑1 (MQ=255)
                                 cccAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTAt                                                                                                                    >  2:4690742/1‑149 (MQ=255)
                                   cAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTAttt                                                                                                                  <  2:370482/149‑1 (MQ=255)
                                    aGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTc                                                                                                                 >  2:3393255/1‑149 (MQ=255)
                                    aGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTc                                                                                                                 >  2:2049547/1‑149 (MQ=255)
                                     gACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCg                                                                                                                <  2:4334680/149‑1 (MQ=255)
                                     gACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCg                                                                                                                >  2:4744090/1‑149 (MQ=255)
                                     gACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCg                                                                                                                <  2:4534986/149‑1 (MQ=255)
                                        aTTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCta                                                                                                             <  2:3517561/149‑1 (MQ=255)
                                         ttGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCtat                                                                                                            >  2:1339089/1‑149 (MQ=255)
                                         ttGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTAACCTACTTAAAGGGGAATAGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTGGCtat                                                                                                            >  2:4156279/1‑149 (MQ=255)
                                             tGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAAc                                                                                                        <  1:3916319/149‑1 (MQ=255)
                                             tGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAAc                                                                                                        <  1:2708679/149‑1 (MQ=255)
                                                  gTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTgg                                                                                                   <  1:3147813/149‑1 (MQ=255)
                                                   ttGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGCATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATGTCGCTATAAACTGTGGc                                                                                                  >  1:4093735/1‑149 (MQ=255)
                                                   ttGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGc                                                                                                  <  2:3565464/149‑1 (MQ=255)
                                                        cATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATc                                                                                             <  2:445947/149‑1 (MQ=255)
                                                             taaAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaata                                                                                        <  1:4680592/149‑1 (MQ=255)
                                                               aaaCGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAaataat                                                                                      >  1:3009490/1‑149 (MQ=255)
                                                               aaaCGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACACGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTTGCAGATCAaataat                                                                                      >  1:1823895/1‑149 (MQ=255)
                                                                aaCGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATc                                                                                     >  1:161621/1‑149 (MQ=255)
                                                                aaCGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATc                                                                                     >  1:2059634/1‑149 (MQ=255)
                                                                aaCGATTAATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATc                                                                                     <  2:2733249/149‑1 (MQ=255)
                                                                 aCGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATcc                                                                                    >  1:1250238/1‑149 (MQ=255)
                                                                    aTTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCAt                                                                                 <  2:1080950/149‑1 (MQ=255)
                                                                     ttCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATc                                                                                <  1:8295/149‑1 (MQ=255)
                                                                      tCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCa                                                                               <  1:4538150/149‑1 (MQ=255)
                                                                       cATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATAGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCAc                                                                              <  2:631187/149‑1 (MQ=255)
                                                                         tAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTg                                                                            >  1:550848/1‑149 (MQ=255)
                                                                           aCGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGAc                                                                          <  1:346118/149‑1 (MQ=255)
                                                                           aCGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGAc                                                                          >  2:4050193/1‑149 (MQ=255)
                                                                            cGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACt                                                                         <  1:4306790/149‑1 (MQ=255)
                                                                            cGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACt                                                                         <  2:3889263/149‑1 (MQ=255)
                                                                              aCCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTaa                                                                       <  2:2323109/149‑1 (MQ=255)
                                                                                cttcttTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAAt                                                                     >  2:2588291/1‑149 (MQ=255)
                                                                                 ttcttTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAAtt                                                                    >  2:1838452/1‑149 (MQ=255)
                                                                                         tGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATAGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCa                                                            <  1:713727/149‑1 (MQ=255)
                                                                                           gTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGc                                                          <  2:1772986/149‑1 (MQ=255)
                                                                                           gTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGc                                                          >  2:749543/1‑149 (MQ=255)
                                                                                            ttGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATAGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCg                                                         <  1:1313293/149‑1 (MQ=255)
                                                                                               tAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAAc                                                      <  1:816026/149‑1 (MQ=255)
                                                                                                 aaaaCCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTg                                                    >  2:1736884/1‑149 (MQ=255)
                                                                                                       ttATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACAt                                              >  2:3106833/1‑149 (MQ=255)
                                                                                                       ttATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACAt                                              >  1:2981272/1‑149 (MQ=255)
                                                                                                          ttCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTaa                                           >  2:2212820/1‑149 (MQ=255)
                                                                                                           tCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAAt                                          <  1:3693664/149‑1 (MQ=255)
                                                                                                           tCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAAt                                          <  2:1954082/149‑1 (MQ=255)
                                                                                                             tttATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGa                                        >  2:3611824/1‑149 (MQ=255)
                                                                                                               tataTGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGAtt                                      >  1:2519244/1‑149 (MQ=255)
                                                                                                                 tatGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGAtttt                                    >  1:1503091/1‑149 (MQ=255)
                                                                                                                     aGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTtgt                                <  2:4402426/149‑1 (MQ=255)
                                                                                                                       tgttgtTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTgg                              >  1:113592/1‑149 (MQ=255)
                                                                                                                       tgttgtTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTgg                              >  1:1068963/1‑149 (MQ=255)
                                                                                                                           gttATTAGGCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGagaa                          <  2:2762271/149‑1 (MQ=255)
                                                                                                                            ttattaCCCGACTTACAGGGGAATAGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGagaag                         >  1:65956/1‑149 (MQ=255)
                                                                                                                             tattaCCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGagaaga                        <  2:2627577/149‑1 (MQ=255)
                                                                                                                             tattaCCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGagaaga                        >  2:3749319/1‑149 (MQ=255)
                                                                                                                              attaCCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGAc                       >  2:387511/1‑149 (MQ=255)
                                                                                                                              attaCCCGACTTACAGGGGAATAGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGAc                       >  2:1576724/1‑149 (MQ=255)
                                                                                                                                  cccGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCgt                   >  2:1167572/1‑149 (MQ=255)
                                                                                                                                  cccGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAAGTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCgt                   <  2:3204642/149‑1 (MQ=255)
                                                                                                                                     gACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGtgtg                >  2:122172/1‑149 (MQ=255)
                                                                                                                                     gACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGtgtg                >  2:819284/1‑149 (MQ=255)
                                                                                                                                      aCTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGa               >  1:1673856/1‑149 (MQ=255)
                                                                                                                                       cTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGAt              <  1:4246396/149‑1 (MQ=255)
                                                                                                                                       cTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGAt              >  2:2855583/1‑149 (MQ=255)
                                                                                                                                        ttACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGAtt             >  1:3984493/1‑149 (MQ=255)
                                                                                                                                        ttACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGAtt             <  2:2059669/149‑1 (MQ=255)
                                                                                                                                               ggAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAAc      >  2:1692196/1‑149 (MQ=255)
                                                                                                                                               ggAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAAc      >  1:2711175/1‑149 (MQ=255)
                                                                                                                                                gAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACg     >  1:3108385/1‑149 (MQ=255)
                                                                                                                                                gAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACg     <  2:3321591/149‑1 (MQ=255)
                                                                                                                                                  aTCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGAc   <  2:1013742/149‑1 (MQ=255)
                                                                                                                                                   tCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGAcc  <  2:3418801/149‑1 (MQ=255)
                                                                                                                                                   tCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTGGTGGAGAAGACGCGTGTGATTGTTAAACGAcc  <  2:547872/149‑1 (MQ=255)
                                                                                                                                                   tCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATAAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGATGACGCGTGTGATTGTTAAACGAcc  <  2:4210610/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCCATGACACAGGCGATTTCACAAGCATGACAACCCAGACATTGCTGACTGTTGGCCATAATAAAACGATTCATAACGACCTTCTTTTTTGGTTGTAAAAACCTTATTCTTTATATGAGTGTTGTTATTACCCGACTTACAGGGGAATCGGCAATGTTCATGTGCCCAGAATAAGTAACTATTTCGCTATAAACTGTGGCAGATCAAATAATCCCATCACTGACTAAATTGCGTTTCAGCGAACTGGAACATTAATGATTTTTTGTGGAGAAGACGCGTGTGATTGTTAAACGACC  >  CP009273/2578811‑2579106

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: