Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,372,002 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,372,0020GC 5.3% 303.6 / 4.2 92intergenic (‑66/+31)carB/carAcarbamoyl‑phosphate synthase large subunit/carbamoyl‑phosphate synthase small subunit
Reads supporting (aligned to +/- strand):  ref base G (39/48);  new base C (3/2);  total (42/50)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTC  >  NC_002947/5371856‑5372149
                                                                                                                                                  |                                                                                                                                                     
gCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGCGCCCGGCAgcg                                                                                                                                                     >  3:338487/1‑147 (MQ=255)
gCAGGCCTGGCAGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGCGCCCGGCAgc                                                                                                                                                      <  5:305255/148‑2 (MQ=255)
 cAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGCGCCCGGCAgcgc                                                                                                                                                    <  4:338487/149‑4 (MQ=255)
    gCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                 >  2:5992/1‑149 (MQ=255)
    gCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                 >  8:266043/1‑149 (MQ=255)
    gCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                 >  2:86822/1‑149 (MQ=255)
    gCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGGGGTCACTCACATCCGGGCTCAGTCGGCAAGCCGCTGTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                 <  5:300760/149‑1 (MQ=255)
      tggcccGATCGTGATCGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGcc                                                                                                                                             >  3:204619/5‑149 (MQ=255)
          ccGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCCGGCACCTGCCgcgccctg                                                                                                                                          <  6:82473/149‑8 (MQ=255)
            gATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgggg                                                                                                                                       >  6:27358/1‑149 (MQ=255)
            gATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgggg                                                                                                                                       <  7:307334/149‑1 (MQ=255)
            gATCACGATCGGGCCAGCGCAGAGAATCAGGATGATTTTGATGTCTTTACGTTTTGGCATGGTGGTCAATCAAATCCGGGGTCAGTCGTCAAGACTATTTGAACAATGTGGGTCATGAGATGCAGGGAGAGGCAGGAGCAGGGGCgggg                                                                                                                                       >  6:256549/1‑149 (MQ=255)
              tCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGcc                                                                                                                                     <  4:190880/149‑1 (MQ=255)
                 cGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTg                                                                                                                                  <  6:230499/149‑1 (MQ=255)
                 cGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTg                                                                                                                                  <  1:386745/149‑1 (MQ=255)
                 cGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCAGCAGGTGCCGGCCCGGGGCCTTg                                                                                                                                  <  2:316290/149‑1 (MQ=255)
                  gATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTg                                                                                                                                  >  6:267589/1‑148 (MQ=255)
                  gATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGa                                                                                                                                 <  7:266043/149‑1 (MQ=255)
                   aTCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTg                                                                                                                                  >  8:399424/1‑147 (MQ=255)
                   aTCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                <  2:253838/149‑1 (MQ=255)
                   aTCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGaa                                                                                                                                <  8:181010/149‑1 (MQ=255)
                        gCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAg                                                                                                                            >  8:28675/1‑148 (MQ=255)
                               ccGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCt                                                                                                                      >  2:63342/1‑147 (MQ=255)
                                cGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAg                                                                                                                   <  8:313937/149‑1 (MQ=255)
                                  agaATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAgcg                                                                                                                 <  1:258782/149‑1 (MQ=255)
                                     aTCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGCGGGCCTTGAAGCAGGATGCTCAgcgg                                                                                                                <  1:5992/148‑1 (MQ=255)
                                      tCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGgcgc                                                                                                             >  8:172247/1‑149 (MQ=255)
                                      tCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCACGGGGCCTTGAAGCAGGATGCTCAGCGgcgc                                                                                                             >  3:317993/1‑149 (MQ=255)
                                      tCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGGCAGTCGGCAAGCCGCTTTGAACAATCTGGGGCAGGCGCTGCCGCGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGgcgc                                                                                                             <  6:329713/149‑1 (MQ=255)
                                         ggATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTg                                                                                                          <  8:343709/149‑1 (MQ=255)
                                         ggATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTg                                                                                                          >  5:318133/1‑149 (MQ=255)
                                         ggATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGACGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTg                                                                                                          <  1:380286/149‑1 (MQ=255)
                                         ggATGCTATTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTg                                                                                                          >  6:322483/1‑149 (MQ=255)
                                          gATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTgg                                                                                                         <  2:278091/149‑1 (MQ=255)
                                          gATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTgg                                                                                                         <  5:80979/149‑1 (MQ=255)
                                             gCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGTCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCa                                                                                                      >  1:336551/1‑149 (MQ=255)
                                              cTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCAt                                                                                                     <  5:256549/149‑1 (MQ=255)
                                               ttttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                    >  8:371171/1‑149 (MQ=255)
                                                tttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCa                                                                                                                    >  8:116611/1‑132 (MQ=255)
                                                tttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCa                                                                                                                    <  7:116611/132‑1 (MQ=255)
                                                 ttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                    <  7:172247/147‑1 (MQ=255)
                                                       tCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAg                                                                                                                               <  8:122255/114‑1 (MQ=255)
                                                       tCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAg                                                                                                                               >  7:122255/1‑114 (MQ=255)
                                                           tACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGaaa                                                                                        <  7:30875/149‑1 (MQ=255)
                                                             cGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACg                                                                                      <  7:399424/149‑1 (MQ=255)
                                                                 ttGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATc                                                                                   <  5:27358/148‑1 (MQ=255)
                                                                      aTGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACag                                                                             <  6:318133/149‑1 (MQ=255)
                                                                             tCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGcgcg                                                                      <  5:157375/149‑1 (MQ=255)
                                                                              cACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGa                                                                     >  8:309836/1‑149 (MQ=255)
                                                                               aCTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAgcggc                                                                                                               >  6:308339/1‑106 (MQ=255)
                                                                               aCTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAgcggc                                                                                                               <  5:308339/106‑1 (MQ=255)
                                                                                 tCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgt                                                                  >  3:336647/1‑149 (MQ=255)
                                                                                   aaaTCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcg                                                                >  6:427608/1‑149 (MQ=255)
                                                                                   aaaTCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcg                                                                >  1:297617/1‑149 (MQ=255)
                                                                                      tCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcggtc                                                             <  1:290270/149‑1 (MQ=255)
                                                                                        cGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGtcggtcgg                                                           <  4:204619/149‑1 (MQ=255)
                                                                                          gggTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACagag                                                                           <  2:150724/131‑1 (MQ=255)
                                                                                          gggTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACagag                                                                           >  1:150724/1‑131 (MQ=255)
                                                                                            gTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGccc                                                       >  2:126171/1‑149 (MQ=255)
                                                                                               aGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCggg                                                    >  7:84242/1‑149 (MQ=255)
                                                                                                 tCGGCAAGTCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCt                                                  <  8:345419/149‑1 (MQ=255)
                                                                                                 tCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                    >  4:146362/1‑99 (MQ=255)
                                                                                                 tCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                    <  3:146362/99‑1 (MQ=255)
                                                                                                 tCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCt                                                  >  4:271346/1‑149 (MQ=255)
                                                                                                   ggCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCg                                                <  2:178102/149‑1 (MQ=255)
                                                                                                     cAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCt                                              >  6:272185/1‑149 (MQ=255)
                                                                                                       aGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGTCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTc                                            <  1:126171/149‑1 (MQ=255)
                                                                                                        gCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGa                                                                                                                          >  5:295559/1‑70 (MQ=255)
                                                                                                        gCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGa                                                                                                                          <  6:295559/70‑1 (MQ=255)
                                                                                                              ttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTgg                                     >  8:182370/1‑149 (MQ=255)
                                                                                                              ttGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTgg                                     >  4:203714/1‑149 (MQ=255)
                                                                                                               tGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTgg                                     <  1:63342/148‑1 (MQ=255)
                                                                                                                 aaCAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGccc                                  <  7:309836/149‑1 (MQ=255)
                                                                                                                  aCAATCTGGTTCAGGGGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCAATGGCATCGTTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGccc                                  <  2:368489/148‑1 (MQ=255)
                                                                                                                  aCAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGccc                                  <  1:86822/148‑1 (MQ=255)
                                                                                                                   cAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTg                                >  1:277074/1‑149 (MQ=255)
                                                                                                                   cAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTg                                <  1:299922/149‑1 (MQ=255)
                                                                                                                      tCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTgg                               >  1:111586/1‑147 (MQ=255)
                                                                                                                             cAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAc                      <  6:379208/149‑1 (MQ=255)
                                                                                                                              aGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAcg                     <  2:277074/149‑1 (MQ=255)
                                                                                                                              aGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAAc                      <  4:93113/148‑1 (MQ=255)
                                                                                                                                gcgcTGCCGGGGGCGGCAGGTGCCGGCCTGTGGCCTTGAAGCAGGATTCTCAGGGGCGGTTGGCCATGGCATGTGTGTTACGATCTTACAGAGGCGCGTCGTGGGGCGGGCGCGGGCTCGCTTCAGGGTGGCCCTGGAATCTGAAcgag                   <  4:317993/149‑3 (MQ=255)
                                                                                                                                 cgcTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACgcgc                  <  7:28675/149‑1 (MQ=255)
                                                                                                                                    ccccGGGCCGGCACCTGCCGCGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGctct               >  7:372996/3‑149 (MQ=255)
                                                                                                                                      ccGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTg             >  6:88886/1‑149 (MQ=255)
                                                                                                                                       cGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGt            >  1:379940/1‑149 (MQ=255)
                                                                                                                                        cgggccgGCACCTGCCGCGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGt            >  1:381211/6‑149 (MQ=255)
                                                                                                                                         ggCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCg          <  2:336551/149‑1 (MQ=255)
                                                                                                                                          ggccgGCACCTGCCGCGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTgg                                     <  7:84621/119‑1 (MQ=255)
                                                                                                                                          ggccgGCACCTGCCGCGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTgg                                     >  8:84621/4‑122 (MQ=255)
                                                                                                                                            gcgGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAgg                                         <  8:201252/115‑1 (MQ=255)
                                                                                                                                            gcgGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAgg                                         >  7:201252/1‑115 (MQ=255)
                                                                                                                                              ggCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACg     >  4:364448/1‑149 (MQ=255)
                                                                                                                                               gCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGc    >  1:408664/1‑149 (MQ=255)
                                                                                                                                                 aGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTc  <  4:336647/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGGCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTC  >  NC_002947/5371856‑5372149

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