Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,373,973 G→C 5.3% intergenic (‑67/+30) carB ← / ← carA carbamoyl‑phosphate synthase large subunit/carbamoyl‑phosphate synthase small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,373,9730GC 5.3% 627.3 / 18.3 188intergenic (‑67/+30)carB/carAcarbamoyl‑phosphate synthase large subunit/carbamoyl‑phosphate synthase small subunit
Reads supporting (aligned to +/- strand):  ref base G (99/78);  new base C (5/5);  total (105/83)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 194 total aligned reads displayed.

TCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGAT  >  NC_002947/5373824‑5374122
                                                                                                                                                     |                                                                                                                                                      
  gCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTg                                                                                                                                                      >  1:1483385/1‑150 (MQ=255)
   cAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGc                                                                                                                                                     <  2:3058634/150‑1 (MQ=255)
    aGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGcc                                                                                                                                                    >  1:576040/1‑150 (MQ=255)
     ggCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                   <  1:4007396/150‑1 (MQ=255)
        tggcccGATCGTGATCGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGccc                                                                                                                                              >  2:1668589/5‑150 (MQ=255)
         tGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGccc                                                                                                                                              <  1:3034858/150‑1 (MQ=255)
         tGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCgcgccc                                                                                                                                              <  2:2815083/150‑6 (MQ=255)
           gCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgg                                                                                                                                            <  1:399567/150‑1 (MQ=255)
           gCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgg                                                                                                                                            <  1:2328760/150‑1 (MQ=255)
           gCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGGTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCACGCCGCGTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgg                                                                                                                                            >  2:2200935/1‑150 (MQ=255)
             cGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgggg                                                                                                                                          <  2:1922740/150‑1 (MQ=255)
              gATCACGATCGGTCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGc                                                                                                                                         <  2:249665/150‑1 (MQ=255)
              gATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGc                                                                                                                                         <  2:5066107/150‑1 (MQ=255)
               aTCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCGCGCCCggcagc                                                                                                                                        <  2:2381181/150‑5 (MQ=255)
                tCACGCCCGGGCCAGCGCCGAGAATCAGGATGCTTTGGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGTGCCt                                                                                                                                       <  2:4817837/150‑1 (MQ=255)
                 cACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCtt                                                                                                                                      <  2:1225485/150‑1 (MQ=255)
                   gtgatcggCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGc                                                                                                                                                     <  2:4363862/128‑1 (MQ=255)
                   gtgatcggCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGc                                                                                                                                                     >  1:4363862/7‑134 (MQ=255)
                   gtgatcggCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGa                                                                                                                                    >  1:4259811/7‑150 (MQ=255)
                    gATCGGGCGAGCGCCGAGTATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCGTTGAACAATCTGGGTCAGGCGCTACCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                   <  2:3890083/150‑1 (MQ=255)
                    gATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                   <  2:3830674/150‑1 (MQ=255)
                    gATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                   >  1:2669543/1‑150 (MQ=255)
                       cGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCTGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCa                                                                                                                                <  1:4657528/150‑1 (MQ=255)
                          gCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGa                                                                                                                             >  1:1017519/1‑150 (MQ=255)
                          gCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGa                                                                                                                             >  1:415647/1‑150 (MQ=255)
                            cAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                   <  2:3329522/142‑1 (MQ=255)
                                gcCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCa                                                                                                                       >  2:4183323/1‑150 (MQ=255)
                                        tCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCt                                                                                                               >  1:3906580/1‑150 (MQ=255)
                                           ggATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTgg                                                                                                            >  2:288830/1‑150 (MQ=255)
                                            gATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGc                                                                                                           >  1:4380207/1‑150 (MQ=255)
                                              tGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCa                                                                                                         <  1:4132707/150‑1 (MQ=255)
                                                 ttttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATgg                                                                                                      >  2:1227299/1‑150 (MQ=255)
                                                 ttttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATgg                                                                                                      >  1:1461457/1‑150 (MQ=255)
                                                 ttttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATgg                                                                                                      >  1:3070578/1‑150 (MQ=255)
                                                 ttttGATGTCTGGACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGc                                                                                                                                 <  2:4139489/123‑1 (MQ=255)
                                                    tGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCgcgccc                                                                                                                                              <  1:3016969/107‑6 (MQ=255)
                                                    tGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCgcgccc                                                                                                                                              >  2:3016969/1‑102 (MQ=255)
                                                          cTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGa                                                                                             <  2:5176303/150‑1 (MQ=255)
                                                           tGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGTGAAGGCGATGCCGGGCGCGCCAGTTGCCGCCCCGGGCCCTGGAGGCGGGATCATCGCCGGCGCTTGGCAATGGCAACGGTGaa                                                                                            >  1:2841518/1‑150 (MQ=255)
                                                            gTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGaac                                                                                           >  2:3706926/1‑149 (MQ=255)
                                                            gTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGaaa                                                                                           <  2:2562126/150‑1 (MQ=255)
                                                            gTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGaaa                                                                                           >  1:592825/1‑150 (MQ=255)
                                                            gTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCGCAGCGGCGCTTGGCCATGGCATCGGTGaac                                                                                           >  2:4379317/1‑149 (MQ=255)
                                                              aCGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACg                                                                                         >  1:2935352/1‑150 (MQ=255)
                                                                gTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGAt                                                                                       >  1:2032760/1‑150 (MQ=255)
                                                                gTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGAt                                                                                       >  1:5334224/1‑150 (MQ=255)
                                                                 ttttGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATc                                                                                      >  1:693320/1‑150 (MQ=255)
                                                                     ggCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAAc                                                                                  >  1:2906532/1‑150 (MQ=255)
                                                                      gCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACa                                                                                 >  2:2738778/1‑150 (MQ=255)
                                                                          ggtggtCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAgg                                                                             >  1:3251524/1‑150 (MQ=255)
                                                                             ggtCACACAAATCCGGGGTCAGTCGGCAAGCCGCTATGAACAATCTGCGACAGGCGCTGCCGGGAGCGGAAGGTGACGGCCCGGGGCCTTGAAGCAGGATGCACAGCGGCGCTTGGCCATGGCATAGGTGAAACGATCGAACAGAGgcgc                                                                          >  1:4656104/1‑150 (MQ=255)
                                                                              gtCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGcgcg                                                                         >  2:3451139/1‑150 (MQ=255)
                                                                               tCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACg                                                                                         <  2:4079152/133‑1 (MQ=255)
                                                                                 aCTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACg                                                                      >  2:2004383/1‑150 (MQ=255)
                                                                                 aCTCAAATCCGGGGTCAGTCGGCAAGCCGCTGTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCCGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGATCCGATCGAACAGAGGCGCGACg                                                                      >  2:3638116/1‑150 (MQ=255)
                                                                                   tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATc                                                                                      >  1:4223195/1‑132 (MQ=255)
                                                                                   tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                    >  1:5138882/1‑150 (MQ=255)
                                                                                   tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                    >  1:3499664/1‑150 (MQ=255)
                                                                                     aaaTCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcgg                                                                  >  1:2962893/1‑150 (MQ=255)
                                                                                        tCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcggtcg                                                               >  2:956989/1‑150 (MQ=255)
                                                                                         ccGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGtcggtcgg                                                              <  2:402324/150‑1 (MQ=255)
                                                                                          cGGGGACAGTCGTCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCTCGACGTCGGTCggg                                                             <  2:4002674/150‑1 (MQ=255)
                                                                                             ggTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGccc                                                          <  2:1480565/150‑1 (MQ=255)
                                                                                             ggTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGccc                                                          <  2:1402684/150‑1 (MQ=255)
                                                                                                 aGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGc                                                      >  1:225133/1‑150 (MQ=255)
                                                                                                      gCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCt                                                 >  2:145292/1‑150 (MQ=255)
                                                                                                       cAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCtt                                                <  1:2290157/150‑1 (MQ=255)
                                                                                                          gCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                    >  1:3389577/1‑127 (MQ=255)
                                                                                                          gCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                    <  2:3389577/127‑1 (MQ=255)
                                                                                                          gCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAg                                             <  2:2450450/150‑1 (MQ=255)
                                                                                                              cTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTg                                         <  1:3512116/150‑1 (MQ=255)
                                                                                                                ttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGc                                       >  1:1316732/1‑150 (MQ=255)
                                                                                                                ttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGc                                       <  1:4782742/150‑1 (MQ=255)
                                                                                                                  gAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGccc                                     <  2:3499664/150‑1 (MQ=255)
                                                                                                                   aaCAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCt                                    <  2:1636986/150‑1 (MQ=255)
                                                                                                                     cAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTgg                                  <  1:1824615/150‑1 (MQ=255)
                                                                                                                     cAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTgg                                  <  2:4259776/150‑1 (MQ=255)
                                                                                                                        tCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAg                               <  2:4656104/150‑1 (MQ=255)
                                                                                                                          tGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCt                             <  2:850514/150‑1 (MQ=255)
                                                                                                                          tGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCt                             <  1:31961/150‑1 (MQ=255)
                                                                                                                          tGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCAAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCt                             <  2:4263406/150‑1 (MQ=255)
                                                                                                                           gggTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTg                            >  2:1096613/1‑150 (MQ=255)
                                                                                                                             gTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGGGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGaa                          <  2:5418173/150‑1 (MQ=255)
                                                                                                                             gTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGaa                          >  1:3329796/1‑150 (MQ=255)
                                                                                                                              tCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAc                         <  1:2230981/150‑1 (MQ=255)
                                                                                                                              tCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAc                         >  1:2308877/1‑150 (MQ=255)
                                                                                                                              tCAGGCGCTGCCGGGCGCGACAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTc                                               <  1:3736475/128‑1 (MQ=255)
                                                                                                                                 ggCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAcgcg                      >  1:820168/1‑150 (MQ=255)
                                                                                                                                   cgcTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGct                    >  2:2973657/1‑150 (MQ=255)
                                                                                                                                   cgcTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGct                    <  2:3973475/150‑1 (MQ=255)
                                                                                                                                    gcTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGctc                   <  1:4817192/150‑1 (MQ=255)
                                                                                                                                    gcTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGctc                   >  2:3747238/1‑150 (MQ=255)
                                                                                                                                       cccGGGCCGGCACCTGCCGCGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCAt                                                                                                        <  2:2115190/61‑1 (MQ=21)
                                                                                                                                        ccGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGt               <  2:566579/150‑1 (MQ=255)
                                                                                                                                         cGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTc              >  1:2581377/1‑150 (MQ=255)
                                                                                                                                                   aGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCg    >  1:2474812/1‑150 (MQ=255)
                                                                                                                                                   aGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCg    >  2:4699501/1‑150 (MQ=255)
                                                                                                                                                    ggTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGa   >  1:54395/1‑150 (MQ=255)
                                                                                                                                                    ggTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGa   >  1:3602666/1‑150 (MQ=255)
                                                                                                                                                     gTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCGTGTCGGTACGCTCGAt  >  2:4436199/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGAT  >  NC_002947/5373824‑5374122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: