Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,807,090 G→A 5.0% I923I (ATC→ATT) ‡ 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,0900GA 5.0% 504.1 / 9.8 159I923I (ATC→ATT) ‡PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (84/67);  new base A (3/5);  total (87/72)
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.89e-01

Only 100 of 160 total aligned reads displayed.

GTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGCG  >  NC_002947/4806942‑4807239
                                                                                                                                                    |                                                                                                                                                     
gTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGa                                                                                                                                                      >  2:4305772/1‑150 (MQ=255)
gTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGa                                                                                                                                                      >  1:1049968/1‑150 (MQ=255)
  gCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATg                                                                                                                                                    <  2:2742017/150‑1 (MQ=255)
    aGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGc                                                                                                                                                  <  1:126621/150‑1 (MQ=255)
        aGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATgc                                                                                                                                              <  1:609459/150‑1 (MQ=255)
            gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGa                                                                                                                                          >  2:5027135/1‑150 (MQ=255)
            gATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGa                                                                                                                                          >  2:5027131/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:2595875/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)
                  cGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCAt                                                                                                                                    <  1:2774496/150‑1 (MQ=255)
                      aGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGcc                                                                                                                                >  2:3812630/1‑150 (MQ=255)
                        gTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGa                                                                                                                              >  1:795196/1‑150 (MQ=255)
                         tCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                                             <  2:3544203/150‑1 (MQ=255)
                         tCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                                             <  2:4850474/150‑1 (MQ=255)
                         tCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGAt                                                                                                                             <  1:1907125/150‑1 (MQ=255)
                           cgcgGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATgg                                                                                                                           >  2:2487784/1‑150 (MQ=255)
                            gcgGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGc                                                                                                                          >  1:4274354/1‑150 (MQ=255)
                                gAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcg                                                                                                                      >  1:4401590/1‑150 (MQ=255)
                                  aCCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgac                                                                                                                    >  2:2296843/1‑150 (MQ=255)
                                    cTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCCTTGAGGGCCTCCTTGATACGGGCGGGGGTGAGGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgacgc                                                                                                                  >  2:4577009/1‑149 (MQ=255)
                                    cTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACcgacga                                                                                                                  >  1:2989108/1‑150 (MQ=255)
                                     tCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCgacgac                                                                                                                 >  1:4103047/1‑150 (MQ=255)
                                     tCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCgacgac                                                                                                                 >  1:4103048/1‑150 (MQ=255)
                                         ccGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCggg                                                                                                             <  1:3789304/150‑1 (MQ=255)
                                         ccGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCAGTGCCGATGGCACCGACGACCggg                                                                                                             <  2:3313042/150‑1 (MQ=255)
                                           ggCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGat                                                                                                           <  2:4593298/150‑1 (MQ=255)
                                            gCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGatc                                                                                                          <  1:2054222/150‑1 (MQ=255)
                                               cGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGatcatc                                                                                                       <  2:4751995/150‑1 (MQ=255)
                                                gACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGAtcatca                                                                                                      <  1:1656421/150‑1 (MQ=255)
                                                     gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACt                                                                                                 >  1:86286/1‑150 (MQ=255)
                                                     gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACt                                                                                                 >  2:1678000/1‑150 (MQ=255)
                                                     gTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACt                                                                                                 >  2:3920697/1‑150 (MQ=255)
                                                        gACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCg                                                                                              >  2:932507/1‑150 (MQ=255)
                                                           cggtGAGCGTCAGTTGCAGTGGACGGTTAGCCAGTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCCTTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCg                                                                                           >  1:3243957/4‑150 (MQ=255)
                                                             ttGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGAt                                                                                         <  2:1910699/150‑1 (MQ=255)
                                                                 gCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGccc                                                                                     >  1:1864994/1‑150 (MQ=255)
                                                                 gCGTCAGTTGCAGTGGACGGTTAGCCAATGAGGGCCTCCATGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGccc                                                                                     <  2:1164427/150‑1 (MQ=255)
                                                                  cGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCt                                                                                    <  2:1206235/150‑1 (MQ=255)
                                                                   gTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCtt                                                                                   <  2:3116899/150‑1 (MQ=255)
                                                                      aGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGa                                                                                <  1:69995/150‑1 (MQ=255)
                                                                       gTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAg                                                                               <  2:1259613/150‑1 (MQ=255)
                                                                         tGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAggg                                                                             <  2:3331956/150‑1 (MQ=255)
                                                                          gCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAgggg                                                                            <  1:3796643/150‑1 (MQ=255)
                                                                          gCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAgggg                                                                            >  2:540728/1‑150 (MQ=255)
                                                                              aggACGGTTAGCCATTGAGGGCATCCTTGATAAGGGCTGGGGTGATGGGTAGATCGCGCAAGCGCTAGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGAAAATCACCACTTCGGCGATGCCCTTGGAGGggtcgg                                                                        >  2:450731/2‑150 (MQ=255)
                                                                                   ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                   <  1:4471442/150‑1 (MQ=255)
                                                                                   ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                   <  1:3270646/150‑1 (MQ=255)
                                                                                   ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                   <  1:4726494/150‑1 (MQ=255)
                                                                                   ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                   <  2:4938927/150‑1 (MQ=255)
                                                                                   ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                   <  2:2919754/150‑1 (MQ=255)
                                                                                   ggTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTcggtcggt                                                                   <  1:4471491/150‑1 (MQ=255)
                                                                                       aGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACa                                                               <  1:939561/150‑1 (MQ=255)
                                                                                            ttGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAgcggc                                                          >  1:3559784/1‑150 (MQ=255)
                                                                                                gggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                                      <  1:1598046/150‑1 (MQ=255)
                                                                                                gggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGgga                                                      <  2:2162540/150‑1 (MQ=255)
                                                                                                 ggCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGggag                                                     <  1:3645387/150‑1 (MQ=255)
                                                                                                     tcctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAtt                                                 >  1:98763/1‑150 (MQ=255)
                                                                                                      cctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGAttt                                                <  1:2478947/150‑1 (MQ=255)
                                                                                                       ctTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTc                                               >  1:4088655/1‑150 (MQ=255)
                                                                                                        ttGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCg                                              >  1:1479127/1‑150 (MQ=255)
                                                                                                        ttGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTAGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCg                                              >  1:514757/1‑150 (MQ=255)
                                                                                                           aTACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTg                                           <  1:1989634/150‑1 (MQ=255)
                                                                                                           aTACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTg                                           <  1:3452629/150‑1 (MQ=255)
                                                                                                           aTACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTg                                           <  2:3637422/150‑1 (MQ=255)
                                                                                                            taagggCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACAGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTgg                                          >  1:4578942/4‑150 (MQ=255)
                                                                                                            tACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCAt                                                                                                                                                >  2:2135813/1‑48 (MQ=255)
                                                                                                            tACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCAt                                                                                                                                                <  1:2135813/48‑1 (MQ=255)
                                                                                                             aCGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGt                                         <  2:949106/150‑1 (MQ=255)
                                                                                                               gggCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTct                                       >  1:5037901/1‑150 (MQ=255)
                                                                                                               gggCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCGGGTcg                                       >  1:5037891/1‑149 (MQ=255)
                                                                                                                  cTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgct                                    <  2:2480182/150‑1 (MQ=255)
                                                                                                                   tGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgctg                                   >  2:3550479/1‑150 (MQ=255)
                                                                                                                   tGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTctgctg                                   >  2:798705/1‑150 (MQ=255)
                                                                                                                     gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGcc                                 >  1:1905671/1‑150 (MQ=255)
                                                                                                                     gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGcc                                 >  1:2549280/1‑150 (MQ=255)
                                                                                                                     gggTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGcc                                 >  1:3580126/1‑150 (MQ=255)
                                                                                                                             gggAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGcc                         >  2:1080104/1‑150 (MQ=255)
                                                                                                                              ggAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACCGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCca                        >  2:3780510/1‑150 (MQ=255)
                                                                                                                              ggAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCca                        >  2:1763151/1‑150 (MQ=255)
                                                                                                                               gagaTCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCcac                       >  1:1700677/1‑150 (MQ=255)
                                                                                                                                     gcgcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcg                 >  1:879258/1‑150 (MQ=255)
                                                                                                                                      cgcACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATcgcgc                <  2:1049968/150‑1 (MQ=255)
                                                                                                                                        gtcgcgCATGCCATCGGCAAGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCgcgcgc              >  1:5036034/3‑150 (MQ=255)
                                                                                                                                        gtcgcgCATGCCATCGGCAAGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCgcgcgc              >  1:5035974/3‑150 (MQ=255)
                                                                                                                                           gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg           >  2:146974/1‑150 (MQ=255)
                                                                                                                                           gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg           >  2:147107/1‑150 (MQ=255)
                                                                                                                                           gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg           <  2:215023/150‑1 (MQ=255)
                                                                                                                                           gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg           >  1:2457770/1‑150 (MQ=255)
                                                                                                                                           gcgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATg           >  1:5028085/1‑150 (MQ=255)
                                                                                                                                            cgcTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGc          <  2:2715685/150‑1 (MQ=255)
                                                                                                                                               ttGCCGATGGCATGCGCGACGGCATTGTCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGc       >  1:3172089/1‑150 (MQ=255)
                                                                                                                                               ttGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGc       >  2:4807952/1‑150 (MQ=255)
                                                                                                                                               ttGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGAGCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGc       >  1:3870113/1‑150 (MQ=255)
                                                                                                                                                 gCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGCTTTCGCTGGTCTGCTGCCACACTGCCACATCGCGCGCATGCGGCAg     >  2:2878343/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTGCAGGTAGTCGATCATCGCCAGGTCCGCGGGAACCTCGACCGGCCCGACCGGTTGACCGTTGAGCGTCAGTTGCAGTGGACGGTTAGCCATTGAGGGCCTCCTTGATACGGGCTGGGGTGATGGGGAGATCGCGCACGCGCTTGCCGATGGCATGCGCGACGGCATTGCCGATGGCACCGACGACCGGGATCATCACCACTTCGGCGATGCCCTTGGAGGGGTCGGTCGGTGACAGCGGCGGGAGGATTTCGCTGGTCTGCTGCCACACCGCCACATCGCGCGCATGCGGCAGGCG  >  NC_002947/4806942‑4807239

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: