Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,372,003 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,372,0030GC 5.3% 305.9 / 4.2 92intergenic (‑67/+30)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 = 8.45e-01

GCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGA  >  NC_002947/5371860‑5372151
                                                                                                                                               |                                                                                                                                                      
gCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                   >  2:86822/1‑149 (MQ=255)
gCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                   >  2:5992/1‑149 (MQ=255)
gCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCg                                                                                                                                                   >  8:266043/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                                                                                                                                         <  7:307334/149‑1 (MQ=255)
        gATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCgggg                                                                                                                                         >  6:27358/1‑149 (MQ=255)
        gATCACGATCGGGCCAGCGCAGAGAATCAGGATGATTTTGATGTCTTTACGTTTTGGCATGGTGGTCAATCAAATCCGGGGTCAGTCGTCAAGACTATTTGAACAATGTGGGTCATGAGATGCAGGGAGAGGCAGGAGCAGGGGCgggg                                                                                                                                         >  6:256549/1‑149 (MQ=255)
          tCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGcc                                                                                                                                       <  4:190880/149‑1 (MQ=255)
             cGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTg                                                                                                                                    <  1:386745/149‑1 (MQ=255)
             cGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTg                                                                                                                                    <  6:230499/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                                                                                                            >  5:318133/1‑149 (MQ=255)
                                     ggATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTg                                                                                                            <  8:343709/149‑1 (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                                                                                                                      <  7:116611/132‑1 (MQ=255)
                                            tttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCa                                                                                                                      >  8:116611/1‑132 (MQ=255)
                                             ttGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                      <  7:172247/147‑1 (MQ=255)
                                                   tCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAg                                                                                                                                 >  7:122255/1‑114 (MQ=255)
                                                   tCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAg                                                                                                                                 <  8:122255/114‑1 (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                                                                  >  1:297617/1‑149 (MQ=255)
                                                                               aaaTCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcg                                                                  >  6:427608/1‑149 (MQ=255)
                                                                                  tCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACgtcggtc                                                               <  1:290270/149‑1 (MQ=255)
                                                                                    cGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGtcggtcgg                                                             <  4:204619/149‑1 (MQ=255)
                                                                                      gggTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACagag                                                                             >  1:150724/1‑131 (MQ=255)
                                                                                      gggTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACagag                                                                             <  2:150724/131‑1 (MQ=255)
                                                                                        gTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGccc                                                         >  2:126171/1‑149 (MQ=255)
                                                                                           aGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCggg                                                      >  7:84242/1‑149 (MQ=255)
                                                                                             tCGGCAAGTCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCt                                                    <  8:345419/149‑1 (MQ=255)
                                                                                             tCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                      <  3:146362/99‑1 (MQ=255)
                                                                                             tCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATg                                                                                                      >  4:146362/1‑99 (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                                       >  8:84621/4‑122 (MQ=255)
                                                                                                                                      ggccgGCACCTGCCGCGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTgg                                       <  7:84621/119‑1 (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)
                                                                                                                                               gTGCCGGCCCGGGGCCTTGAAGCAGGATGCACAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGa  >  1:41422/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
GCCTGGCCGATCACGATCGGGCCAGCGCCGAGAATCAGGATGCTTTTGATGTCTGTACGTTTTGGCATGGTGGTCACTCAAATCCGGGGTCAGTCGGCAAGCCGCTTTGAACAATCTGGGTCAGGCGCTGCCGGGCGCGGCAGGTGCCGGCCCGGGGCCTTGAAGCAGGATGCTCAGCGGCGCTTGGCCATGGCATCGGTGAAACGATCGAACAGAGGCGCGACGTCGGTCGGGCCCGGGCTCGCTTCAGGGTGGCCCTGGAAGCTGAACGCGCTCTTGTCGGTACGCTCGA  >  NC_002947/5371860‑5372151

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: