Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,807,097 T→A 5.2% H921L (CAT→CTT)  PP_4234 ← bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,807,0970TA 5.2% 514.6 / 12.1 153H921L (CAT→CTT) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (80/65);  new base A (3/5);  total (83/70)
Fisher's exact test for biased strand distribution p-value = 4.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 155 total aligned reads displayed.

TAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGCGGTA  >  NC_002947/4806949‑4807242
                                                                                                                                                    |                                                                                                                                                 
tAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATg                                                                                                                                                  <  2:1689859/150‑1 (MQ=255)
   tCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATgcgc                                                                                                                                               >  2:3346429/1‑150 (MQ=255)
     gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGa                                                                                                                                             >  2:5027131/1‑150 (MQ=255)
     gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGa                                                                                                                                             >  2:5114218/1‑150 (MQ=255)
     gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgt                                                                                                                                             <  2:2267291/150‑2 (MQ=255)
     gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgt                                                                                                                                             <  1:432217/150‑2 (MQ=255)
     gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgt                                                                                                                                             <  2:4006404/150‑2 (MQ=255)
     gATCAGCGCCAGGTCCGCGGGAACCTCGCCCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCATGCCATCGGCAAGcgcgt                                                                                                                                             <  1:3481240/150‑2 (MQ=255)
      atcatcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGAc                                                                                                                                            <  2:5020061/150‑1 (MQ=255)
        catcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACgg                                                                                                                                          >  2:712457/1‑150 (MQ=255)
         atcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCTCACGCGCTTGCCGATGGCATGCGCGACGGc                                                                                                                                         <  2:880251/150‑1 (MQ=255)
          tcGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCa                                                                                                                                        <  1:3842171/150‑1 (MQ=255)
               aGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGcc                                                                                                                                   >  2:3812630/1‑150 (MQ=255)
                  tCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                                                <  1:1907125/150‑1 (MQ=255)
                  tCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                                                <  2:4850474/150‑1 (MQ=255)
                             cTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCCTTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgacga                                                                                                                     >  2:2230592/1‑150 (MQ=255)
                             cTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgacga                                                                                                                     >  1:2989108/1‑150 (MQ=255)
                             cTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgacga                                                                                                                     >  1:4478127/1‑150 (MQ=255)
                              tCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCgacgac                                                                                                                    >  1:4103048/1‑150 (MQ=255)
                              tCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCgacgac                                                                                                                    >  1:4727442/1‑150 (MQ=255)
                                  ccGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCggg                                                                                                                <  1:3789304/150‑1 (MQ=255)
                                    ggCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGat                                                                                                              <  2:4593258/150‑1 (MQ=255)
                                    ggCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGat                                                                                                              <  2:4593298/150‑1 (MQ=255)
                                       ccGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGatcat                                                                                                           <  1:2764951/150‑1 (MQ=255)
                                         gACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGAtcatca                                                                                                         <  1:1656421/150‑1 (MQ=255)
                                              gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACt                                                                                                    >  1:86286/1‑150 (MQ=255)
                                              gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACt                                                                                                    >  1:86344/1‑150 (MQ=255)
                                              gTTGACCGTTGAGCGTCAGTGGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACt                                                                                                    >  2:2533838/1‑150 (MQ=255)
                                               ttGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACtt                                                                                                   <  2:2075281/150‑1 (MQ=255)
                                                 gACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCg                                                                                                 >  2:932507/1‑150 (MQ=255)
                                                    cggtGAGCGTCAGTTGCAGTGGACGGTTAGCCAGTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCCTTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCg                                                                                              >  1:3243957/4‑150 (MQ=255)
                                                    cGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCg                                                                                              >  1:3821040/1‑150 (MQ=255)
                                                        gAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGc                                                                                          >  1:3260941/1‑150 (MQ=255)
                                                          gCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGccc                                                                                        >  1:1864994/1‑150 (MQ=255)
                                                           cGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCt                                                                                       <  2:1206235/150‑1 (MQ=255)
                                                           cGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCt                                                                                       <  2:1206205/150‑1 (MQ=255)
                                                            gTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCtt                                                                                      <  2:3116916/150‑1 (MQ=255)
                                                            gTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCtt                                                                                      <  2:3116899/150‑1 (MQ=255)
                                                               aGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGa                                                                                   <  1:69995/150‑1 (MQ=255)
                                                                  tGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAggg                                                                                <  2:3331956/150‑1 (MQ=255)
                                                                   gCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAgggg                                                                               >  2:540728/1‑150 (MQ=255)
                                                                       tGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGggtcgg                                                                           >  1:4093637/1‑150 (MQ=255)
                                                                       aggACGGTTAGCCATTGAGGGCATCCTTGATAAGGGCTGGGGTGATGGGTAGATCGCGCAAGCGCTAGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGAAAATCACCACTTCGGCGATGCCCTTGGAGGggtcgg                                                                           >  2:450731/2‑150 (MQ=255)
                                                                            ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                      <  1:4471442/150‑1 (MQ=255)
                                                                            ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                      <  2:1104480/150‑1 (MQ=255)
                                                                            ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                      <  1:716360/150‑1 (MQ=255)
                                                                            ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                      <  2:3759408/150‑1 (MQ=255)
                                                                            ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                      <  2:2919754/150‑1 (MQ=255)
                                                                                 gCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAg                                                                 >  1:3591500/1‑150 (MQ=255)
                                                                                      tGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAgcggcg                                                            <  2:3821040/150‑1 (MQ=255)
                                                                                         gggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                                         <  1:3066148/150‑1 (MQ=255)
                                                                                         gggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                                         <  1:540728/150‑1 (MQ=255)
                                                                                         gggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                                         <  2:2162364/150‑1 (MQ=255)
                                                                                         gggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                                         <  1:4326408/150‑1 (MQ=255)
                                                                                         gggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                                         <  1:4016246/150‑1 (MQ=255)
                                                                                             ctcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAt                                                     >  1:36830/1‑150 (MQ=255)
                                                                                              tcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAtt                                                    >  1:98763/1‑150 (MQ=255)
                                                                                              tcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAtt                                                    >  1:98863/1‑150 (MQ=255)
                                                                                              tcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAtt                                                    >  1:4199133/1‑150 (MQ=255)
                                                                                              tcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGAAGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAtt                                                    >  2:1316882/1‑150 (MQ=255)
                                                                                               cctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAttt                                                   <  1:2478947/150‑1 (MQ=255)
                                                                                                ctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTc                                                  >  2:2332893/1‑150 (MQ=255)
                                                                                                ctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTc                                                  >  1:1635724/1‑150 (MQ=255)
                                                                                                 ttGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCg                                                 >  1:1479127/1‑150 (MQ=255)
                                                                                                 ttGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTAGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCg                                                 >  1:514757/1‑150 (MQ=255)
                                                                                                  tGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGc                                                <  1:4434039/150‑1 (MQ=255)
                                                                                                    aTACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTg                                              <  2:3637422/150‑1 (MQ=255)
                                                                                                    aTACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTg                                              <  1:1989634/150‑1 (MQ=255)
                                                                                                     taagggCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACAGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTgg                                             >  1:4578942/4‑150 (MQ=255)
                                                                                                     tACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCAt                                                                                                                                                   <  1:2135813/48‑1 (MQ=255)
                                                                                                      aCGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGt                                            <  2:949106/150‑1 (MQ=255)
                                                                                                        gggCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTct                                          >  1:5037880/1‑150 (MQ=255)
                                                                                                        gggCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCGGGTcg                                          >  1:5037891/1‑149 (MQ=255)
                                                                                                           cTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgct                                       <  2:2480182/150‑1 (MQ=255)
                                                                                                            tGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgctg                                      >  2:798705/1‑150 (MQ=255)
                                                                                                            tGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgctg                                      >  2:3550479/1‑150 (MQ=255)
                                                                                                              gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGcc                                    >  1:3580126/1‑150 (MQ=255)
                                                                                                              gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGcc                                    >  1:1905671/1‑150 (MQ=255)
                                                                                                              gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGcc                                    >  1:3303578/1‑150 (MQ=255)
                                                                                                              gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGcc                                    >  1:2549280/1‑150 (MQ=255)
                                                                                                                      gggAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGcc                            >  2:1080104/1‑150 (MQ=255)
                                                                                                                      gggAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGCCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGcc                            <  1:910362/150‑1 (MQ=255)
                                                                                                                       ggAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCca                           >  2:1763151/1‑150 (MQ=255)
                                                                                                                        gagaTCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCcac                          >  1:1700677/1‑150 (MQ=255)
                                                                                                                                 gtcgcgCATGCCATCGGCAAGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCgcgcgc                 >  1:5035974/3‑150 (MQ=255)
                                                                                                                                 gtcgcgCATGCCATCGGCAAGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCgcgcgc                 >  1:5035968/3‑150 (MQ=255)
                                                                                                                                  aCGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCa                >  2:1660878/1‑150 (MQ=255)
                                                                                                                                    gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg              >  2:146974/1‑150 (MQ=255)
                                                                                                                                    gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg              >  1:2457770/1‑150 (MQ=255)
                                                                                                                                    gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg              <  2:215023/150‑1 (MQ=255)
                                                                                                                                     cgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGc             <  2:2715685/150‑1 (MQ=255)
                                                                                                                                        ttGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGc          >  2:4807952/1‑150 (MQ=255)
                                                                                                                                        ttGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGc          >  2:4807889/1‑150 (MQ=255)
                                                                                                                                        ttGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGc          >  1:1305988/1‑150 (MQ=255)
                                                                                                                                          gCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGCTTTCGCTGGTCTGCTGCCACACTGCCACATCGCGCGCATGCGGCAg        >  2:2878343/1‑150 (MQ=255)
                                                                                                                                          gCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAg        >  2:2684826/1‑150 (MQ=255)
                                                                                                                                            cGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGc      >  1:2012175/1‑150 (MQ=255)
                                                                                                                                             gATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGCg     >  1:3511566/1‑150 (MQ=255)
                                                                                                                                                ggCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGCGGTa  <  1:3920697/150‑1 (MQ=255)
                                                                                                                                                ggCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGCGGTa  <  2:2012175/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGCGGTA  >  NC_002947/4806949‑4807242

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: