Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,737,921 G→A 5.7% H501H (CAC→CAT hutU ← urocanate hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,737,9210GA 5.7% 897.4 / 26.8 282H501H (CAC→CAThutUurocanate hydratase
Reads supporting (aligned to +/- strand):  ref base G (99/167);  new base A (7/9);  total (106/176)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 285 total aligned reads displayed.

CCAGGCCCTGCTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCT  >  NC_002947/5737773‑5738070
                                                                                                                                                    |                                                                                                                                                     
 cAGGCCCTGCTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGtg                                                                                                                                                     >  2:2257951/1‑150 (MQ=255)
   ggCCCTGCTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGtgct                                                                                                                                                   >  1:866182/1‑150 (MQ=255)
   ggCCCTGCTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGtgct                                                                                                                                                   >  1:3412110/1‑150 (MQ=255)
    gCCCTGCTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGCGCGAAGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGtgctg                                                                                                                                                  >  1:746392/1‑150 (MQ=255)
          cTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGaa                                                                                                                                            <  2:9153483/150‑1 (MQ=255)
          cTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGaa                                                                                                                                            <  2:4500093/150‑1 (MQ=255)
            ccTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGc                                                                                                                                          <  1:6734144/150‑1 (MQ=255)
            ccTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGc                                                                                                                                          <  1:8855815/150‑1 (MQ=255)
              ttGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGccc                                                                                                                                        <  2:3377932/150‑1 (MQ=255)
                 gcgcAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATg                                                                                                                                     <  1:9050110/150‑1 (MQ=255)
                         gATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGcc                                                                                                                             >  1:3969308/1‑150 (MQ=255)
                           tGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCAc                                                                                                                           <  1:5631144/150‑1 (MQ=255)
                                 tGTCATAACCGGCCTCGGCGTAACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGAAGCGCTCGGCGGCCTCAGCGGTACCGTGGCAGACGTTGACCCTGCCCGAGTGCTGCGAGTAGCCGATGCCCACGCCACCGCCAt                                                                                                                     >  2:1595900/1‑150 (MQ=255)
                                  gTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATg                                                                                                                    <  1:2831898/150‑1 (MQ=255)
                                  gTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATg                                                                                                                    >  2:6371657/1‑150 (MQ=255)
                                  gTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATg                                                                                                                    <  2:5924842/150‑1 (MQ=255)
                                   tCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATgg                                                                                                                   <  1:7906470/150‑1 (MQ=255)
                                   tCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTAGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATgg                                                                                                                   >  2:2039034/1‑150 (MQ=255)
                                    cATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCAGCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGt                                                                                                                  >  1:8282390/1‑150 (MQ=255)
                                       aaCCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCa                                                                                                               <  1:6632408/150‑1 (MQ=255)
                                         ccGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGc                                                                                                             <  2:1901998/150‑1 (MQ=255)
                                         ccGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGc                                                                                                             >  2:7557249/1‑150 (MQ=255)
                                          cGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCg                                                                                                            >  1:3711725/1‑150 (MQ=255)
                                               tCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATAcc                                                                                                       >  2:5590759/1‑150 (MQ=255)
                                                cGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATAccc                                                                                                      <  1:1854487/150‑1 (MQ=255)
                                                   cGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAgg                                                                                                   <  1:379825/150‑1 (MQ=255)
                                                    gTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGt                                                                                                  <  1:8086214/150‑1 (MQ=255)
                                                      gACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGtgctgc                                                                                                                                                 <  2:4452364/101‑1 (MQ=255)
                                                      gACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTgg                                                                                                >  1:5783160/1‑150 (MQ=255)
                                                         gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                             >  1:8215390/1‑150 (MQ=255)
                                                         gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                             >  1:5236799/1‑150 (MQ=255)
                                                         gCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGgcgc                                                                                             >  1:6243929/1‑150 (MQ=255)
                                                          cATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGcgcc                                                                                            <  1:3264651/150‑1 (MQ=255)
                                                            tGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGTGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGcgccgc                                                                                          >  2:5283441/1‑150 (MQ=255)
                                                              aCGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCt                                                                                        <  1:5062659/150‑1 (MQ=255)
                                                                 ccAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCg                                                                                     >  2:2866790/1‑150 (MQ=255)
                                                                      cccTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGAGGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGtt                                                                                <  2:8914767/150‑1 (MQ=255)
                                                                       ccTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTc                                                                               <  2:7014118/150‑1 (MQ=255)
                                                                        cTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcc                                                                              >  2:3422476/1‑149 (MQ=255)
                                                                          gggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagc                                                                            <  1:7075967/150‑1 (MQ=255)
                                                                           ggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagca                                                                           >  2:5712107/1‑150 (MQ=255)
                                                                           ggTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTcagca                                                                           <  1:4258295/150‑1 (MQ=255)
                                                                               gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                       >  1:3039777/1‑150 (MQ=255)
                                                                               gTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGc                                                                       >  2:5812313/1‑150 (MQ=255)
                                                                                  ggTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGtt                                                                    <  1:5985719/150‑1 (MQ=255)
                                                                                  ggTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGtt                                                                    <  2:3276678/150‑1 (MQ=255)
                                                                                   gTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTg                                                                   <  1:2564506/150‑1 (MQ=255)
                                                                                   gTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTg                                                                   <  1:4680596/150‑1 (MQ=255)
                                                                                    tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                  <  2:1323604/150‑1 (MQ=255)
                                                                                    tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                  <  2:1689175/150‑1 (MQ=255)
                                                                                    tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                  >  1:7657320/1‑150 (MQ=255)
                                                                                    tCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGa                                                                  <  2:7934016/150‑1 (MQ=255)
                                                                                         aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                             <  2:7328920/150‑1 (MQ=255)
                                                                                         aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                             <  2:382076/150‑1 (MQ=255)
                                                                                         aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                             <  1:6121518/150‑1 (MQ=255)
                                                                                         aCCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGagcagc                                                             <  2:8670348/150‑1 (MQ=255)
                                                                                          cccGGGCGATGCTCTCGGCGGCCTCATCTGTACCGTCGCAGACGATGTCCATGCCCGAGTGCTGCGAGAAGCCCAGGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCCGCGGTGTTCAGCAGGGCGTTGAGCAGCg                                                            <  1:1595900/150‑1 (MQ=255)
                                                                                          cccGGGCGATGCGCTCGGCGGCCTCATCGGTCCCGTCGCAGCCGAGGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCAGGGTGCAGCGATACCCCGGTGGCGCCGCCGGCGGTGTTCAGCAGGGCGTTGAGCAGCg                                                            >  1:5246045/1‑150 (MQ=255)
                                                                                          cccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCg                                                            <  1:6399343/150‑1 (MQ=255)
                                                                                           ccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCgg                                                           <  2:267059/150‑1 (MQ=255)
                                                                                           ccGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCgg                                                           <  2:3972070/150‑1 (MQ=255)
                                                                                                 gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                     <  2:1783195/150‑1 (MQ=255)
                                                                                                 gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                     <  2:4719570/150‑1 (MQ=255)
                                                                                                 gATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                     <  1:3879961/150‑1 (MQ=255)
                                                                                                 gATGCGCTCGGCGGCCTCATCGGGACCGTCGCCGACGATGACCAGGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCCGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTc                                                     <  2:3677024/150‑1 (MQ=255)
                                                                                                   tGCGCTGGTCGGCATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCgg                                                   <  1:2305646/150‑1 (MQ=255)
                                                                                                    gcgcTGGTCGGCATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGa                                                  <  1:4150565/150‑1 (MQ=255)
                                                                                                     cgcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaa                                                 <  2:8305501/150‑1 (MQ=255)
                                                                                                      gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                <  2:2424623/150‑1 (MQ=255)
                                                                                                      gcTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGaaa                                                >  1:7975173/1‑150 (MQ=255)
                                                                                                        tCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACa                                              >  1:4575131/1‑150 (MQ=255)
                                                                                                           tcggcATCTTCCGAACCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCg                                           >  2:2892307/2‑150 (MQ=255)
                                                                                                                   ttccgaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGcc                                   >  1:1125367/7‑150 (MQ=255)
                                                                                                                   aTCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGGGAAGCCCATGCCCACGCCACCGCCATGGTGCAg                                                                                                              <  2:3110702/75‑1 (MQ=255)
                                                                                                                    tCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCa                                  >  1:6596038/1‑150 (MQ=255)
                                                                                                                      ggTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATc                                <  2:8776183/150‑1 (MQ=255)
                                                                                                                       gaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTct                               >  2:9014121/3‑149 (MQ=255)
                                                                                                                       gaaCCGTCGCAGCAGATGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTct                               >  1:8831925/3‑149 (MQ=255)
                                                                                                                            gTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCAt                          <  2:7978908/150‑1 (MQ=255)
                                                                                                                             tCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATg                         >  2:8254738/1‑150 (MQ=255)
                                                                                                                             tCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATg                         <  1:3588422/150‑1 (MQ=255)
                                                                                                                              cGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATgg                        <  2:3765858/150‑1 (MQ=255)
                                                                                                                                cAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCGCCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGcc                      >  2:3531098/1‑150 (MQ=255)
                                                                                                                                  gACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTc                    <  2:1113309/150‑1 (MQ=255)
                                                                                                                                  gACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTc                    >  1:6716450/1‑150 (MQ=255)
                                                                                                                                     gATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGt                 >  1:2500874/1‑150 (MQ=255)
                                                                                                                                      aTGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGtt                <  1:5414446/150‑1 (MQ=255)
                                                                                                                                      aTGACCATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGtt                >  1:4856852/1‑150 (MQ=255)
                                                                                                                                       tGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGttt               >  2:4844578/1‑150 (MQ=255)
                                                                                                                                          ccATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgc            <  2:5070405/150‑1 (MQ=255)
                                                                                                                                          ccATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgc            >  2:1973632/1‑150 (MQ=255)
                                                                                                                                          ccATGCCAGAGCGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGAGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTGTcgc            >  2:3880671/1‑150 (MQ=255)
                                                                                                                                           cATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTcgcg           >  1:8575649/1‑150 (MQ=255)
                                                                                                                                           cATGCCCGAATGTTGGGAGAAGCCCATGCCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGTTTcacg           <  2:5461871/150‑1 (MQ=255)
                                                                                                                                            aTGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCgg          <  2:611073/150‑1 (MQ=255)
                                                                                                                                            aTGCCCGAATGTTGGGAGAAGCCCATACCGACGCCGCCACCATGATGGAGCGATACCCATGTCGCACCCGATGCGGTGTTGAGCAACGCATTGAGCAGCGGCCAGTCGGAGACTGCGTCCGAGCCGTCTTTCATGGATTCGGTTTCACgg          <  2:6760592/150‑1 (MQ=255)
                                                                                                                                               cccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTc       <  1:2733935/150‑1 (MQ=255)
                                                                                                                                                ccGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCg      <  2:2540543/150‑1 (MQ=255)
                                                                                                                                                  gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCggg    <  1:2107970/150‑1 (MQ=255)
                                                                                                                                                  gAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCggg    >  1:1327361/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCAGGCCCTGCTCCTTGGCGCAGTCGATGGCGATGTCATAACCGGCATCGGCGTGACGCATGACGCCAGTCCCTGGGTCGTTGGTCAGTACCCGGGCGATGCGCTCGGCGGCCTCATCGGTACCGTCGCAGACGATGACCATGCCCGAGTGCTGCGAGAAGCCCATGCCCACGCCACCGCCATGGTGCAGCGATACCCAGGTGGCGCCGCCTGCGGTGTTCAGCAGGGCGTTGAGCAGCGGCCAGTCGGAAACAGCGTCCGAACCATCACGCATGGCCTCGGTTTCGCGGTTCGGGCT  >  NC_002947/5737773‑5738070

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: