Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,373,972 G→C 5.3% intergenic (‑66/+31) 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,9720GC 5.3% 630.2 / 18.3 188intergenic (‑66/+31)carB/carAcarbamoyl‑phosphate synthase large subunit/carbamoyl‑phosphate synthase small subunit
Reads supporting (aligned to +/- strand):  ref base G (100/78);  new base C (5/5);  total (105/83)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 197 total aligned reads displayed.

TTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGA  >  NC_002947/5373823‑5374121
                                                                                                                                                     |                                                                                                                                                      
ttCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAg                                                                                                                                                        >  1:566579/1‑150 (MQ=255)
 tCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAgg                                                                                                                                                       <  2:3512980/150‑1 (MQ=255)
   gCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTg                                                                                                                                                     >  1:1483385/1‑150 (MQ=255)
     aGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGcc                                                                                                                                                   >  1:576040/1‑150 (MQ=255)
         cTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGcc                                                                                                                                              >  1:516882/1‑150 (MQ=255)
          tGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACACTCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGccc                                                                                                                                             >  1:5080073/1‑150 (MQ=255)
          tGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCgcgccc                                                                                                                                             <  2:2815083/150‑6 (MQ=255)
            gCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgg                                                                                                                                           <  1:2328760/150‑1 (MQ=255)
             ccGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCggg                                                                                                                                          <  2:1940979/150‑1 (MQ=255)
               gATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGc                                                                                                                                        <  2:5066107/150‑1 (MQ=255)
                aTCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGcc                                                                                                                                       <  2:2510826/150‑1 (MQ=255)
                 tCACGCCCGGGCCAGCGCCGAGAATCAGGATGCTTTGGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGTGCCt                                                                                                                                      <  2:4817837/150‑1 (MQ=255)
                 tCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCt                                                                                                                                      >  1:4745996/1‑150 (MQ=255)
                  cACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCtt                                                                                                                                     <  2:1225485/150‑1 (MQ=255)
                    gtgatcggCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGc                                                                                                                                                    >  1:4363898/7‑134 (MQ=255)
                    gtgatcggCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGc                                                                                                                                                    <  2:4363862/128‑1 (MQ=255)
                    gtgatcggCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGa                                                                                                                                   >  1:4259776/7‑150 (MQ=255)
                     gATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                  >  1:2669543/1‑150 (MQ=255)
                     gATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                  >  2:3434816/1‑150 (MQ=255)
                       tCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCgcgcc                                                                                                                                              <  1:3371603/136‑5 (MQ=255)
                       tCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCgcgcc                                                                                                                                              >  2:3371603/1‑132 (MQ=255)
                        cGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCTGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCa                                                                                                                               <  1:4657528/150‑1 (MQ=255)
                           gCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGa                                                                                                                            >  1:1017519/1‑150 (MQ=255)
                             cAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                  >  1:3329522/1‑142 (MQ=255)
                                 gcCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCa                                                                                                                      >  2:3999291/1‑150 (MQ=255)
                                 gcCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCa                                                                                                                      >  2:4183323/1‑150 (MQ=255)
                                                gCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCAt                                                                                                       <  2:3775281/150‑1 (MQ=255)
                                                 cTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                      >  1:306169/1‑150 (MQ=255)
                                                  ttttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATgg                                                                                                     >  1:1461457/1‑150 (MQ=255)
                                                  ttttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATgg                                                                                                     >  2:1860221/1‑150 (MQ=255)
                                                  ttttGATGTCTGGACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGc                                                                                                                                <  2:4139489/123‑1 (MQ=255)
                                                  ttttGATGTCTGAACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGc                                                                                                                                >  1:4139489/1‑123 (MQ=255)
                                                         gTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGt                                                                                              >  2:2167387/1‑150 (MQ=255)
                                                           cTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGa                                                                                            <  2:5176303/150‑1 (MQ=255)
                                                            tGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGTGAAGGCGATGCCGGGCGCGCCAGTTGCCGCCCCGGGCCCTGGAGGCGGGATCATCGCCGGCGCTTGGCAATGGCAACGGTGaa                                                                                           >  1:2841518/1‑150 (MQ=255)
                                                             gTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGaaa                                                                                          >  1:592825/1‑150 (MQ=255)
                                                             gTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCGCAGCGGCGCTTGGCCATGGCATCGGTGaac                                                                                          >  2:4379317/1‑149 (MQ=255)
                                                               aCGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACg                                                                                        >  1:2935352/1‑150 (MQ=255)
                                                                 gTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAg                                                                                                                                                        <  2:4724855/85‑1 (MQ=255)
                                                                 gTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGAt                                                                                      >  1:2032760/1‑150 (MQ=255)
                                                                  ttttGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATc                                                                                     >  1:693320/1‑150 (MQ=255)
                                                                      ggCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAAc                                                                                 >  1:2906532/1‑150 (MQ=255)
                                                                       gCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCAAGGCCTCGGTGAAACGATCGAACa                                                                                >  2:2738797/1‑150 (MQ=255)
                                                                         aTGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACaga                                                                              >  1:3611002/1‑150 (MQ=255)
                                                                          tggtggTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACagag                                                                             <  1:3895973/150‑1 (MQ=255)
                                                                              ggtCACACAAATCCGGGGTCAGTCGGCAAGCCGCTATGAACAATCTGCGACAGGCGCTGCCGGGAGCGGAAGGTGACGGCCCGGGGCCTTGAAGCAGGATGCACAGCGGCGCTTGGCCATGGCATAGGTGAAACGATCGAACAGAGgcgc                                                                         >  1:4656104/1‑150 (MQ=255)
                                                                                tCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACg                                                                                        >  1:4079152/1‑133 (MQ=255)
                                                                                tCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGa                                                                       >  1:3576985/1‑150 (MQ=255)
                                                                                  aCTCAAATCCGGGGTCAGAAGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACg                                                                     >  2:1659361/1‑150 (MQ=255)
                                                                                    tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCGTTGAAGCTGGATGCTCAGCGGCGCTTGGCCTTGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                   >  2:2667249/1‑150 (MQ=255)
                                                                                    tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATc                                                                                     >  1:4223195/1‑132 (MQ=255)
                                                                                    tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                   >  2:3032255/1‑150 (MQ=255)
                                                                                    tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                   >  2:4067922/1‑150 (MQ=255)
                                                                                    tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                   >  1:302694/1‑150 (MQ=255)
                                                                                    tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtc                                                                   >  1:4435043/1‑150 (MQ=255)
                                                                                     cAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcg                                                                  <  2:4623274/150‑1 (MQ=255)
                                                                                      aaaTCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcgg                                                                 >  1:3973475/1‑150 (MQ=255)
                                                                                        aTCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcggtc                                                               <  2:2453936/150‑1 (MQ=255)
                                                                                         tCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcggtcg                                                              >  2:956989/1‑150 (MQ=255)
                                                                                          ccGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGtcggtcgg                                                             <  2:402007/150‑1 (MQ=255)
                                                                                                  aGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGc                                                     >  1:225133/1‑150 (MQ=255)
                                                                                                   gTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGaa                                                                                           >  2:134908/1‑111 (MQ=255)
                                                                                                       gCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCt                                                >  2:145647/1‑150 (MQ=255)
                                                                                                        cAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCtt                                               <  1:2290157/150‑1 (MQ=255)
                                                                                                           gCCGCTTTGAACATTCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAg                                            <  2:2450339/150‑1 (MQ=255)
                                                                                                               cTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTg                                        <  1:3211333/150‑1 (MQ=255)
                                                                                                               cTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTg                                        <  2:1538910/150‑1 (MQ=255)
                                                                                                               cTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTg                                        <  1:3511584/150‑1 (MQ=255)
                                                                                                                tttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTgg                                       <  1:3999291/150‑1 (MQ=255)
                                                                                                                 ttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGc                                      >  2:3487128/1‑150 (MQ=255)
                                                                                                                 ttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGc                                      >  1:1316732/1‑150 (MQ=255)
                                                                                                                 ttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGc                                      <  1:4782781/150‑1 (MQ=255)
                                                                                                                  tGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGcc                                     <  2:2986010/150‑1 (MQ=255)
                                                                                                                   gAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCAAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGccc                                    >  1:4263406/1‑150 (MQ=255)
                                                                                                                      cAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTgg                                 <  1:1824615/150‑1 (MQ=255)
                                                                                                                      cAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTgg                                 <  2:4259776/150‑1 (MQ=255)
                                                                                                                        aTCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGaa                               <  2:2284694/150‑1 (MQ=255)
                                                                                                                         tCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAg                              <  2:4656176/150‑1 (MQ=255)
                                                                                                                           tGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCt                            <  2:850514/150‑1 (MQ=255)
                                                                                                                           tGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCAAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCt                            <  2:4263406/150‑1 (MQ=255)
                                                                                                                            gggTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTg                           >  2:1096613/1‑150 (MQ=255)
                                                                                                                            gggTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTg                           <  2:3865852/150‑1 (MQ=255)
                                                                                                                              gTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGGGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGaa                         <  2:5418173/150‑1 (MQ=255)
                                                                                                                               tCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAc                        <  1:3144854/150‑1 (MQ=255)
                                                                                                                               tCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAc                        >  1:2308877/1‑150 (MQ=255)
                                                                                                                               tCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAc                        <  1:2230981/150‑1 (MQ=255)
                                                                                                                               tCAGGCGCTGCCGGGCGCGACAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTc                                              >  2:3736475/1‑128 (MQ=255)
                                                                                                                                  ggCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAcgcg                     <  1:5373940/150‑1 (MQ=255)
                                                                                                                                  ggCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAcgcg                     >  1:820168/1‑150 (MQ=255)
                                                                                                                                   gcgcTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACgcgc                    >  2:2914641/1‑150 (MQ=255)
                                                                                                                                   gcgcTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACgcgc                    >  2:3681477/1‑150 (MQ=255)
                                                                                                                                         ccGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGt              <  2:566579/150‑1 (MQ=255)
                                                                                                                                          cGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTc             >  1:2581377/1‑150 (MQ=255)
                                                                                                                                            gggccgGCACCTGCCGCGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCg            >  2:314112/5‑150 (MQ=255)
                                                                                                                                              cgcgGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTa         <  1:3434816/150‑1 (MQ=255)
                                                                                                                                                 ggCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGc      >  2:1421478/1‑150 (MQ=255)
                                                                                                                                                  gCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCt     >  2:2017100/1‑150 (MQ=255)
                                                                                                                                                    aGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCg   >  1:2474812/1‑150 (MQ=255)
                                                                                                                                                     ggTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGa  <  1:2973657/150‑1 (MQ=255)
                                                                                                                                                     ggTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGa  <  2:3906580/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TTCGCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGA  >  NC_002947/5373823‑5374121

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: